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PREFACE 


Beginning  in  late  1978,  the  United  States  Air  Force  through  the  Office  of 
Aerospace  Research  initiated  a  major  science  of  fracture  project  at  the  Univer¬ 
sity  of  Pittsburgh.  The  rationale  for  this  special  effort  emerged  from  the  re¬ 
cognition  that,  in  some  form,  fracture  analysis  was  becoming  more  and  more  gene¬ 
rally  employed  as  an  engineering  design  tool  to  assess  structural  integrity. 

Fracture  technology  was  increasingly  sophisticated  and  complex.  With  the  many 
static,  fatigue,  and  spallation  problems  being  attacked  using  different  analytical, 
numerical,  and  experimental  approaches,  and  the  many  theories  -  Griffith,  Wiebull, 
maximum  strain,  and  hydrodynamic  -  appropriate  for  metals,  plastics,  elastomers, 
ceramics,  and  composites,  it  seemed  timely  to  inquire;  is  there  sufficient  experi¬ 
ence  to  generalize  the  state-of-the-art,  i.e.  to  provide  the  basis  for  a  comp¬ 
rehensive  fracture  analysis  discipline  which  is  fundamentally  consistent  and 
complete? 

The  University  of  Pittsburgh  program  was  designed  to  answer  this  question. 

The  effort  was  organized  in  terms  of  five  program  categories  each  containing  diffi¬ 
cult  but  fundamental  problems  whose  solution  would  significantly  contribute  to  the 
desired  unification.  The  five  areas,  to  be  construed  in  a  broad  sense,  were 

1.  Analytical  efforts  and  basic  mechanism 

2 .  Structure-property  relat ionships 

3.  Numerical  analysis  methods 

4.  Instrumentation  science 

5.  Technology  transfer  and  applications 

Primary  emphasis  during  the  early  phases  was  upon  the  first  two  categories,  with, 
simultaneously,  a  lower  level  of  effort  devoted  to  the  latter  three. 

In  a  research  sense,  the  program  is  high-risk  and  for  this  reason  was  ini¬ 
tially  proposed  for  a  five-year  period.  Subsequently  as  a  result  of  an  adminis¬ 
trative  review  by  the  Office  of  Scientific  Research  in  May  1978,  the  technical 
effort  was  reduced  and  redirected  to  emphasize  specific  tasks  structured  in  the 
more  traditional  manner.  Consequently  this  report,  covering  the  period  15  August 
1979  to  30  September  1980,  reflects  a  transition  to  more  concentrated  attention 
upon  the  first  two  of  the  above  five  areas. 


M.L.  Williams 
C.C.  Yates 

October  1980 

Pittsburgh,  Pennsylvania 
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1.0  INTRODUCTION 


The  initial  thrust  of  this  program  was  directed  toward  the  analysis  of 
continuum  fracture  and  the  structural  integrity  as  influenced  by  materials 
science,  mechanics,  and  the  loading  environment.  Five  components  of  the  pro¬ 
gram  were  identified. 

1.  Analytical  Efforts  and  Basic  Mechanisms 

2.  Structure-Property  Relationships 

3.  Numerical  Analysis  Methods 

A.  Instrumentation  Science 

5.  Technology  Transfer  and  Applications 
As  a  result  of  a  USAF  Office  of  Scientific  Research  administrative  re¬ 
view,  the  technical  work  will  be  reduced  and  redirected  to  emphasize  tasks 
structured  in  a  more  traditional  task  oriented  manner.  Consequently,  during 
this  current  transition  year,  the  work  will  reflect  a  concentration  primarily 
upon  the  first  two  of  the  original  five  component  parts . 

2.0  STATUS  OF  RESEARCH 

2.1  Analytical  Efforts  and  Basic  Mechanisms 

The  keystone  for  understanding  fracture  rests  upon  knowledge  of  mech¬ 
anisms  which  control  the  process  and  upon  the  accompanying  analytical  work. 
This  approach  generally  proves  more  cost-effective  than  an  unguided  testing 
program  and  engineering  history.  The  method  chosen  is  to  study  a  class  of 
idealized  problems  which  are  analytically  tractable  yet  reveal  key  elements 
and  mechanisms  of  fracture  behavior. 

Four  illustrations  were  given  in  the  original  proposal.  They  were  (i) 
the  three-dimensional  stress  singularity,  (ii)  the  elemental  "toothpick"  pro¬ 
blem,  (iii)  resolution  and  comparison  of  fracture  criterion  for  metals,  poly¬ 
mers,  and  ceramics,  and  (lv)  crack  growth  and  environmental  effects. 


Nearly  all  of  modern  fracture  analysis  Is  based  upon  the  characteristic 
inverse  square  root  elastic  stress  singularity  which  was  deduced  for  plane 
stress  (thin  plate  stock)  or  plane  strain  (infinitely  thick  plates).  Prac¬ 
tically  speaking,  all  plates  have  a  finite  thickness  and  when  they  fracture 
they  usually  exhibit  important  features  which  cannot  be  predicted  using  two- 
dimensional  theory.  Consequently,  with  success  obtained  using  the  simple 
plane  stress  result,  it  is  appropriate  to  encourage  continuing  work  on  the 
still  unsolved  3-D  stress  distribution.  Our  internal  effort  has  been  devoted, 
as  planned,  to  maintaining  contact  with  others  attempting  to  solve  this  pro¬ 
blem,  particularly  Dr.  E.S.  Folias  (University  of  Utah),  and  Dr.  T.  Kawai 
(Institute  of  Industrial  Science,  University  of  Tokyo).  In  addition.  Dr. 
Atkinson  who  has  followed  this  work  has  prepared  a  conmentary  included  here¬ 
in  as  Appendix  4.1.  At  best  the  present  situation  can  be  described  as  unset¬ 
tled;  there  is  no  present  analytical  solution  which  can  be  used  confidently. 
Dr.  Atkinson  has  also  contributed  a  position  paper  (Appendix  4.2)  dealing  with 
the  use  of  strain  energy  release  rates  in  connection  with  the  mathematically 
important  problem  of  discontinuous  elastic  stress  singularities  which  occur, 
for  example,  when  a  crack  suddenly  appears  at  a  90  degree  corner.  The  order 
of  singularity  in  this  latter  case  jumps  discontinuously  from  its  uncracked 
90  degree  geometry  to  the  stronger  -  1/2  value  characteristic  of  the  zero  de¬ 
gree,  free-sided  crack  geometry.  The  present  conclusion  is  that  there  is  no 
viable  answer  to  the  problem:  the  transition  from  no-crack  to  crack  remains 
fundamentally  discontinuous.  The  analogy  to  a  circular  hole  in  a  plane  stress 
specimen  can  be  drawn.  The  stress  at  a  circular  hole,  regardless  how 
small  -  but  finite  -  it  may  be,  varies  as  the  inverse  square  of  the  distance 
from  the  flaw  radius.  If,  however,  there  is  no  hole  in  the  plate  in  the 
first  place,  the  stress  is  bounded  by  a  constant.  Hence,  there  is  a  jump  in 
singularity  with  geometry.  Appendix  4.3,  also  by  Atkinson,  summarizes  some 


recent  work  on  non-local  elastic  media  and  its  application  to  the  line  crack 
problem.  He  concludes  that  recent  claims  that  such  theories  lead  to  a  finite 
stress  at  a  crack  tip  are  unjustified. 

Atkinson,  Smelser  and  Sanchez  (Appendix  4.4)  have  combined  to  report 
continuing  work  utilizing  the  cracked,  diametrically  loaded  Brazil  disk  to 
deduce  the  Mode  II  fracture  toughness.  This  shear  mode,  which  is  difficult 
to  generate  experimentally,  is  characteristic  of  a  pure  shear  failure  as  might 
occur  in  a  spinning  disk  or  turbine  rotor.  One  of  the  virtues  of  this  Brazil 
test  is  the  possibility  of  assessing  mixed  mode  fracture  toughness,  perhaps  in 
a  better  way  than  using  angled  cracks  in  a  plane  stress  tensile  specimen. 

In  a  manuscript  being  prepared  for  journal  publication,  the  Avila  and 
Betz  work  on  the  adhesive  debonding  of  circular  rods  from  the  matrix  will  be 
reported.  This  class  of  free-surface  debond  geometries,  as  well  as  encapsulated 
(fiber)  rods  are  important  in  composite  materials.  The  primary  characteristic 
controlling  the  debonding  is  believed  to  be  the  specific  adhesive  fracture 
energy  along  the  interface.  One  way  of  measuring  this  interfacial  energy 
to  engineering  accuracy  is  to  use  a  variation  of  the  blister  test.  Somayaji 
(Appendix  4.5)  has  assembled  the  requisite  experimental  assembly  for  a  plexi¬ 
glass-steel  combination  using  a  point  loading  for  these  relatively  stiff  materi¬ 
als.  A  pressure  loading  is  also  contemplated  for  possible  use  with  a  glass- 
polymer  or  glass-epoxy  bimaterial  combination.  Such  an  independent  measure¬ 
ment  of  the  adhesive  energy  permits  one  to  then  use  its  (now  known)  value  in 
analytical  work  using  the  same  bimaterial  combination,  e.g.  fiberglass  rein¬ 
forced  epoxy  resins.  Our  experiments  would  hopefully  use  the  "standard  materi¬ 
als"  being  developed  by  Kelley  at  the  University  of  Akron.  To  prepare  for  such 
analytical  work  a  literature  survey  of  various  contributions  to  the  behavior. 
Including  fracture,  of  composite  materials  has  been  assembled  for  later  use 
(Appendix  4.6). 
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Fracture  criteria  for  various  materials,  preferably  in  some  universal 

form,  is  being  investigated  analytically  and  experimentally  (see  University  of 

Akron  sub-contract).  In  principle,  the  idea  may  be  understood  by  considering 

2 

the  usual  Griffith  fracture  criterion  expressed  in  the  form  a  a=(cEy,  in  which 

K  is  a  numerical  constant  depending  upon  geometry,  E  is  the  modulus  and  Y  the 

2 

fracture  energy.  If  it  is  written  in  the  form  a  /(2E)=(</2) (y/a) ,  and  the 

2 

(uniaxial)  strain  energy  density,  W,  is  recognized  as  W =0  /(2E)  which  at  fail¬ 
ure,  say  Wf,  is  the  area  under  the  brittle  stress-strain  curve  to  failure, 
then  one  has  y=K'aWf,  with  k'  a  numerical  constant  and  "a"  the  crack  size. 

One  might  then  jump  to  a  hypothesis  that  in  any  material,  the  fracture  energy 
for  a  given  geometry  and  crack  size  (obtained  say  from  NDE)  may  be  deduced 
from  the  area  under  the  uniaxial  tensile  stress-strain  curve  to  failure.  Kelley 
at  the  University  of  Akron  has  used  essentially  this  result,  as  proposed  by 
Andrews  for  linear  viscoelastic  polymers.  The  strength  of  materials  approach 
recognized  this  principle  many  years  ago  in  denoting  the  area  under  any  material 
stress-strain  curve  to  failure  as  "toughness".  It  is  hoped  that  the  results 
of  the  numerical  analysis  being  conducted  by  Becker  (Appendix  5.2)  can  be  used 
in  checking  a  general  fracture  criterion  through  a  knowledge  of  the  actual 
large  strain  at  a  crack  tip. 

Another  comparative  fracture  analysis  has  been  pursued  by  Batdorf 
(Appendix  4.7)  during  his  one  month  visit  with  us.  His  contribution  is  a 
position  paper  dealing  with  the  micro-mechanics  of  failure  as  applied  to  mixed 
media,  in  particular  to  fibrous  composites.  The  eventual  aim,  which  may  not 
be  achieved  during  the  next  year,  is  to  see  if  there  is  any  common  ground  or 
criteria  for  which  both  the  energy  balance  (Griffith)  criterion  and  statistical 
(Weibull)  criterion,  applicable  principally  to  ceramics,  can  be  combined  or 
compared. 


2.2  Structure  -  Property  Relationships 


The  relation  between  the  chemical  structure  of  a  material  and  its  mech¬ 
anical  behavior  is  perhaps  the  most  fundamental  problem  of  all.  It  affects 
not  only  the  constitutive  law,  but  controls  processing  characteristics  and 
the  failure  mode. 

We  contemplated  the  development  of  an  Interaction  Matrix  which  relates 
parameters  describing  the  material  structure  (Burger's  vector,  cross-link 
density,  Weibull  modulus,  etc.)  with  those  describing  its  mechanical  proper¬ 
ties.  In  the  subcontract  arrangement  with  Dr.  F.N.  Kelley  and  Dr.  P.  Dreyfuss, 
University  of  Akron,  the  intent  is  to  synthesize  a  "standard  elastomer"  and 
a  "standard  epoxy"  having  controlled  and  known  chemical  compositions.  The  cur¬ 
rent  status  of  the  work  is  given  in  Appendix  5.1.  During  this  first  year  the 
principal  effort  has  been  concentrated  upon  designing  the  experimental  setups 
and  obtaining  reproducible  synthesized  material  from  systems  which  are  chemical¬ 
ly  clean  and  as  environmentally  insensitive  as  practicable.  (Subsequent  doping 
for  environmental  effects  was  originally  contemplated,  but  this  effort  will  not 
now  be  conducted  on  this  program.)  Two  systems  were  proposed:  the  polymeric 
because  there  is  a  already  linearly  viscoelastic  analysis  capability  in  the  com¬ 
munity  which  can  therefore  be  used  to  relate  theory  and  experiment,  and  the  epoxy 
resin  because  as  a  usually  non-linear  material  the  analysis  capability  is  very 
limited  and  definitive  material  property  data  is  needed  to  guide  the  theoretical 
assumptions.  So,  in  both  cases  material  characterization  Is  needed  -  in  the 
first  to  verify  the  analysis  techniques,  and  in  the  second  to  guide  the  stress- 
strain  assumptions  incorporated  into  possible  analyses.  As  the  University  of 
Akron,  (Appendix  5.1)  reinforced  by  internal  support  from  Plazek  (Appendix  4.8), 
preliminary  machanical  and  thermal  characterization  of  rubbery  and  glassy  network 
polymers  has  been  completed  with  an  empahsis  on  fracture  and  tear  energies  as  a 
function  of  temperature  reduced  strain  rate.  Current  chemical  structure  variables 
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are  controlled  cross-link  density  and  chain  back-bone  flexibility.  In  the  model 
epoxy  system,  a  method  has  been  found  to  describe  local  homogeneities  by  plac¬ 
ing  accurately  known  spherical  particles  in  the  polymer  matrix. 

One  of  the  important  factors  affecting  the  mechanical  properties  of  a 
material  is  the  effect  of  the  fabrication  process  and  precursor  material  varia¬ 
bility.  Dr.  J.T.  Lindt  has  completed  a  preliminary  investigation  of  "melt  frac¬ 
ture"  (Appendix  4.9).  This  work  will  not  be  continued,  however,  due  to  funding 
limitations. 

Two  other  topics  complete  our  present  activities  in  structure-property 
relationships.  The  first,  by  Dr.  H.  Ruhn  (Appendix  4.10)  is  related  somewhat 
to  Dr.  J.T.  Lindt's  polymer  processing  experiment.  Dr.  Kuhn  is  investigating 
the  strain  threshold  at  fracture  during  forming,  especially  in  compression,  such 
as  during  a  forging  operation.  He  has  also  made  micro-analyses  of  a  cobalt- 
nickel  alloy,  known  as  MP  35,  and  found  that  void  nucleation  is  the  mechanism 
which  controls  the  first  stage  of  the  fracture  process.  It  occurs  by  grain 
boundary  sliding  induced  by  dynamic  recrystallization.  Finally,  the  only  lab¬ 
oratory  work  dealing  with  ceramics  has  been  conduced  by  Porter  (Appendix  4.11). 
The  work  he  has  completed  to  date  has  concentrated  upon  obtaining  consistent 
powders  and  evaluating  the  reproducibility  of  simple  shaped  fabricated  speci¬ 
mens.  The  only  processing  variable  considered  so  far  is  surface  contamination 
of  silicon  carbide  powder  particles.  Low  temperature  plasma  heating  during 
the  sintering  was  shown  to  produce  apparently  clean,  active  surfaces  in  the  pow¬ 
der  compacts.  Continuing  tests  will  probably  be  directed  toward  varying  the 
initial  green  densities  of  the  powder. 
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3.0  Statement  of  Work 


The  following  research  is  to  be  conducted: 

a)  Investigate  a  class  of  idealized  problems  which  are  analytically 
tractable  but  key  to  the  understanding  of  fracture  mechanisms  and 
their  control.  This  class  of  problems  will  include,  but  not  be 
limited  to,  the  following: 

1.  Study  the  three  dimensional  stress  singularities  along  a  crack 
surface  and  its  interface  with  free  surfaces  or  other  material 
interfaces,  as  well  as  the  engineering  importance  of  these 
singularities.  Investigate  the  nature  of  this  stress  sing¬ 
ularity  in  terms  of  typical  crack  geometric  parameters  and  the 
nature  of  the  mathematical  discontinuity  as  these  parameters 
approach  zero . 

2.  Investigate  the  fracture  behavior  of  an  elastic  fiber  in  a  matrix 
material  considering  typical  interfacial  defects  and  uniaxial 
and  combined  stress  states  on  the  fiber.  Investigate  the  adhe¬ 
sive  and  cohesive  strength  mechanisms  of  the  fiber  matrix  inter¬ 
face  , 

3.  Investigate  the  extension  of  models  of  fracture  mechanisms  of 
the  micro/macro  scale  in  order  to  predict  fatigue  crack  growth, 
including  as  appropriate,  strain  rate  effects,  temperature, 
environmental  conditions,  and  nonlinear  behavior. 

4.  For  the  classes  of  problems  to  be  addressed  under  "a"  -  metals, 
polymers  and  ceramics,  individually  and  in  combination,  will  be 
considered  in  the  development  of  fracture  criteria  in  which 
particular  significance  can  be  attributed  to  constitutive  equa¬ 
tions,  thermodynamics,  kinetics  or  material  processing, 

b)  Investigate  the  structure-property  characteristics  of  metals,  poly- 
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mers  and  ceramics  with  respect  to  both  adhesive  and  cohesive  fracture 
behavior. 

1.  Review  the  existing  literature  and  develop  an  interaction  matrix 
between  the  chemical  and  crystalline  or  noncrystalline  structure 
of  a  material  and  its  mechanical  properties. 

2.  Investigate,  both  theoretically  and  experimentally,  the  quanti¬ 
tative  nature  of  this  interaction  and  its  significance  to  frac¬ 
ture. 

3.  Coordinate  with  AFOSR  the  planned  investigation  to  develop  the 
quantitative  relationship  of  the  interaction  matrices. 

c )  Review  the  current  state-of-the-art  in  numerical  analysis  methods 
which  are  applicable  to  the  science  of  fracture,  considering,  in 
particular,  the  micro/macro  modeling  levels,  kinetic  phenomena,  and 
thermodynamic  considerations. 

d)  Investigate  the  state-of-the-art  in  the  field  of  instrumentation  and 
measurement  science  which  has  application  to  the  science  of  fracture 
and  structure/property  relationships.  In  particular,  include  investi¬ 
gation  of  techniques  which  may  be  developed  to  study  in  a  nondestruc¬ 
tive  manner  internal  material,  i.e.,  nonsurface,  behavior. 

e)  Develop  an  active  program  to  transfer  fracture  research  results  to 
the  user  community  and  to  maintain  as  assessment  of  critical  fracture 
technology  requirements  on  a  national  level. 

Although  some  work  was  continued  during  the  second  year  in  the  last  three 
areas,  i.e.  numerical  analysis,  instrumentation  science  and  technology  transfer, 
it  was  at  a  considerably  reduced  level  in  view  of  the  changed  emphasis  in  the 
overall  program. 
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V  V 


1.  Introduction 


The  way  In  which  the  stress  field  varies  at  notches  and  angular 
corners  in  plates  is  obviously  of  importance  to  failure  analysis.  In 
particular  the  square  root  stress  singularity  at  the  tip  of  a  line  crack 
in  plane  stress/plane  strain  plays  a  central  role  in  fracture  mechanics. 

An  enormous  amount  of  work  has  been  done  in  this  subject,  much  of  the 
recent  work  being  reviewed  in  (1). 

As  a  representative  of  early  work  note  the  work  of  Knein  (2,3) 
who  considered  an  elastic  wedge  of  angle  n/2  in  plane  strain,  one  side 
having  zero  tractions,  the  other  having  zero  displacement.  Assuming 
displacement  components  of  the  form 

u±  -  rA  fi(X,  6),  i  -  1,  2  (1.1) 

with  (r,  6)  polar  coordinates,  then  the  eigenvalues  X  are  solutions  of  the 
equation 

(Sin(irX/2)  }2  -  (4(1  -  v)2  -  X2}/(3  -  4v)  (1.2) 

Later  Williams  (4,5)  considered  various  boundary  conditions  in  angular 
corners  for  plates  under  bending  or  extension.  For  a  wedge  of  angle  a  in 
plane  strain  with  traction  free  sides  he  looked  for  separable  solutions  of 
type  (1.1)  and  deduced  an  eigenvalue  equation  for  X,  i.e., 

(Sin(Xa))2  -  X2(Sin  a)2  -  0  (1.3) 

Many  other  results  of  this  type  have  been  found,  some  of  which  are  referred 
to  in  (1) . 

The  method  used  in  (3)  and  (4)  is  to  assume  a  solution  separable  in 
polar  coordinates  and  then  to  deduce  the  eigenvalue  equation  for  X  such 
that  zero  boundary  conditions  are  satisfied  on  the  infinite  wedge  sides. 


13 


The  usefulness  of  such  expressions  Is  that  close  to  the  corner  the  terns 
with  A  small  should  dominate,  the  associated  stress  field  often  being 
singular  (i.e.,  leading  to  infinite  stresses  as  r  -*•  0) .  Furthermore  the 
smallest  A  allowed  by  considerations  of  finite  energy  may  be  complex. 

This  happens,  for  example,  in  the  case  of  a  crack  at  an  interface  between 
two  dissimilar  elastic  media  {(6),  see  also  (1)}. 

For  plane  elastic  systems  the  above-mentioned  results  are  relatively 
simple  to  deduce,  although  some  care  is  still  required.  For  example  it 
may  not  be  sufficient  to  assume  a  solution  of  the  form  u  ■  r*  f(9),  the 
r  dependence  may  involve  log  r  terms.  To  illustrate  this  we  give  in 
Appendix  A  a  solution  to  a  plane  potential  problem  (anti-plane  elasticity) 
where  this  is  easily  seen.  The  approach  we  use  in  Appendix  A  is  to  employ 
the  Mellln  transform  and  in  this  way  the  eigenvalue  results  appear  quite 
naturally  as  the  poles  of  the  resulting  transform.  Furthermore,  arguments 
based  on  the  Mellin  transform  are  potentially  useful  in  discussing  the 
question  of  completeness  of  the  eigenfunctions  for  the  wedge.  It  is  worth 
enlarging  on  the  completeness  question  here  since  the  reader  will  see  it 
being  mentioned  in  the  context  of  the  '3-d  singularity'  discussed  later. 

Given  the  set  of  eigenfunctions  (1.1)  and  their  associated  eigenvalues.  A, 
a  natural  question  to  ask  is  whether  they  are  complete  in  the  sense  that 
they  can  represent  fairly  general  boundary  data  on  the  arc  r  ■  1  say.  It 
is  not  sufficient  to  show  that  infinite  sums  of  such  eigenfunctions  converge 
on  r  «  1,  it  is  also  necessary  to  show  that  they  converge  to  the  data. 
Considerations  of  this  type  have  been  made  most  recently  by  Gregory  (7)  for 
plane  elastic  problems.  In  particular  he  shows  that  William's  eigenfunctions 

4 

for  V  $  •  0  in  the  wedge  are  complete  in  the  sense  that  they  may  be  used  to 


expand  the  data 


where  £(8),  g(6)  are  'arbitrary'  independent  functions  prescribed  on 
(-a^,  a^) ,  if  2a^  is  the  wedge  angle  (a  of  (1.3)  corresponds  to  2a^  here). 

See  (7)  and  (8)  for  more  details  of  these  results.  Thus  stress  singular¬ 
ities  in  the  two  dimensional  elastic  situation  are  reasonably  well  under¬ 
stood.  However  it  should  be  noted  that  the  nature  of  the  singularity 
(i.e.,  the  exponent  X)  and  its  angular  variation  (f(X,6)}  although 
essential  are  only  the  first  step  for  a  realistic  stress  analysis.  The 
coefficient  of  the  eigenfunction  which  will  result  from  the  influence  of 
the  applied  load  is  equally  important,  and  the  determination  of  this 
coefficient  (e.g.,  the  stress  Intensity  factor  when  a  -  2ir,  ■  ir)  may 
involve  extensive  numerical  analysis. 

The  situation  regarding  three  dimensional  stress-singularities  where 
cracks  intersect  free  surfaces,  etc.,  is  much  more  difficult  to  analyze  and 
the  nature  of  the  singularity  has  'perhaps'  only  recently  been  illuminated. 

We  use  the  adverb  'perhaps'  because  there  is  not  yet  total  agreement  on  the 
results  and  also  because  naturally  there  is  not  yet  complete  mathematical 
rigor  associated  with  the  calculations .  Nevertheless  a  plausible  and 
fairly  comprehensive  investigation  of  the  singularity  at  the  quarter-infinite 
crack  in  a  half-space  has  been  made  by  Benthem  (9),  (10),  and  (11).  His 
work  has  the  additional  virtue  that  he  has  repeated  his  earlier  calculations 
by  a  different  method  in  his  later  paper.  Furthermore,  there  is  some 
measure  of  agreement  between  Benthem's  results  and  results  obtained  by 


BazanC  and  Estenssoro  (12)  who  used  a  finite  element  method  of  determining 
critical  eigenvalues.  Each  of  the  above-mentioned  authors  (9-12)  use 
spherical  polar  coordinates  (r,  6,  <p)  with  origin  0  at  the  terminal  point 
at  which  the  crack  front  edge  meets  the  free  surface  (see  Fig.  1) . 

They  then  assume  the  solution  to  be  separable  in  the  form  r^  F(A,  0,  $) 
and  look  for  the  eigenvalues  A  and  associated  eigenfunctions  F.  The  methods 
used  for  doing  this  will  be  discussed  later. 

Before  discussing  the  results  of  the  above  work  in  more  detail  it  is 
appropriate  to  indicate  other  work  which  has  appeared  in  the  literature. 

In  (12)  these  other  works  are  referred  to  by  the  sentence:  -  "Other 
analytical  solutions  for  orthogonal  cracks  (13,14)  have  also  been  presented, 
but  they  disagree  with  the  results  of  Benthem  (10)  as  well  as  the  present 
paper".  We  have  changed  the  reference  number  to  agree  with  the  numbering 
used  here.  In  (11),  Benthem  referring  to  his  papers  (9)  and  (10)  states 
"The  need  was  felt  to  confirm  these  results  independently  by  quite  another 
method.  In  some  private  correspondences  and  concept-papers  kindly  sent  to 
the  author  the  results  of  (9,10)  were  questioned".  He  does  not  refer  to 
any  of  the  opposing  work  specifically.  Note,  however,  that  Benthem  and 
Koiter  (15)  discuss  the  paper  by  Folias  (14)  and  in  his  reply  Folias  raises 
questions  concerning  Benthem's  report  (9).  One  of  these  questions,  as  to 
whether  the  problem  admits  a  solution  in  separable  form ^particularly  in 
6  and  <t>,  is  avoided  by  the  approach  used  in  Benthem’s  later  paper  (11). 

No  reference  at  all  to  (13)  and  (14)  is  made  by  Benthem  in  (10)  and  (11) 
and  only  the  sentence  quoted  above  made  by  Bazant  and  Estenssoro  (12). 

This  is  unfortunate  since  constructive  criticism  would  have  been  useful  in 
forming  a  balanced  view  (such  criticism  may  have  occurred,  of  course,  in 
references  we  have  overlooked) .  There  is  available  however  a  reported 


discussion  by  Kawai  and  Folias  which  is  useful  insofar  as  it  gives  another 
view  of  the  problem  and  highlights  possible  criticisms  of  the  other  work. 

The  report  of  this  discussion  which  was  prepared  by  Dr.  R.E.  Smelser  is 
included  as  Appendix  2.  It  is  clear  from  this  discussion  that  these  authors 
are  not  entirely  satisfied  with  the  approach  of  (9)  and  (10)  but  are  not 
really  ready  to  offer  a  viable  alternative. 

In  this  paper  we  shall  adopt  the  position  that  the  results  of 
Benthem  (11)  are  the  best  available  at  the  moment,  while  leaving  our 
options  open  just  in  case  more  exotic  behavior  at  the  corner  point  is 
discovered  in  the  future.  Even  if  the  singularity  at  the  corner  point 
were  more  exotic,  involving  log  r  terms  say,  the  exponents  calculated  by 
Benthem  would  still  probably  be  a  good  approximation  since  log  r  tends 
to  infinity  very  slowly  as  r  tends  to  zero  making  r*  log  r  and  r^  hardly 
discernible  as  r  -*■  0. 
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2.  Benthem's  Results  for  the  Quarter-Infinite  Crack  in  a  Half-Space  (Fig.  1) 
The  crack  geometry  is  shown  in  Figure  1.  Cartesian  coordinates  (x,  y,  z) 


and  spherical  coordinates  (r,  6,  <p)  are  used  where 

x  *  r  Sin  9  Cos  <(> 

y  *  r  Sin  9  Sin  $  (2.1) 

z  *  r  Cos  9 

The  half-space  is  the  region  z  >  0  (free-surface  z  •  0)  and  the  crack  occupies 
the  quarter-plane  y  =  0,  x>0,  z>0.  Thus  the  region  to  be  investigated  is 
a  conical  region  (with  vertex  at  r  =  0)  between  the  three  planes  9  **  ir/2, 

<j>  *  40  and  $  *  2ir  -  0.  These  planes  are  taken  to  be  stress  free.  In  this 
conical  region  Benthem  considers  a  state  of  stress  which  in  terms  of  the 
cartesian  stress-components  is  written 

°xx  *  r*  fxx(A*  9’  axy  m  *X  fxy(X>  6’  ^  etc*  (2*2) 

If  these  stress  components  are  required  to  leave  the  half-space  surface  and 
the  crack  surfaces  stress  free,  then  this  will  only  be  possible  for  an 
infinite  enumerable  set  of  values  of  X.  These  X  are  the  eigenvalues  and  the 
corresponding  functions  f^  the  eigenfunction^  the  numerical  determination 
of  which  is  the  aim  of  Benthem's  papers  (9-11).  Assuming  for  the  moment  that 
the  numerical  method  is  capable  of  determining  the  X's,  etc.,  with  sufficient 
accuracy,  there  still  remains  the  possibility  of  other  kinds  of  functions  of 
r  in  addition  to  rA  which  might  appear  in  a  complete  set.  Consider  the  illus¬ 
tration  the  simple  example  of  Appendix  1,  where  r  log  r  terms  appear.  In 
a  paper  on  conical  regions  in  elasticity  theory  (16),  Benthem  discusses  a 
theorem  which  says  that  if  X^  is  an  eigenvalue  so  is  -1  -X^  and  in  the  course 
of  that  work  he  is  aware  that  occasionally  terms  like  r^  In  r  appear  instead 
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of  r^, his  definition  of  in  terms  of  displacements,  i.e.t 


in 

such  a!ul 


u  f(X^,  9,  41).  If  such  alcl  *ccur  in  the  problem  under  discussion 

then  the  assumption  (2.2)  would  miss  them.  Also  it  is  clear  from  the 
work  of  Kondratev  (17)  that  rigorous  analysis  q^such  problems  can  be 
very  involved. 

Possible  values  of  X  are  constrained  by  the  condition  that  the  strain 
energy  be  finite,  this  means  that  in  a  neighborhood  of  the  vertex,  only  ** 
eigenfunctions  with 

Re  X  >  -3/2  (2.3) 

are  allowed.  Hence  within  the  framework  of  linear  elasticityf displacements 
u»«\rY,  Re  y  >  -1/2  which  become  infinite  at  the  vertex  are  not  ruled  out  by 
such  a  finite  energy  condition.  In  the  results  of  Folias  (14)  such  infinite 
displacements  were  reported,  i.e.,  he  reported  u  w*  r1^2  2v  which  becomes 
infinite  as  r  +  0  if  v  >  1/4.  Such  a  result  is  disturbing  to  physical 
intuition  and  was  criticized  for  this  reason  by  Benthem  and  Koiter  in  (15) 
(see  also  Appendix  2) . 

To  find  the  smallest  possible  eigenvalue  X  subject  to  the  constraint 
(2.3)  the  equations  of  elasticity  theory  are  applied  to  the  demi-spherical 
surface  r  =  1,  there  is  no  loss  of  generality  in  doing  this  since  along  a 
given  radius  all  stresses  are  assumed  proportional  to  r^  and  all  displace¬ 
ments  to  r^+^.  On  the  surface  r  *  1  the  range  of  the  two  remaining 
spherical  coordinates  is 

0  <  9  <  tt/2,  0  <  4  <  2v  (2.4) 


The  eigenvalue  problem  is  to  find  X,  fxx(X,  ®*  >  etc.,  as  solutions  of  the 

equations  of  elasticity  such  that  traction  free  boundary  conditions  hold  on 
the  crack  surfaces  at  r  -  1,  i.e.,  4*0,  and  2tt,  0  <  9  _<  x/2  and  on  the 
free  surface  9  *  tt/2,  0  <_  <j>  <_  2n.  In  (9)  and  (10)  Benthem  applies  a  method 
based  on  separation  of  variables  (9,  4).  In  (11)  he  uses  finite  differences 


and  in  (12)  and  (18)  Bazant  and  Estenssoro  use  a  finite  element  method 
In  (11)  Benthem  defines  a  new  variable 


C  -  01/2  so  0  <  ;  <  (ir/2)1/2  (2. 

and  divides  the  #  range  into  M  X  N  equal  meshes  AC, 

Tables  1,  2  and  Fig.  2  are  abstracted  from  his  paper  (with  our 
numbered  reference  notation) .  Symmetrical  (Table  1)  and  antisymmetrical 
(Table  2)  states  of  stress  refer  to  symmetry  (antisymmetry)  about  the 
plane  -  it. 


Fij.  I.  The  half-space  *  >  0  with  the  quaiter-inftnite  crack  }*0,j>0,;>0.  Cartesian  coordinates  z.  y,  i 

Spherical  coordinates  r,  0,  <t>. 


Table  I.  Eigenvalues  A  for  symmetrical  states  of  stress  <«r  —  r*! 


Poisson's 

ratio 

Obtained  with  Af  *  N 

10  x  30  14  x  42  16  x  40 

Benthem 

Refs.  (9,tO) 

Baiaaff 

Refs,  (u,!  »> 

Exact 

values 

-1 

-0.146 

-OIJJ 

-0.965 

-0.223 

-0241 

-0.75 

-0.4404 

-0.4516 

-0L5 

-0.4990 

-0i007  -030064 

-0.500*9 

0 

-049997 

-OS 

-OS 

MS 

-0.4035 

-04036 

-04*4 

OJ 

-04519 

-0.4519 

-0.4523 

-0452 

M 

-04141 

-04129 

-041J2 

-0413 

IS 

-03452 

-4-33*1  -03J67 

-0-33 11 

os 

0 

6 

03 

0.217 

0210 

0 

0.430 

0.410 

• 

0 

05 

05 

0.3 

0.664 

0.611 

OlS 

1 

• 

1 

tlaiant  and  Estenssoro  do  not  actually  give  values,  but  state  (bat  their  estrapoiadon  to  an  in  Orate 
number  of  meshes  gives  values  within  tH'*  from  those  of  previous  column  it)  s  »  s  0.4). 
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Tabk  2.  Eigenvalues  A  for  antimetrieal  states  of  stress 


PoissoA's 

Obtained  with  A l*N 
10x30  14*42 

Baiaat* 

Ref.  (u> 

Exact 

values 

• 

-l4wr 

-0.4999 

-15 

US 

-1S44S 

-05440 

-1545 

u 

>14074 

-05073 

-1591 

14 

-14291 

-052*4 

-4504 

15 

-16442 

-05442 

0 

-0.271 

US 

—0.157 

03 

-0.037 

14 

-0.030 

15 

• 

13 

0492 

0l3 

• 

0492 

os 

-1 

04*2 

OS 

tlaOOarooeously  is  stated  that  these  eigenvalues  are  double. 


The  results  shown  in  Tables  1  and  2  demonstrate  a  measure  of  agreement 
with  those  of  Bazant  and  Estenssoro.  Of  course,  both  sets  of  authors  use 
the  same  basic  assumption  (2.2),  the  subsequent  methods  used  are  however 
quite  different.  The  agreement  shown  in  the  tables  lends  support  to  the 
accuracy  of  the  results. 

-54 

X 

—  1 

-54 
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Fig.  I  Tkee-duneasional  stress  singularities  of  the  quart  er-m/lmte  crack  m  the  half-space  Solid  line. 
A -values  belonging  to  symmetrical  states  of  stress  (Table  I).  Broken  line.  A -values  belonging  to  antunetneai 
states  of  stress  iTable  21.  A.  B.  C.  D.  exact  poiou. 


i.  P.  Bomrai 


e,«f* •,  i  A  9,*t 

p.0.#  tphtficol  cooromotvs 


21 


Concluding  Remarks 

We  have  fixed  attention  on  the  Quarter-infinite  crack  shown  in  Fig.  1, 
since  this  is  perhaps  the  simplest  configuration  in  which  the  three 
dimensional  nature  of  the  singularity  when  a  crack  meets  a  free  surface  is 
present.  It  is  expected,  of  course,  that  the  nature  of  the  singularity  will 
be  precisely  the  same  for  a  finite  plate  provided  the  through  crack  is  at  a 
right  angle  to  the  plate  surface.  If  the  assumption  (2.2)  is  correct  then 
finding  the  lowest  eigenvalue  becomes  a  problem  of  doing  accurate  numerical 
analysis.  There  is  still  a  measure  of  disagreement  however,  c.f.,  the  paper 
by  Fujitani  (19)  as  well  as  the  previously  cited  papers  (13), (14).  We  have 
chosen  the  results  presented  in  section  2  for  the  reasons  given  there, 
however  the  reader  is  again  reminded  of  Appendix  2  for  the  other  side  of  the 
discussion.  As  well  as  the  eigenvalue  X  for  which  tables  of  results  are 
presented  here  there  are  graphical  presentations  of  the  eigenfunction  given 
in  (11)  and  much  more  detail  than  we  have  presented  here.  Furthermore, 
Bazant  and  Estenssoro  (12)  have  applied  their  numerical  method  to  the 
problem  of  a  through  crack  which  meets  the  surface  obliquely  to  find  the 
eigenvalues  there.  In  this  case  there  can  be  complex  eigenvalues. 

In  conclusion,  there  now  exist  numerical  results  for  the  eigenvalue 
X  in  assumed  stress  distributions  of  the  form  a  r^  (r  spherical  radius) 
for  the  three  dimensional  singularity  of  a  through  crack  at  a  free  surface. 
Furthermore,  there  is  some  measure  of  numerical  agreement  between  two  sets 
of  authors,  i.e.,  (11)  and  (12).  The  associated  eigenfunctions  with  this 
dominant  singularity  Re  X  >  -3/2  have  been  obtained  numerically  and  pre¬ 

sented  graphically  in  (11)  for  example.  Note  that  other  related  three 
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dimensional  singularities  have  also  been  considered, 
inclusions  of  a  wedge-plan  form  see  (1)  and  (12)  for 
of  these  papers.  Finally  in  references  (20)-(27)  we 
attempts  at  modelling  three  dimensional  cracks. 


e.g.,  cracks  and 
references  to  some 
cite  less  sophisticated 
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Appendix  1 

Consider  Laplace’s  equation  in  plane  polar  coordinates, 

r  £  <r  +  A  - 0  (A.D 

3r  ar  30z 

in  the  region  0  <_  0  <_  3it/2,  0  <  r  <  =  of  the  (r,9)  plane  subject  to  the 
boundary  conditions 


*  0  on  8  =  0  for  all  r 


34>  _ 

ae 


r 


for  0  <  r  <  1 


=  0  for  r  >  1 
and  i(  +  0  as  r  + 


on  8  =  3tt/2 


(A.  2) 


g-1 

Multiplying  (A.l)  by  r  and  integrating  over  r  from  o  to  »  gives  the 
equation 


+  s 


(A. 3) 


for  the  Mellin  transform  of  ^  $  rS  1  dr.  Solving  (A. 3)  subject  to 

the  transform  of  the  boundary  conditions  (A. 2)  and  using  the  inversion 
theorem  of  the  Mellin  transform  leads  to  the  result 


=  _i_  7*°°  Sin  s  6  .  r~S  dg 

2lT*  y-i°°  s(s  +  1)  Cos(-^p-) 


(A. 4) 


The  condition  that  $  be  finite  at  the  origin  r  ■  0  together  with  ♦  +  0  as 
r  -*•  oo  designates  the  strip  of  regularity  of  <f>(s)  so  that 


-1/3  <  y  <  1/3 


(A. 5) 


The  solution  near  r  «  0  can  be  deduced  from  (A. 4)  by  closing  the  contour 
in  the  left  half-plane  and  evaluating  the  residues.  These  poles  are  simple 
with  the  exception  of  a  double  pole  at  s  *  -1.  The  final  result  for  r  <  1,  is 
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(A.6) 


-  3(— 1)  (2N+l)/3 

N^0  ir(2N+l)(N-l) 


(2N+l)/3  Sin{(~±i)e} 


+  (0  Cos  0  +  log  r  .  Sin  9  -  Sin  6) 

2 

That  the  last  expression  is  a  solution  of  7  <}>  =  0  can  be  easily  checked. 
Tne  interesting  feature  from  the  present  point  of  view  is  the  presence  of 
the  log  r  term  which  may  not  have  been  suspected  a  priori. 


Appendix  2 


Summary  of  Discussion  of  T.  Kawai 
and  E.S.  Folias,  Austin,  Texas 
March  29,  1979 

The  discussion  centered  on  two  points.  The  first  was  the  discrepancy 
between  J.P.  Benthem's  and  Kawai' s  solution  and  the  second  was  an  inquiry 
into  whether  the  3-D  problem  possessed  a  separable  solution  in  a  form  similar 
to  Williams  2-D  eigenfunction  expansion.  Benthem  and  Kawai  both  chose  similar 
starting  points  for  their  attempts  to  solve  the  3-D  crack  problem.  They  both 
assumed  the  solution  to  be  separable  in  the  form  r*F(x,  0,  <j>)  where  r,  6, 4>  are 
spherical  coordinates  and  X  is  an  eigenvalue.  While  Benthem  constructed  a 
solution  which  satisfied  the  stress  free  boundary  conditions  on  the  crack 
flanks  and  then  satisfied  the  stress  free  surface  conditions,  Kawai  began  by 
satisfying  the  stress  free  surface  conditions  and  then  satisfying  the  condi¬ 
tions  that  the  crack  flanks  be  stress  free.  Kawai  felt  it  was  easiest  to  find 
functions  which  were  zero  on  the  plate  surface.  The  approach  of  Kawai  pro¬ 
duced  eigenfunctions  which  Benthem  did  not  find.  Kawai  noted  that  he  had  a 
student  do  additional  work  on  the  problem  but  was  unable  to  resolve  this 
discrepancy  and  had  not  published  this  work. 

Folias  pointed  out  that  Benthem's  1972  paper  contained  an  incomplete  set 
of  functions  which  Benthem  had  added  to  in  1976.  However,  the  question  of 
completeness  of  Benthem's  solution  has  not  been  resolved.  Folias  had  at  an 
earlier  meeting  in  San  Francisco  suggested  that  Kawai  show  Benthem's  solution 
was  incomplete.  By  showing  one  of  the  eigenfunctions  which  Benthem's  method 
had  not  produced  satisfied  both  the  differential  equation  and  boundary  condi¬ 
tions,  Kawai  would  be  able  to  demonstrate  that  the  solution  of  Benthem  was 
still  incomplete.  Kawai  has  not  attempted  this  and  seemed  reluctant  to  pursue 

29 


'fc  \ 


r 


the  point.  The  reason  for  this  reluctance  was  two-fold;  (i)  lack  of  a 
capable  graduate  student  to  follow  up  the  initial  work,  and  (ii)  the  adverse 
reaction  of  the  fracture  community  to  ideas  challenging  Benthem's  hypothesis. 
It  was  also  noted  that  the  singularity  which  Folias  had  found  was  reproduced 
by  Kawai's  technique.  The  result  was  found  to  be  caused  by  a  singular  condi¬ 
tion  of  the  matrix  in  the  eigenvalue  determination.  This  indicates  that  the 
functions  of  Folias  are  most  probably  zero  on  the  surface.  Folias 

and  Kawai  agreed  that  Benthem's  method  was  not  systematic  and  Folias  con¬ 
cluded  that  on  the  basis  of  Kawai's  adopting  a  more  systematic  approach  to  the 
solution  his  results  were  more  likely  to  be  complete. 

Further,  work  of  J.  Newman  also  furnishes  indirect  evidence  that  the  work 
of  Folias  and  Kawai  is  correct.  Newman  found  that  near  the  surface  of  a  plate 
the  stress  intensity  factor  rises  slightly  as  a  function  of  the  plate  thick¬ 
ness  for  a  crack  length  to  thickness  ratio  much  less  than  one.  This  is  in 
agreement  with  the  work  of  Folias  and  Kawai.  Kawai  also  had  found  that  the 
first  eigenvalue  changes  as  the  plate  thickness  increases.  Elasto-plastic 
finite  element  work  of  Kawai  (presented  at  Swansea  in  January)  on  the  cracked 
plate  also  indicated  that  the  yielding  begins  on  the  plate  surface  and  spreads 
into  the  interior.  From  this  he  inferred  that  the  strength  of  the  singularity 
changes  near  the  surface.  A  final  problem  was  proposed  by  Kawai  as  an  alter¬ 
native  method  for  determining  the  strength  of  the  singularity.  The  problem 
is  equivalent  to  the  Inglis  solution  of  an  elliptic  hole.  The  solution  would 
be  sought  for  the  elliptic  hole  in  a  plate  and  the  minor  axis  would  be  allowed 
to  approach  zero.  In  this  way  some  insight  may  be  gained  concerning  the  3-D 
crack  singularity. 


The  problem  of  separability  of  the  solution  was  discussed.  Kawai  pointed 


to  the  fact  that  the  problem  of  the  conical  cavity  in  a  half  space  was 
separable  and  by  inference  the  crack  solution  should  be  also.  Folias  had 
some  reservations  and  pointed  out  that  the  eigenvalue's  slow  rate  of  con¬ 
vergence  in  the  numerical  solution  technique  is  possible  evidence  that  the 
solution  may  not  be  separable.  This  Kawai  admitted  may  indeed  be  true  and 
thought  that  while  the  solution  for  the  stress  analysis  may  be  separable, 
separability  may  not  apply  to  the  solution  of  the  singularity.  Folias 
pointed  out  that  spherical  coordinates  were  not  the  natural  coordinates  for 
the  problem,  and  while  the  solution  may  not  be  separable  in  spherical  coor¬ 
dinates,  it  may  be  in  some  other  system.  Folias  noted  that  it  may  be  possible 
to  obtain  the  dominant  singularity  by  a  separation  technique.  Kawai  agreed 
that  expansions  in  Legendre  polynomials  may  not  be  appropriate  but  another 
special  function  such  as  Bessel  functions  or  a  generalized  function  may  work. 
Agreement  was  reached  that  the  final  solution  may  not  be  obtainable  through 
a  generalization  of  the  Williams  solution  techniques.  Thus  to  date  the 
dominate  singularity  is  not  known  for  the  3-D  problem  and  may  not  be 
obtainable  through  standard  techniques. 


31 


SETEC  # 
October 


Appendix  4.2 


A  Note  on  Stress  Concentrations 
and  Energy  Release  Rates 


! 

i 


Dr.  Colin  Atkinson 
Imperial  College,  London 
Visiting  Professor,  University  of  Pittsburgh 


i 


j 


i 


i 


Science  of  Fracture  Research  Project 
AFOSR  Contract  //F49620-78-C-0101 


i 


\  ' 

I  ' 


M 


f 

I 


School  of  Engineering 
University  of  Pittsburgh 
Pittsburgh,  PA  15261 


i 


■f 


1 


80-049 

1980 


mCMDlM  PAfl I  HiiB 


1.  Introduction 


There  are  many  situations  where  pre-existing. stress  concentrators  exist 
in  a  solid,  either  because  of  particular  design  requirements  or  because  of 
defects  which  occur  during  fabrication  or  subsequent  handling.  One  particular 
kind  of  defect  is,  of  course,  a  crack  and  much  attention  has  been  given  to 
this,  leading  to  the  well  developed  subject  of  fracture  mechanics.  For  a 
recent  review  of  stress-concentrations  and  fracture  mechanics  see  (1) .  Our 
concern  here  is  not  with  the  crack  defect  per  se  but  rather  with  its  relation 
and/or  comparison  with  other  stress  concentrations  such  as  corners  or  holes. 
However,  since  the  concept  of  'energy  release  rate'  is  useful  in  discussing 
crack  extension  it  is  perhaps  natural  to  ask  whether  an  energy  balance 
approach  can  in  some  analogous  way  quantify  the  effect  of  other  stress 
singularities  and  'perhaps'  explain  the  transition  from  a  corner  (without  a 
crack  starter)  to  a  crack-like  singularity.  We  use  the  adverb  perhaps  here, 
because  it  is  natural  when  using  energy  arguments  to  discuss  the  change  in 
energy  with  a  change  in  configuration  of  the  defect,  and  it  may  be  that 
configuration  changes  other  than  crack  initiation  are  energetically  favorable. 
A  useful  tool  for  such  energy  calculations  is  the  elastic  energy  momentum 
tensor  which  has  been  used  for  expressing  the  force  on  a  dislocation  or 
point  defect  (Eshelby  (2)  and  (3))  and  for  calculating  the  energy  flow  into 
a  moving  crack  tip  (Atkinson  and  Eshelby  (1)).  We  shall  give  a  brief  account 
of  this  in  Section  2  which  is  based  to  some  extent  on  a  paper  by  Eshelby  (5). 
In  Section  3  we  discuss  the  implications  of  the  formulae  derived  in  Section  2 
when  applied  to  cracks  or  blisters  forming  on  comers  or  holes.  Our  conclu¬ 
sion  will  be  that  the  initiation  of  a  small  crack  from  a  hole  will  lead  to  an 
energy  release  rate  or  crack  extension  force  which  usually  tends  to  zero  as 
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the  crack  length  tends  to  zero.  On  the  other  hand,  there  will  usually  be  a 
finite  effective  normal  force  (per  unit  surface  area)  on  a  hole  boundary 
(c.f.  equation  (2.12)).  However,  although  this  finite  force  can  be  used 
to  give  a  finite  energy  release  rate  in  erecting  a  blister  (c.f.  for  example 
(2.10))  this  must  be  balanced  by  the  energy  required  to  produce  a  blister. 
This  required  energy  may  be  large  since  removal  of  material  is  required 


for  blister  formation 


2.  Energy  Change  Calculations  and  The  Energy  Momentum  Tensor 

For  completeness,  we  give  here  results  (abstracted  from  (5))  which  apply 
to  finite  deformation  with  a  general  stress-strain  relation,  though  our 
subsequent  applications  will  be  to  infinitesimal  elasticity.  We  use  rectan¬ 
gular  coordinates  X^  to  label  particles  in  the  reference  configuration 
(initial  unstrained  state).  The  cartesian  coordinates  (in  the  same 
coordinate  system)  refer  to  the  current  configuration  to  which  the  particle 
X^  is  moved  on  deformation  by  a  vector  displacement  such  that 


W  =  xi<V  -  Xi  (2-1) 

The  nominal  (Piola-Kirchof f  or  Boussinesq)  stress  p  is  defined  so  that 
p^j  dS  is  the  component  parallel  to  the  X^  axis  of  the  force  on  an  element  of 
area  which,  before  deformation,  was  an  element  of  area  dS  perpendicular  to 
the  Xj  axis.  Thus  if  W  is  the  strain  energy  in  unit  initial  volume  then 


3W(u  ) 


ij 


3u 


Ui 


(2.2) 


where  u.  ,  =  3u. /3X; 
1,  J  1 


In  the  absence  of  body  forces  the  equations  of  equilibrium  are 


(2.3) 


Note  that  for  finite  deformation  the  p^  are  not  symmetric. 

In  terms  of  the  energy  density  W,  the  energy  momentum  tensor  is  defined  as 


?li  ’  W5£j  "  Pij  ui ,1  (2-4) 

being  the  dirac  delta  function.  It  can  be  shown  by  direct  calculation  that 
wherever  W  does  not  depend  on  X^  its  divergence  vanishes,  i.e.. 
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(c.f.  (5))  This  latter  property  still  holds  if  the  energy  momentum  tensor 
is  defined  as 


3L 


Ij  3u 


i  ,t 


u.  *  “  1*6  a  . 


(2.6) 


where  L  is  a  Lagrangian.  The  Euler  equations,  such  that  the  integral  of  L 

over  a  region  of  X  space  is  an  extremum,  then  correspond  to  the  equilibrium 
m 

equations  (2.3)  when  L  =  -W.  With  this  formalism  more  general  situations 
can  be  considered  see  (5 — 10)  for  some  application. 

We  are  primarily  interested  here  with  the  elastic  situation  and  hence 
equation  (2.4)  and  its  connection  with  the  force  on  a  defect.  In  (5)  the 
situation  shown  in  Fig.  1  is  considered;  Fig.  1  represents  a  body  subject  to 
surface  loading  and  containing  a  defect  together  with  other  defects  D^. 

If  suffers  a  small  vector  displacement  6?  it  is  shown  that  the  change 
in  total  energy  of  the  system  (elastic  energy  plus  potential  energy  of  the 
loading  mechanism)  5Etot  is  given  as: 


6E 


tot 


-*h  h 


(2.7) 


where 


(2.8) 


S  is  any  surface  surrounding  and  isolating  it  from  D2.  F^  is  thus  the 
force  on  the  defect,  where,  following  (5),  the  rate  of  decrease  of  the  total 
energy  of  the  system  with  respect  to  a  parameter  is  defined  as  the  generalized 
force  acting  on  that  parameter.  Furthermore  from  (2.5)  S  may  be  deformed  with¬ 
out  altering  F^  provided  it  remains  in  material  where  W  does  not  depend 
explicitly  on  (S  still  isolates  from  D^,  of  course). 
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It  is  also  possible  to  use  the  energy-momentum  tensor  to  calculate  the 
force  on  an  interface.  Consider  the  situation  shown  in  Fig.  2  where  we 
assume  the  material  inside  B  is  a  cavity  with  surface  S.  S  may  be  closed 
05  as  indicated  by  the  dotted  lines,  may  extend  to  the  specimen  surface.  It 
is  now  supposed  that  S  has  developed  a  blister  changing  it  to  S',  geometrically 
this  can  be  specified  by  erecting  a  small  vector  <55^  at  each  point  of  S.  It 
can  then  be  shown  see  (5)  that  the  change  in  the  total  energy  due  to  the 


formation  of  the  blister  is 


--  f 


dS. 


(2.9) 


where  P  ^  is  the  energy  momentum  tensor  on  the  A  side  of  S  and  the  normal  in 
dSj  =  n^  dS  is  directed  outwards  from  B  to  A.  The  results  derived  in  (5)  also 
apply  to  the  situation  when  B  may  be  'transformed'  elastic  region  which  may  have 
different  elastic  properties  but  such  that  the  displacement  is  continuous  across 
S  and  then  (2.9)  is  written  as 


65  £  (00  6^.  -  PytUi^)  dSi 


(2.10) 


where  (w]  means  the  jump  in  W  across  S.  Using  the  result 


J  U  . 

Lui,J  -  flTJ  nl 


(2.11) 


where  3/3n  denotes  differentiation  along  the  normal  it  is  possible  to  interpret 
(2.10)  as  saying  that  there  is  an  effective  normal  force 


F  =  [W]  -  T  . 


(2.12) 


per  unit  area  of  intarface  where  T  is  the  surface  traction  at  the  interface. 
Note  that  for  a  stress  free  cavity  this  result  says  that  F  becomes  the  energy 
density  at  the  surfact  (a  necessarily  positive  quantity),  hence  any  increase 


in  volume  leads  to  a  decrease  of  E  (this  result  has  been  derived  in  the 

tot 

small-displacement  case  by  Rice  and  Drucker  (11)). 
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Thus  we  have  two  Important  expressions  for  the  change  in  total  energy 
with  respect  to  different  kind  of  configurational  change.  The  first  given 
by  (2.7)  and  (2.8)  applies  directly  to  a  crack  like  defect.  When  the  crack 
extends  in  its  plane  in  the  1  direction  say  then  given  by  (2.8)  is 
equivalent  to  G^the  crack  extension  force.  Moreover,  by  its  nature  this 
crack  extension  involves  no  removal  of  material.  Expression  (2.9)  on  the 
other  hand  extends  a  cavity  by  forming  a  blister  which  does  remove  material. 
Even  though  a  physical  mechanism  can  be  involved  for  this,  e.g.,  diffusion  or 
perhaps  multiple  cracking  it  does  seem  to  have  a  different  character.  Never¬ 
theless,  applying  the  blister  argument  to  a  thin  cavity  can  be  shown  to  give 
the  same  result  as  that  produced  by  applying  (2.8)  in  the  limit  that  the 
cavity  is  very  thin  (c.f.  (5)). 

In  addition  to  the  above  uses  of  the  energy-momentum  tensor  in  calculating 
energy  release  due  to  defect  motion  or  blister  formation,  it  is  also  possible 
to  calculate  energy  release  due  to  cavity  expansion  or  rotation.  To  do  this 
some  useful  tools  are  path-independent  integrals  first  discovered  by  Gunther 
(12)  and  rediscovered  more  recently  by  Knowles  and  Sternberg  (13).  These 
can  be  written  as 


hi  *  J(XkP£j  '  X£Pkj  +  \  p£j  ■  ul  pkj}  dSj 

and 

M  -  -  i  u£  p£j>  dSj 


(2.13) 


(2.14) 


The  are  path-independent  for  a  non-linear  homogeneous  isotropic  medium, 
whereas  the  integral  M  is  path  independent  for  a  homogeneous  linear  anisotropic 
medium  (c.f.  (9)).  In  particular  for  plane  strain  it  can  be  shown  that  the 
expression  (2.14)  reduces  to 


M 


dS 


(2.15) 
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See  Eshelby  (14)  for  some  applications  of  this  result.  In  (15),  Atkinson 
and  Leppington  show  how  an  amended  M  integral  (which  is  not  path-independent) 
can  be  deduced  for  materials  with  couple-stresses  and  how  it  can  be  used  to 
prove  a  useful  theorem  about  energy  release  rates  in  such  materials. 

Furthermore,  the  L  and  M  integrals  defined  in  (2.13)  and  (2.14)  have 
been  shown  to  be  related  to  energy-release  rates  associated  with  cavity  or 
crack  rotation  or  expansion  by  Budiansky  and  Rice  (16) . 
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3.  Simple  Applications 

In  the  last  section  we  have  outlined  how  the  energy  momentum  tensor 
can  be  used  to  give  a  fairly  unified  approach  to  calculations  of  the  force 
on  a  defect,  or  change  in  total  elastic  energy,  under  a  variety  of  con¬ 
figuration  changes.  In  this  section  we  shall  explore  by  means  of  simple 
examples  the  energy  released  when  (i)  a  cavity  changes  position  slightly, 
due  to  a  slight  displacement  of  a  position  coordinate  or  (ii)  a  small 
blister,  such  as  discussed  in  Section  2,  forms  on  a  cavity  surface.  For 
comparison  we  shall  also  consider  the  situation  when  a  crack  initiates  from 
a  cavity.  The  important  point  to  notice  here  is  that  these  calculations 
tell  us  only  what  the  release  of  elastic  energy  is,  whether  or  not  such 
changes  can  occur  in  practice  depends  on  the  amount  of  energy  required  to 
effect  them. 

Ex.  1  The  circular  hole  in  an  infinite  plate 

Consider  a  small  circular  hole  of  radius  a  In  an  infinite  plate  (under 
plane  stress  or  plane  strain  conditions)  and  a  tensile  loading  at  infinity 
in  the  1  direction  such  that  at  r  » 


and 


Pll  a  S*  p12  *  0  "  P22 


(3.1) 


p  »  p  ,  *  0  on  the  hole  boundarv 
*rr  rr 9 


A  well  known  result  in  this  case  is  that  the  stress-field  given  in  polar 
coordinates  (r,9),  with  origin  at  the  hole  center,  is 


i 
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prr  -  |  (1  -  a2/r2)  +  |  (1  +  3a4/r4  -  4a2/r2)  Cos  28 

Pg 0  -  |  (1  +  a2/r2)  -  -|  (1  +  3a4/r4)  Cos  26  (3.2) 

Pre  -  -  |  (1  -  3a4/r4  +  2a2/r2)  Sin  28 


Thus  on  the  boundary  of  the  hole  {>  -  p  -  0  and 

rr  r  u 

Pgg  *  S(1  ~  2  COS  2  0) 


(3.3) 


Calculating  as  defined  in  Equation  (2.8),  by  evaluating  the  integral 
around  a  closed  contour  enclosing  the  hole,  gives  the  force  on  the  hole 
(recall  that  5Etot  =  F^,  Equation  (2.7))  when  the  origin  of  the  hole 

is  moved  by  6?^.  The  value  turns  out  to  be  zero  which  is  not  so  surprising 
since  any  position  of  the  hole  in  the  infinite  plate  is  as  good  as  any  other. 
However,  if  we  imagine  that  a  blister  forms  on  the  hole  we  can  calculate  the 
energy  change  using  Equations  (2.9)  or  (2.10)  of  the  last  section.  An 
elementary  calculation  gives 


f2  aS2  1 

\  a  W  d 8  =  -  (30  -  2  Sin  26  +  -4  Sin  4  9) 


(3.4) 


and  hence  if  a  blister  forms  on  -ir/4  <8  <  ir/4  say  we  have  a  release  of  energy 

(aS^/2E)  (3ir/2  -  4)  |  ^  |  whereas  if  a  blister  forms  on  tt/4  <  8  <  3w/4  the  energy 
released  is  (aS^2E)(3ir/2  +  4)  |  |  making  it  energetically  favorable  for  the 

blister  to  form  on  the  side  at  right  angles  to  the  applied  stress,  an 
intuitively  obvious  result.  The  important  point  here  is  that  these  results 
give  rise  to  a  finite  force  per  unit  area  of  interface.  Opposing  this,  of 
course,  is  the  amount  of  energy  required  to  form  a  blister  which  may  be  very 
large. 
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Suppose  we  consider  an  infinitesmal  'rectangle*  of  blister  to  be  formed 


near  0  ■  it/ 2,  then  integrating  (3.4)  with  9^  *  n/2  ±  59  and  applying  (2.9)  gives 
the  change  in  total  energy  as  ^ 


5E 


tot 


(3.5) 


5A 


where  6A  *  2a  60  .  5?  is  the  area  of  the  blister  which  has  been  formed.  It  is 
of  interest  to  compare  this  result  with  that  of  a  crack  of  length  e  initiated 
on  the  hole  boundary  at  9  *  ir/2  at  right  angles  to  the  applied  stress.  For 
this  case  it  is  not  difficult  to  show  that  if  e  <<  a  then  the  energy  release 
rate  or  crack  extension  force,  G,  is 


G  -  e(1.1215)2  it  .  ^ 


(3.6) 


Thus  as  the  crack  length  e  tends  to  zero  so  does  the  energy  release  rate. 
This  result  is  to  be  contrasted  with  the  result  (3.5)  recalling  that  5Etot 
*  -F6£  where  F  is  the  force  on  the  defect  (the  blister  in  this  case).  Thus 


F 


9S2  a69 
2E 


(3.7) 


and  as  the  blister  gets  very  thin  (almost  crack  like)  F  tends  to  zero  with 
a6 0.  Note  that  for  plane  strain  E  (Young's  modulus)  is  replaced  by  E/(l  -  v^) 


Ex.  2  A  step  on  the  surface  of  a  solid 

For  simplicity,  we  consider  anti-plane  deformation  in  a  half-space  x,,  >  0 
with  a  step  of  height  h  situated  on  X2  ■  h,  x^  >  0  and  an  external  shear  stress 
p^3  ■  o  applied  at  infinity. 
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This  problem  can  be  solved  by  conformal  mapping  and  has  been  treated  by 


Smith  (17).  The  following  results  are  obtained  near  the  corner  (0,h). 


P13  “  Cos  e^3 

(3.8) 

p23  *  o(3^)  Sin  9/3 


(p^2  and  p23  are,  of  course,  the  only  non-2ero  stresses  in  anti-plane  strain), 

r  and  0  are  polar  coordinates  centered  on  the  comer  (0,h),  0  ■  0  corresponds 

to  the  surface  x2  *  h,  x^  >  0  and  0  *  3ir/2  to  the  edge  x^  =  0,  0  <  x2  <  h. 

2  2 

The  energy  density  is  now  &/2y)(p*  +  P23)  and  is  infinite  at  the  comer  r  =  0. 

However,  integrating  along  0=0,  with  respect  to  r  will  give  a  finite  energy 
release  rate  (per  unit  surface  area)  for  formation  of  a  blister  in  the  x2 
direction.  Similar  conclusions  apply  to  the  edge  x^  =  0,  0  <  x2  <  h.  Note 
that  in  either  case  the  effective  force  per  unit  surface  area  of  the  step  will 
be  proportional  to  oV/3.  It  is  interesting  to  compare  this  with  the  energy 
release  rate,  G,  associated  with  a  crack  of  length  e  emitting  from  (0,h)  in 


the  direction  x„  to  the  point  (0,'n  +  e).  For  e  <<  h  it  is  shown  in  (17)  that 
£ 

this  is 


a  W/3  h2/3  e1/3/2„ 


(3.9) 


1/3 

This  tends  to  zero  as  e  tends  to  zero  like  e  so  the  presence  of  the  comer 
although  increasing  dG/dt  compared  with  the  corresponding  result  for  a  hole 
(Equation  3.6)  still  does  not  provide  a  finite  crack  extension  force  for  crack 
initiation.  The  best  we  can  do  in  this  regard  is  to  employ  the  force  on  an 
interface  idea  discussed  in  Section  2  (e.g..  Equation  2.10).  The  result  is 
often  a  finite  effective  force  per  unit  area,  but  as  indicated  earlier  this  is 


offset  by  the  energy  required  to  remove  material  for  such  blister  formation  to 
be  possible. 
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Perhaps,  we  must  be  content  when  discussing  fracture  from  holes  or 
corners  to  discuss  the  increased  likelihood  of  fracture  due  to  results  like 
(3.9).  Results  of  this  general  type  have  been  derived  for  the  more  realistic 
plane  strain  situation  by  Atkinson  (18)  and  in  work  in  preparation  by 
Atkinson  and  Avila  (19) . 
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ABSTRACT 

Further  evidence  is  given  to  show  that  the  classical  line  crack  problem 
formulated  for  non-local  elastic  media  with  ’smooth'  moduli  may  have  no 
solution  with  non-infinite  crack  opening  displacements.  In  particular^  claims 
of  finite  crack  tip  stresses,  etc.,  with  this  theory  are  likely  to  be  incorrect. 
We  give  numerical  results  which  support  our  earlier  theoretical  conclusions 
(1)  and  (2).  Computations  are  made  for  a  modulus  which  has  an  elastic 
component  as  well  as  the  ’smooth’  non-local  moduli  referred  to  above. 

INTRODUCTION 

In  two  recent  papers  (1)  and  (2)  we  have  discussed  calculations  made  by 
Eringen  and  co-workers  of  stresses  around  line  cracks  in  a  non-local  elastic 
continuum  (3,4,5).  Our  main  conclusions  in  (1)  were  (i)  that  the  approxi¬ 
mation  used  (3,4,5)  had  a  non-uniform  character  and  hence  would  probably  give 
unreliable  results  and  (ii)  that  problems  of  the  type  solved  numerically  by 
Eringen  et  al.  in  (3)  could  be  shown  to  have  no  solution  with  finite  (i.e., 
non-infinite)  crack  opening  displacements.  In  (2)  we  extended  these  results 
to  plane  strain  situations  and  constructed  solutions  for  a  non-local  theory 
with  moduli  which  had  a  delta  function  component  in  addition  to  the  'smooth' 
moduli  used  by  Eringen  et  al.  With  these  amended  non-local  moduli  it  is  not 
difficult  to  see  that  solutions  should  exist.  However,  the  line  crack  problem 
will  have  the  usual  square  root  stress  singularity  at  its  tip.  In  principle 
it  is  possible  to  do  calculations  with  such  moduli  and  hence  deduce  stress 
intensity  factors  for  media  with  such  non-local  constitutive  laws.  Thus 
stress  intensity  factors  could  then  be  used  just  as  in  classical  fracture 
mechanics.  The  question,  of  course,  is  whether  sensible  (specific)  constitu¬ 
tive  relations  can  be  obtained  which  are  relevant  to  real  materials.  The 
constitutive  laws  chosen  by  Eringen  et  al.  did  attempt  to  connect  with  lattice 
theory  through  dispersion  relations,  but  as  we  have  suggested  above  and  in  (1) 
and  (2)  there  are  difficulties  with  the  application  of  such  a  theory  to  the 
classical  line  crack  problem. 
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In  a  private  communication  concerning  (1)  Dr.  Eringen  has  claimed  that 
for  the  non-local  theory,  boundary  conditions  cannot  be  applied  in  the 
classical  sense.  We  do  not  know  precisely  what  this  means  but  he  appears 
(in  retrospect)  to  be  trying  to  avoid  just  those  regions  of  non-uniformity 
that  we  show  exists  in  his  approximation  scheme.  He  also  makes  the  comment 
that  our  proofs  in  (1)  are  not  rigorous  mathematical  ones.  To  counter  this 
criticism  we  present  here  in  more  rigorous  form  one  of  these  proofs;  filling 
in  the  gaps  in  the  other  proofs  of  (1)  and  (2)  is  left  to  the  reader.  It  may 
be  thought  that  the  application  of  non-local  theories  to  the  fracture  problem 
is  only  of  academic  interest,  however  some  publicity  has  been  given  to 
non-local  calculations  which  predict  finite  stresses  at  crack  tips  (6)  and  it 
would  seem  useful  to  attempt  to  clarify  the  matter,  perhaps  this  conference 
can  contribute  something  in  this  direction? 

The  plan  of  this  paper  is  as  follows:  in  Section  2  we  state  equations  of 
equilibrium  and  constitutive  equations  for  the  linear  non-local  theory,  and 
show  how  mathematical  rigor  can  be  added  to  the  non-existence  proofs  of  (1) . 
Finally  in  Section  3  we  describe  some  numerical  results  for  a  crack  in  a 
non-local  elastic  solid  with  modulus  of  the  form 

a(|x|)  =  b  6 ( |x| )  +  (1  -  b)  a1(|x|)  (1.1) 


(compare  Section  2  for  a  fuller  description  of  such  constitutive  laws)  b  is  a 
constant  0  <  b  <  1  and  a-j(|xJ)  is  a  smooth  function  of  x  depending  on  a  small 
parameter  a  (the  lattice  parameter).  We  have  shown  in  (2)  that  for  moduli  of 
this  general  type  when  a  (crack  length)  tends  to  zero,  the  crack  tip  displace¬ 
ment  and  stress  field  become  respectively,  the  classical  crack  tip  elastic 
displacement  field  multiplied  by  b~%  and  the  classical  crack  tip  stress  field 
multiplied  by  b'S.  When  b  -*■  0  equation  (1.1)  is  a  modulus  of  the  type  con¬ 
sidered  by  Eringen.  In  this  case  our  asymptotic  results  say  that  the  crack 
opening  displacement  tends  to  infinity  like  b~ ^  which  supports  the  non¬ 
existence  results  of  (1)  and  (2).  As  can  be  seen  in  Section  3,  the  numerical 
results  confirm  this  asymptotic  result. 

2  ANALYSIS 


The  basic  equations  of  linear,  homogeneous,  isotropic,  non-local 
elasticity  with  vanishing  body  and  inertia  forces  can  be  written  as: 


tk€,  k 
tkL 


f 


UT(|x-x’  j)  err(x*)Sj+.  2p’(|x  -  x’|) 


(2.1) 

(2.2) 


with 

V.  =  1  (uk,£  +  u(,k5 


(2.3) 
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In  the  above  equations  the  summation  convention  with  respect  to  repeated 
indices  is  used,  and  the  indices  run  from  1  to  3  in  general.  Clearly  from 
(2.2)  the  stress  tk£(x)  at  a  point  x  depends  on  the  strains  eit£(x')  at  all 
points  of  the  body.  The  integral  in  (2.2)  is  over  the  volume  V  of  the 
body  enclosed  within  a  surface  3V. 

Various  expressions  for  the  non-local  moduli  X'(|^|)  and  y'(|x|)  have 
been  used  in  the  cited  papers.  Each  of  these  is  of  the  form 

(X',y')  =  (X,u)  a(|x'  -  x|)  (2.4 

where  (A,p)  are  constants  and  the  function  a  is  normalized  so  that 


J"a(|x’  -  x|)  dV(x')  =  1 


the  integration  being  over  the  whole  body.  Expressions  for  a  which  have  been 
obtained  by  Eringen  et  al*by  comparison  with  (one-dimensional)  dispersion 
relations  of  lattice  dvnamics  are: 


(i)  «(|x'  -  xj)  =  aQ(a  -  |x’  -  x|)  for 


x'  -  x  <  a 


(2.6) 


if  lx*  -  x|  >  a 


a  is  the  lattice  parameter  and  a  is  a  normalization  constant  deduced  from  (2.5) 


(ii)  a(  |x'  -  x|)  =  aQ  exp{-  (|)  (x£  -  xk)  (x^.  -  xk)  } 


where  the  summation  convention  applies  to  k,  6  is  a  constant.  See  (3)  and  (4) 
for  a  justification  of  the  use  of  these  expressions. 

Other  examples  for  a  were  considered  in  (1)  for  illustrative  purposes. 

One  of  these,  for  example,  is 


(iii)  a(|*|)  a~  ~ 

(  l&i  +  a  ) 


(2.8) 


Note  that  each  of  (2.6),  (2.7)  or  (2.8),  with  the  correct  normalizing  factor 
a0,  have  the  property  that  formally  letting  a  -*•  0  implies  oc C  | -  x  | )  -*■ 

5(|x’  -  x|).  This  limiting  case  is,  of  course,  classical  elasticity,  however 
with  a  ^  0  moduli  such  as  (2.6),  (2.7)  and  (2.8)  behave  quite  differently. 
Later  as  a  contrasting  example  to  the  moduli  of  (2.6),  (2.7)  or  (2.8)  we  shall 
consider  moduli  of  the  form 


(iv)  a(|x'  -  x|)  =  b  <5 ( | x '  -  x|)  +  (1  -  b)  cx^lx’  -  x|) 


(2.9) 


where  b  is  a  constant  0  <  b  <  1  and  is  a  function  of  the  ’smooth'  type 
discussed  above  {(2.6),  (2.7)  or  (2.8)}.  Expressions  of  type  (2.9)  have  been 
suggested  by  Kroner  (7). 
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Using  equations  (2.1),  (2.2)  and  moduli  (2.6)  or  (2.7)  in  the  two 
dimensional  form  applicable  to  plane-strain,  Eringen  and  co-workers  in  a 
series  of  papers  (3,4,5)  have  treated  boundary  value  problems  of  line  cracks 
under  various  loading  conditions,  plane  strain  shear,  plane  strain  tension  and 
anti-plane  shear.  In  (3)  they  treat  a  model  problem  via  numerical  solution 
of  an  integral  equation  and  then  develop  an  approximation  procedure  to  which 
they  attribute  good  accuracy.  On  the  other  hand  we  show  in  (1)  and  (2)  that 
their  approximation  procedure  is  a  non-uniform  one  and  hence  is  likely  to 
lead  to  unreliable  results.  Moreover  in  (1)  and  (2)  we  show  that  no  solution 
with  finite  (i.e.,  non-infinite)  crack  opening  displacements,  can  exist  for 
the  problems  posed  and  treated  numerically  by  Eringen  et  al  (3,4,5).  The 
way  this  was  done  was  to  reduce  the  problem  to  an  integral  equation  and  then 
use  analytic  continuation  arguments  to  show  that  no  solution  with  non-infinite 
displacements  could  exist.  For  example  in(l)  we  reduced  the  model  problem 
considered  in  (3)  to  the  integral  equation 

L 

1=  f  K(x-5)  f(5)  «,  | x |  <_  L  (2.10) 

(for  convenience  we  have  replaced  a  constant  by  unity)  f(£)  is  a  dislocation 
density  (for  laplaces  equation)  so  that  the  integrability  of  f(£)  is  necessary 
for  non-infinite  crack  opening  displacements.  The  problem  is  to  find  integrable 
f (£ )  satisfying  (2.10)  for  a  specified  kernel  K  resulting  from  a  given  non-local 
modulus.  The  reader  is  referred  to  (1)  and  (2)  for  more  details,  here  we  con¬ 
sider  again  the  simple  example  given  by  (2.8).  In  this  case 

K(x)  °  2  X  2  (2.11) 

x  +  a 

Note  that  when  a  equals  zero  (2.10)  is  a  singular  integral  equation  which  has 
a  well  known  solution  and  is  similar  to  that  derived  by  representing  a  crack 
by  a  distribution  of  dislocations.  However  for  a  ^  0  the  following  non¬ 
existence  proof  can  be  devised. 

Non-Existence  Proof 

For  complex  z,  define 
L 

♦(2)  =  ^  K(z  -  £)  f(£)  d£  (2.12) 

-L 

This  makes  sense  in  the  strip  |lm  zj  <  a  and  is  analytic  in  z.  Furthermore, 

<f>(z)  =  1  for  z  =  x,  |xj  _<  L  hence  <j>(z)  r  1  in  the  whole  region  |  Im  z  |  <  a,  and 

in  particular  for  z  =  x  (real)  for  all  x.  (Here  the  theorem  used  is  that  if 
ip(z)  is  analytic  in  region  D,  and  i/i(z)  =  0  at  a  sequence  of  points  z  =  zn  with 

zn  -*■  z*  in  D  then  i|i(z)  =  0  in  D.)  However,  for  real  x  >  L,  and  |£  |  _<  1 
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I«*  -MlirHr-rrT  <2a3) 

SO 

,  L 

|$(x)|  i  ^  f(°  dC 

This  gives  <j>(x)  +  0  as  x  +  >  contradicting  $(x)  i  1. 

This  rigorises  the  non-existence  argument  given  in  (1) .  In  the  next 
section  we  consider  the  numerical  solution  of  line  crack  problems  with  moduli 
of  the  form  given  by  (2.9),  i.e.,  where  there  is  an  'instantaneous'  elastic 
component  to  the  modulus. 


3  NUMERICAL  SOLUTION  OF  LINE  CRACK  PROBLEMS  IN  NON-LOCAL  MEDIA  WITH  MODULI  OF 
THE  FORM  a(jx|)  =  b  6(|x|)  +  (1  -  b)  oti(|x|) 


In  this  section  we  consider  a  line 
shear  where  the  non-local  modulus  is  of 
there  is  only  one  non-zero  displacement 

t13,l  +  t23, 2  =  0 


crack  problem  of  length  2  in  anti-plane 
the  form  given  in  (2.9).  In  this  case 
component  u3,  equation  (2.1)  reduced  to 

(3.1) 


and 


(2.2) 


t 


k3 


reduces  to 


S 


x’|)  u3>k  dS(x') 


(3.2) 


with  x  =  (xi,x2)  and  the  integral  in  (3.2)  being  over  the  whole  (x^,X2)  space. 
We  use  the  usual  line  crack  formulation  and  consider  the  stress  boundary 
condition  on  the  crack  line  as 

t23(x]LP)  =  -tQ(constant) ,  (x-J  <1  (3.3) 

Because  of  the  normalizing  condition  (2.5)  on  the  modulus  u'  it  is  clear  that 
u3,2  =  co  in  (3.2)  corresponds  to  t23  =  tQ  so  a  uniform  stress  could  be  added 
to  the  problem  giving  a  stress  free  crack.  A  formulation  like  this  is  common 
to  the  cited  papers.  For  the  numerical  analysis  given  here,  we  consider  a 
specific  form  for  y'(|xj),  this  is  y'(|xj)  -  pa(|x|)  where 

a(|x|)  =  b  6(|x|)  +li2^—  «(|feM|)  (3-4) 

Note  that  iot(lxl)  dS(x)  =  1  and  6  is  the  delta  function.  Defining 
S  v' 

°k3  =  u  u3,k  (3*5) 


arguments  are  given  in  (3),  (4)  and  (5)  to  show  that  equation  (3.1)  implies 
the  equation 
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I 


°k3.k  ‘  0 


(3.6) 


Using  screw  dislocation  solutions  of  (3.6)  with  (3.5)  and  thus  automatically 
satisfying  the  condition  of  continuity  of  displacement  outside  the  crack 
the  following  integral  equation  is  derived  in  (2). 


t23(Xl,0)  =  \  f3(ri)  k(X;l  -  n)  dn 


(3.7) 


where  f3(n)  is  the  unknown  dislocation  density,  f 3( n)  dn  ■  0  is  a  necessary 
condition  for  the  crack  to  close  at  xi  *  ±1  and 


7  1  r  a(!xil»lx2l)  (xl  +  t)  dxl 

-  2  i  S  ,  i  .'2,12 - - 


/  1  _\2  .  12 
(X;L  +  t)  +  x2 


Using  (3.4)  equation  (3.7)  can  be  written  as 


(3.8) 


t23^Xl*°) 


1  f3(n)dn 


j  -y- +  (1  -  b)  ]  f3(n)  K2(Xl  -  n)  dn  (3.9) 


The  kernel  comes  from  using  the  second  part  of  (3.4)  in  (3.8).  We  have 
solved  the  singular  integral  equation  (3.9)  by  turning  it  into  a  Fredholm 
equation  in  the  manner  described  in  (8)  and  then  solving  this  equation 
numerically  (the  derivation  is  too  lengthy  to  include  here). 


Writing 

f2(fi)  =  -N(n)/(1  -  n2)'2 


(3.10) 


we  have  calculated  N(l)  for  various  t  and  b  values.  Note  that  the  crack  2  i 
opening  displacement  at  the  tip  X^  =  1  will  be  proportional  to  N(l)(l  -  x^)^ 
so  the  square  root  singularity  common  to  the  classical  elasticity  problem 
occurs.  The  calculation  is  a  sensitive  one  especially  for  small  values  of 
b  so  results  are  shown  for  increasing  values  of  M  the  number  of  unknowns  in 
the  discretization  of  the  Fredholm  integral  equation.  When  b  =  1,  the 
elastic  case,  the  exact  result  is  N(l)  =  .31831.  Now  according  to  the  results 
derived  in  (2)  when  t  tends  to  zero  the  displacement  at  the  crack  tip  should 
have  N(l)  =  (.31831)  b-^,  i.e.,  the  elastic  crack  tip  displacement  multiplied 
by  b-'*.  As  can  be  seen  in  the  table  the  agreement  with  this  limiting  result 
is  very  good. 


l  -  .5 

b  -  .9 

.6 

.7 

.6 

.5 

.4 

M  -  40 

.34142 

.36866 

.40143 

.44192 

.49385 

.56406 

M  -  60 

.34169 

.36923 

.40232 

.44307 

.49500 

.56433 

M  -  100 

.34186 

.36966 

.40315 

.44458 

.49768 

.56926 

l  -  .2 

M  -  40 

.33661 

.35782 

.38284 

.41292 

.4500 

.4973 

M  -  60 

.33730 

.35938 

.38551 

.41716 

.4574 

.5085 

M  -  100 

.33794 

.36091 

.38832 

.42181 

.4641 

.5198 

t  -  .1 

M  •  50 

.33348 

.35074 

.3706 

.3938 

.4215 

.4552 

M  »  60 

.33435 

.3528 

.3743 

.4001 

.4328 

.4596 

M  •  100 

.33549 

.3554 

.3788 

.4070 

.4416 

.4857 

( .31831) 

.33553 

.3559 

.3805 

.4109 

.4502 

.5033 

TABLE  1  Values  of  N(l) 


4  CONCLUDING  REMARKS 

In  the  earlier  sections  and  in  references  (1)  and  (2)  we  have  specifi¬ 
cally  addressed  the  line  crack  problem  and  what  certain  non-local  (continuum) 
calculations  have  to  say  about  it.  We  hope  to  have  supplied  sufficient 
argument  to  suggest  that  the  claims  in  (6)  are  unsubstantiated.  Whether 
some  modification  of  the  boundary  conditions,  etc.,  may  resolve  the  dif¬ 
ficulties  we  have  referred  to,  is  at  present  unproven.  Note,  however  that 
if  one  were  satisfied  with  a  singular  stress  at  a  line  crack  tip  then  one 
could  still  have  continua  with  some  non-local  behavior  such  as  considered 
in  Section  3.  Finally,  we  mention  more  abstract  approaches  such  as  those 
of  Rogula  (9).  In  principle  we  believe  they  should  confirm  the  spirit  of 
our  conclusions. 
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STATUS  OF  BRAZILIAN  DISK  TEST  INVESTIGATION 

The  preliminary  investigation  of  the  cracked  Brazilian  disk  test  was 
completed.  The  results  of  this  work  are  summarized  in  the  attached  manu¬ 
script  which  has  been  accepted  for  publication  in  the  International  Journal 
of  Fracture.  Further  analytical  and  experimental  work  on  the  test  geometry 
are  warranted  to  assess  the  applicability  of  the  specimen  for  evaluation  of 
mixed  mode  fracture. 

In  addition,  J.  Avila  has  completed  an  M.S.  thesis  on  the  rod  pull  out 
problem.  The  major  finding  of  this  work  was  that  debonding  of  the  rod  from 
the  surface  was  favored  for  most  initial  flaw  configurations.  The  work  also 
presented  a  simple  scheme  for  obtaining  energy  release  rates  using  only  nu¬ 
merical  data  from  an  unflawed  specimen.  This  obviates  the  need  for  numerous 
computer  solutions  of  differing  flaw  lengths  or  the  use  of  complicated  data 
reduction  schemes  to  obtain  energy  release  data.  The  findings  of  this  work 
are  currently  being  summarized  in  a  manuscript  to  be  submitted  for  publication. 

The  parametric  study  of  this  problem  is  being  pursued  by  M.  Klienowsky. 
Some  numerical  difficulties  have  been  encountered  in  the  finite  element  scheme 
for  certain  of  the  geometries.  These  problems  are  currently  being  studied 
and  should  be  resolved  shortly.  The  results  of  this  study  should  be  available 
within  the  year. 
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1.  INTRODUCTION 


Recently,  a  series  of  papers  by  Eftis  et  al.  [1,2,3]  has  re-examined  the 
problem  of  an  inclined  crack  in  a  plate  under  remote  biaxial  loading.  Their 
purpose  was  to  demonstrate  the  effect  of  remote  stresses  and  load  biaxiality 
on  Kj  and  K^.^.  as  well  as  the  effect  on  various  fracture  criteria.  A  critical 
review  of  the  various  fracture  criteria  has  been  given  by  Swedlow  [4].  Here 
it  is  pointed  out  that  the  existing  test  data  do  not  allow  any  selection  to 
be  made  between  the  various  competing  criteria,  and  it  is  recommended  that 
alternative  tests  be  considered  as  well  as  an  expansion  of  the  data  base  for 
combined  mode  fracture. 

It  is  in  this  later  area  where  the  current  specimen  would  seem  to  be 
most  useful.  The  specimen  geometry  is  shown  in  Figure  1  and  is  a  variant 
of  the  early  work  of  Libatskii  and  Kovchik  [5]  and  Yarema  and  Krestin  [6]. 

The  specimen  of  Figure  1  has  been  investigated  by  Ajawi  and  Sato  [7]  and 
by  Sanchez  [8]  as  an  alternative  to  the  compact  shear  specimen  of  Jones  and 
Chisholm  [9,  10].  In  addition  to  the  determination  of  K  ,  Awa j i  and  Sato  [7] 

1  lw 

also  examined  the  possibility  of  the  existence  of  a  fracture  criteria  in  the 
form 


K  u  K  u 

(ith  +  (r1-) 

KIC  KIIC 


(1.1) 


In  contrast  to  the  numerical  analyses  of  [7,8],  we  present  explicit 
formulae  for  Kx  and  K^j.  valid  for  any  angle  of  crack  orientation  in  the 
cracked  Brazilian  disk  test,  Eq.  (3.1).  Also,  an  approximate  analysis  valid 
for  short  cracks  is  given,  and  the  effect  of  friction  is  included  when  crack 
closure  is  presented.  Finally,  experimental  results  are  presented  which 
demonstrate  the  feasibility  of  producing  short  cracks  in  moderate  size  speci¬ 
mens,  thus  indicating  the  usefulness  of  the  approximate  analysis  and  of  the 
specimen  geometry  in  expanding  and  illucidating  the  problem  of  combined  mode 
fracture.  A7 
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2.  ANALYSIS 


To  give  a  complete  numerical  treatment  of  the  problem,  we  follow  a 
method  outlined  by  Atkinson  [11].  The  crack  is  represented  by  a  continuous 
distribution  of  edge  dislocations,  simulating  discontinuities  in  the  normal 
and  tangential  components  of  displacement.  Because  the  crack  lies  on  a 
diameter  of  the  disc,  there  is  no  interaction  between  dislocations  repre¬ 
senting  discontinuities  in  the  normal  displacement  and  those  representing 
discontinuities  in  the  tangential  displacement.  The  tractions  on  the  axis 
y=o  due  to  the  unknown  distribution  of  dislocations  can  be  written  in  the  form 


The  precise  form  for  H^(c,  x)  can  be  deduced  from  [12]  where  the  case  of  a 
circular  inclusion  was  considered.  Note,  however,  a  misprint  in  equation 

(2.2.2)  of  [12].  The  factor  in  the  last  term  in  the  curly  brackets  should 

•> 

2  X*”  2G 

be  —  not  ^  .  In  (2.1),  D  =  — ^  +  ^ ,  G  is  the  shear  modulus,  and  <  =  3-4v 
(plane  strain)  or  <  =  ^  +  ^  (plane  stress)  v  being  Poisson's  ratio.  For 
convenience,  we  have  taken  the  Burger's  vector,  to  be  unity.  The 
corresponding  opening  of  the  crack  faces  can  be  written 

(l 

Aui  =  l  Nt(c)dc  i  =  1,2  (2.2) 

where  (u^,u0)  are  the  (x,y)  components  of  the  displacement,  and  is  the 

usual  stress  tensor.  The  subsidiary  conditions 

£ 

N  (c )dc  =0  i  *  1,2  (2.3) 

-l  1 

ensure  that  the  crack  closes  at  the  ends  (+  £,o). 


:  Vc>dc  Ni(c)Hi(c,x)  dc  i  =  1,2 

-l  X_C  J  -t 
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The  distributions  N^(c)  are  to  be  determined  from  the  condition  that  the 
crack  be  stress  free  and  that  it  closes  at  the  ends  (condition  (2.3)).  The 
first  of  these  conditions  is  met  by  replacing  a  ^  in  (2.1)  by  minus  the  stress 
induced  on  the  crack  line  (|x]“£,  y«o)  by  the  applied  stresses  set  up  in  the 
disc  in  the  absence  of  the  crack.  Complete  expressions  for  these  stresses  are 
given  in  Appendix  1.  For  the  moment,  we  assume  the  stress  is  represented  by 


ai2(x)  *  -o  U) 


i  =  1,2 


given  on  y*o,  |x|<£.  Equation  (2.1)  then  becomes  an  integral  equation  for 
the  unknown  densities  N^(c). 

Two  methods  of  treating  this  equation  numerically  are  given  in  [11]. 
One  of  these  methods  consists  of  solving  (2.1)  directly  by  expanding  (2.1) 
in  terms  of  Chebyschev  polynomials.  The  other  method  consists  of  first 


transforming  (2.1)  into  a  Fredholm  equation  and  then  solving  it  numerically. 
We  use  the  latter  method  here  since  it  is  possible  to  express  the  kernal  in 
closed  form. 

Using  (2.4)  in  (2.1)  the  resulting  Fredholm  equation  can  be  written 


(see  Atkinson  [11]) 


Nu  (n) 


n2(l-nV5 


!  yy  (1— xjv  dXi  + 

-1  X1  -n 


Nil^cl^  Fi(c1’r,)dc1 


i  =  1,2 


where  c  »  ic ^  ,  a  ■  <ia^  ,  x  *  £x^  , 

F1<c1,n)-f1  VfVfV  (1~*y  d. 


and  N^(r))  ■  N^(£r)).  Making  these  coordinate  transformations  in  (2.2),  gives 


Au1(x,o)  -  i 


Nil(cl)dcl 
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We  expect  N^Cn)  to  have  inverse  square  root  singularities  at  n  =*  +  1, 
so  we  write 


2  -^5 

Nil  (n)  “  Nio  (n)  (1_n  ) 


(2.7) 


Now  make  the  substitutions  n  ■  cos  a  ,  c^  -  cos  ip  and  x^  *  cos  ij^  so  that  equations 


(2.5)  become 


N.  (cosa)  ■  — -r 

io  n2 


f(cosa)  +1  N.  (cosi|>)  F  (cosi^,  cosa)  d ^ 
I  io  i 


(2.8) 


for  0  <_  a  <_  it.  f(cosa)  is  to  be  derived  from  the  first  integral  on  the  right 
of  equation  (2.5).  Note  that  (cos  4>,  cos  a)  can  be  expressed  in  closed 
form  by  completing  the  integral  following  Equation  (2.5). 


The  next  step  in  the  solution  of  the  problem  is  to  replace  the  integral  in 
(2.8)  as  a  sum  and  solve  the  resulting  system  of  simultaneous  equations.  Once 
this  is  done,  substituting  into  (2.6)  shows  that 

Au^x.o)  -  (2l)h  N.o(l)  Oe-x)**  (2.9) 

as  x  ->  .  Furthermore,  from  (2.2),  the  crack  opening  displacements  along 

the  crack  flanks  can  be  computed  as 


Au^(x  ,o) 


where 


(x/D 


N.  (cost(i)dtf> 
io 


(2.10) 
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3.  RESULTS 


The  analysis  of  section  2  may  be  used  to  determine  the  stress  intensity 
factors  associated  with  cracks  of  varying  length  and  orientation.  Such  a 
determination  requires  the  evaluation  of  f (cos  41)  in  Equation  (2.8). 
f(cos  b)  may  be  evaluated  through  an  expansion  of  o^(x)  in  terms  of 
Chebyschev  polynomials  and  then  numerically  integrating  to  determine 
f(cos  $).  This  procedure  is  numerically  sensitive.  Thus  an  alternate 
method  of  evaluating  f (cos  <p)  was  used.  Each  term  in  the  approximate 
stress  distributions  given  in  Appendix  1  may  be  independently  used  to 
solve  Equation  (2.8).  The  resulting  N^q(1)'s  yield  series  approximations 
to  the  stress  intensity  factors  in  the  form 

n  t  2i-2 
N  =  1  T  (V1  A  (e) 

1  i=l  1  a  i 

(3.1) 

n  n  r\  .  r* 

Ntt  =  2  sin  26  I  S.  (-)  B.  (0) 
ii  .  ,  l  a  1 

i=l 


where  Nj  = 


0  >'  7T-C 

o 


,  N 


'II 


II 


a  /rr 
o 


are  normalized  stress  intensities  and 


a  =  —  .  Tlie  T.'s  and  S.'s  in  (3.1)  are  numerical  factors  determined  through 
o  na  i  i  6 

the  solution  of  (2.8).  The  first  five  values  of  T^  and  and  the  corresponding 
A^(6)'s  and  B^(0)'s  are  given  in  Tables  1  and  2  respectively. 


The  normalized  uncracked  disk  stress  distributions,  cr 


00 


a90  ,  ap0 

—  *nd  °o0  =  ~ 
°0  °0 


for  —  =  0.6  are  shown  in  Figure  2.  As  can  be  seen,  the  two  term  approximation 
a 

yields  a  stress  distribution  greatly  differing  from  the 
analytical  solution.  However,  the  approximation  improves 
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for  small  Figure  3. 


The  stress  distributions  of  Figure  2 


produce  only  small  differences  in  stress  intensity  factors  as  evidenced 

in  Figure  4.  Here  and  for  the  two  and  five  term  approximations  for 

l 

—  =  0.5  are  compared  with  the  numerical  results  of  Awaji  and  Sato  [7] .  For 
the  five  term  approximations,  no  discernible  differences  can  be  detected  in 
the  results,  while  the  two  term  approximation  shows  differences  of  less  than 
ten  percent.  This  small  difference  for  the  two  term  approximation  is  due  to 
the  dominance  of  the  first  term  of  the  expansion.  In  addition.  Figure  4 
also  contains  a  series  of  finite  element  calculations  due  to  Sanchez  [8). 

These  results  were  obtained  using  a  special  crack  tip  element  of  Tong  e£  al . 

[13]  in  a  mesh  of  constant  strain  elements.  The  computations  are  seen  to  be 
in  agreement  with  the  other  numerical  results  for  the  coarse  mesh  of  Figure  5. 

Figure  4  shows  that  as  the  crack  changes  its  orientation  from  the  vertical 
to  the  horizontal  the  normalized  stress  intensity  factor  changes  from 
positive  to  negative  indicating  crack  closure.  At  the  point  of  crack  closure, 
the  model  of  the  crack  must  be  altered.  This  point  was  not  considered  in  [7] 
and  in  general  has  not  beer  incorporated  in  problems  of  mixed  mode 
fracture  [1-3].  Condition  (2.2)  must  be  modified  so  that 

Au2  =  0  (3.2) 

for  £*  <  | x |  <  i  where  t*  is  the  position  at  which  the  crack  first  closes  and 

>  >  * 

0  _  i*  t  on  y  =  0.  Now  if  Aug  is  zero,  an  additional  traction,  is 

developed  along  the  crack  surface  due  to  friction 

\n  (x)  1  y  °22  (x)  (3.3) 

In  the  following  section,  an  analysis  for  short  cracks  is  presented  and  the 
effect  of  friction  is  incorporated. 
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4.  SMALL  CRACK  APPROXIMATION 


If  the  crack  length  is  small  relative  to  the  radius,  then  the  second 
term  on  the  right  of  (2.1)  may  be  neglected,  and  the  problem  treated  as  a 
crack  in  an  infinite  media.  Equation  (2.1)  becomes 


i2_  f  *  Vc)  dc 

D  “  J  -t  x  -  c 


i  =  1,  2 


(4.1) 


Now  from  Appendix  1,  the  normal  and  shear  stresses  along  the  crack  line  are 
approximately 


«  -h  1*1  +  *2  <I>21 


—  ,  s  2P  sin  20  r„  ,  „  ,x(2, 

a12  (x)  — 77 —  [B1  +  B2  (7>  1 


(4.2) 


since  x  =  p  on  y  =  0.  Now,  022(x)  anc*  c>i2(x)  are  even  functions  of  x,  and  the 
normalized  stress  intensity  factors  may  be  written  as 
1 


Nt  = 

I  or 

o 


t  c?2 2  (x)  dx 


'll 


a  it 
o 


-t  a2  -  x2)4 


t  0^2  (x)  dx 

■i  a2  -  x2)H 


(4.3) 


[Note,  N j  and  Njj  are  not  to  be  confused  with  the  dislocation  densities  Ni(c).] 
Substituting  (4.2)  into  (4.3)  and  integrating  yields  explicitly 


Ao  f  2 

NT  =  A.  +  ~  (— ) 


n 

Nn  Bi  +  T  (7} 


(4.4) 
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A  comparison  of  expressions  (4.4)  with  expressions  (3.1)  and  the  finite 

£ 

element  results  are  given  in  Figure  6  with  —  =  0.3.  The  approximations  of 

SL 

(4.4)  are  seen  to  be  quite  good. 

From  Equation  (4.4)^  we  may  obtain  an  explicit  equation  for  the  angle 
at  which  the  crack  just  begins  to  close 

Ao  /  9 

NT  =  A.  +  -r-  er  =  o  (4.5) 

I  12a 

£ 

This  expression  has  been  evaluated  for  fixed  values  of  — ,  and  the  results 
are  compared  with  those  of  [7,8]  in  Table  3.  For  angles  of  inclination 
greater  than  those  in  Table  3,  the  crack  closes  over  some  portion  of  its 

surface  and  a  frictional  stress  is  developed. 

£ 

If  in  (4.5)  9  is  fixed  in  deference  to  — ,  we  obtain  an  equation  for  the 
position,  £*,  at  which  the  crack  closes.  Now  equations  (4.4)  obtain  for 
p  <_  l*.  If  £  _>  £*,  we  have  in  addition  to  the  stress  intensity  factor  N ^ 
given  by  (4 . 4 ) 0  an  additional  contribution  due  to  the  frictional  stress. 

We  note  that  this  is  true  even  if  Nj  is  constrained  from  being  negative  by 
imposing  AU2  =  0  on  v  =  0  for  0  <_  l*  |x|  £.  Additionally,  this  is  true 

in  the  general  case  since  there  is  no  interation  between  the  normal  and 
shear  dislocations. 

To  evaluate  the  additional  contribution  to  N^j,  we  must  evaluate  the 
stresses  in  the  region  L*  <_  |x|.  The  stresses  ahead  of  the  crack  in  this 
small  crack  approximation  are  given  by 


°22  (x)  f  (  (c)  dc 

°  J  -t*  C  "  x 


(4.6) 


'M 

The  ^(c)  can  be  determined  from  the  first  term  on  the  right  hand  side  of 
Equation  (2.5)  provided  we  replace  t  by  (.*,  n  by  c  and  use  expression  (4.2)^ 
for  -o^(x)  with  the  condition  that  Au2  =  0  on  v  =  0  |x|  >  l *. 
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1 


I 


Equation  (4.6)  may  then  be  evaluated  to  give,  for  |x|  >  &* 


o22  (x) 


a  x 
o 


2  2  n*2 
A.  -  -=■  (x  -  £  Z) 

1  a2 


(4.7) 


and  the  friction  stress  within  the  contact  zone  is  given  by 


°12  =  -u  °22 


(4.8) 


p  being  the  coefficient  of  friction. 

So  in  addition  to  the  stress  intensity  factor  of  (4.4)2,  we  have 
from  (4.3)2 


*  =_2p 

II  a  it 


£  o22  (x)  dx 

,*  a2" -  x2)14 


(4.9) 


in  the  zone  l*  <_  |x|  <_  l.  Inserting  c>22  (x)  from  (4.7),  we  find 


*  "U  A2  2  *2 
11  ~  2 


a  - 1  n 


which  must  be  added  to  (4.4)2  to  account  for  the  friction  developed  within 


the  contact  zone. 


5.  EXPERIMENTAL  RESULTS 

Tests  were  conducted  on  a  commercial  grade  of  polymethyl  methacrylate 

6  2 

(PMMA)  with  an  elastic  modulus  of  0.5  x  10  16/in  and  a  tensile  strength 

3  2 

of  8  x  10  lb/in  .  The  specimens  were  1.5  in.  in  diameter  and  approximately 
0.125  in.  thick.  Details  of  the  specimen  preparation  and  testing  procedures 
may  be  found  in  [8J.  In  addition,  a  series  of  three  point  bend  tests  were  also 
conducted.  To  avoid  errors  which  may  be  produced  in  specimen  preparation,  the 
beam  specimens  were  cut  from  3.0  inch  disks  which  were  prepared  in  a  fashion 
identical  with  the  1.5  in  disks.  The  results  for  the  opening  mode  failure  are 
given  in  Fig.  7.  In  Fig.  7,  Kj  has  been  normalized  by  using  the  average  Kjc 
values  from  both  the  bend  specimen  and  the  disk  test.  As  is  evident,  the 
experimental  results  follow  the  analytical  curve  of  Libatskii  and  Kovchik  [5] 

quite  closely  with  maximum  deviation  at  short  crack  lengths. 

Figure  8  presents  results  fcr  the  mode  II  fracture  of  PMMA.  The  results 
here  show  greater  scatter  than  in  mode  I.  This  is  most  probably  due  to  the 
precision  required  in  orienting  the  specimen.  Nevertheless,  it  is  felt  that 
the  above  results  demonstrate  the  feasibility  of  using  the  disk  test  as  a 
reliable  method  for  obtaining  K  and  K  fracture  parameters.  Additionally, 

llw 

the  testing  also  demonstrates  that  the  production  of  specimens  of  modest  size 

l 

with  small  —  ratios  is  feasible.  This  coupled  with  analysis  of  Section  4  will 
allow  a  detailed  study  of  the  effects  of  combined  mode  fracture  as  well  as  the 
role  of  friction  in  the  failure  event. 
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6.  CONCLUSIONS 


The  foregoing  has  presented  an  analysis  of  the  cracked  Brazilian  disk 
test.  The  results  of  a  numerical  investigation  yielded  explicit  formulae 
for  the  mixed  mode  stress  intensity  factors  valid  for  any  crack  orientation. 

In  addition,  it  was  noted  that  for  some  angle  the  opening  mode  stress  intensity 
factor  became  zero,  and  beyond  this  angle,  the  crack  was  closed.  An  approxi¬ 
mate  analysis  was  developed  for  short  cracks  and  a  scheme  for  the  inclusion 
of  friction  was  presented.  Experiments  were  conducted  on  PMMA  specimens,  and 
it  was  demonstrated  that  modest  size  specimens  could  be  produced  with  short 
crack  lengths  which  yield  reliable  results.  Thus,  it  is  possible  to  analvze 
and  investigate  the  problem  of  combined  mode  fracture  using  the  cracked  Brazil 
disk  test. 
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APPENDIX 


Following  [14,15],  the  stresses  relative  to  a  cylindrical  polar 


coordinate  system  p,  0,  located  at  the  center  of  the  disk  are 


j  2P 

pp  =  rr 
ira 


.  (1-  —  cos0)  (cos0-  — )2 

_ a _ a 

(1+  -  2  —  cos0)2 

z  a 


(1+  —  COS0)  (COS0+  —  )2 
a  a 


(1+  *=•  +  2  ^  cos0)  J 
2  a 

a 


00 


2P 

ira 


2 

(1-  —  cos0)  sin  0 
_ a _ 

2 

(1+  -  2  ^  cos©)2 

Z  a 
a 


(1+  —  cos0 )  sin20 
_ a _ 

(1+  4  +  2  J  cos  )2 
a 


(A.l) 


p0  ira 


(1-  —  cos0)  (cos0-  — )sin0  (1+  —  cos0)  (cos0+  — )  sin0 

a _ a _ _ a _ a _ 


L  (1+  ^  -  2  fi-  COS0)‘ 
z  a 
a 


(1+  ^  +  2  cos0 )  2 
z  a 
a 


These  may  be  expanded  in  a  power  series  in  —  to  give  for  o  and  v  _ 

a  yy  pu 


n 


°ee  '  h  ill  V6>  'I*2”'2 


(A.  2) 


P0 


2Psin20  "  B  L  2n-2 
ira  i=  1  l  a 


where  A^(0)  to  A^(0)  and  B^(0)  to  B^(0)  are  given  in  Table  2. 


These  expressions  are  due  to  Prof.  R.J.  Erdlac,  Dept,  of  Mechanical 
Engineering,  University  of  Pittsburgh 
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Angles  for  Pure  Mode  II  Behavior  for  Various  Crack 
Lengths 


TABLE  1 


S1 

T2 

S2 

T3 

S3 

T4 

S4 

1.014998 
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0.376991 

0.376991 
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1.257488 
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0.318086 

0.316245 

0.325033 

0.320834 

0.334831 

0.327411 

0.347941 

0.336447 

0.365559 

0.349219 


TABLE  2 


2  2 
8s  (l-4c  ) 

2  2  4 

-4s  (3-36c  +48c  ) 

-16s2  (-1+24c2-80c4+64c6) 

-20s2  (1-40c2+240c4-448c6+256c8) 


-5+8c 

-3+8(1-2c2)  (2-3c2) 

3+16(l-2c2)  -  12  (l-2c2)2  -  32  (l-2c2)3 
5-16(l-2c2)  -  60  (l-2c2)2  +  32  (l-2c2)3  +  80  (l-2c2)4 


s  =*  sin0;  c  =  cos0;  cos20  “AC¬ 
TABLE  3 


Sanchez  [8] 

Awa j i  and 

Saco  [7] 

27.7 

27.2 

25.2 
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23.2 

22.9 
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20.1 
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Figure  1:  Geometry  of  the  cracked  Brazilian  disk  test. 

Figure  2:  Stress  distributions  in  the  uncracked  disk  for  -^  =  0.6. 

Figure  3:  Stress  distributions  in  the  uncracked  disk  for  ~  =  0.3. 

I 

Figure  4:  Stress  intensity  factors  for  —=0.5 
a)  Nr  b)  Nir 

Figure  5:  Finite  element  mesh  for  stress  intensity  factor  determination. 

Figure  6:  Stress  intensity  factors  for  the  small  crack  approximation 
a)  Nj,  b)  Njj. 

Figure  7:  Critical  Mode  I  stress  intensity  factor  for  PMMA. 

Figure  8:  Critical  Mode  II  stress  intensity  factor  for  PMMA. 
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EFFECTS  OF  CHANGE  IN  CRACK  DIAMETER  ON 
CRITICAL  LOAD-PRELIMINARY  STUDY 


By 


Shan  Somayaji 


Introduction 

Williams  showed  that  from  the  viewpoint  of  continuum 
mechanics,  and  particularly  the  energy  concept  of  fracture 
that  adhesive  and  cohesive  failure  are  similar  (1),  Experi¬ 
mental  investigations  such  as  blister  tests  and  modified 
blister  tests  were  carried  out  in  the  last  decade,  the  results 
of  which  tend  to  prove  the  above  hypotheses,  Malyshev  and 
Salganik  (3)  used  a  concentrated  load  to  tear  off  a  thin  plate 
bonded  flatwise  to  a  rigid  base  and  due  to  the  axial  symmetry, 
the  crack  was  assumed  to  be  penny-shaped.  They  tested  a 
plexiglass-steel  combination,  loaded  at  the  center  to  a  point 
where  the  load  stopped  changing.  The  fracture  energy  so 
determined  was  found  to  be  in  close  agreement  with  that 
determined  by  other  types  of  tests.  Jones  and  Williams  used 
a  method  similar  to  that  of  Malyshev,  et  al.  ,  except  that  the 
central  point  load  was  replaced  by  a  uniform  pressure  (2). 

They  found  that  the  pressure  loading  was  simpler  than  the 
central  load,  and  resulted  in  a  smaller  deflection.  The  test 
configuration  consisted  of  a  glass  disk  with  a  central  hole  bonded 
to  a  thin  sheet  of  polyurethene  rubber.  Simultaneous  measurements 
were  made  of  the  diameter  of  the  unbonded  area  at  that  critical 
imposed  pressure  at  which  the  unbond  radius  would  increase. 
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In  the  Malyshev  -  Salganik  investigations  the  debond  radius 
was  assumed  to  occur  at  a  critical  load  but,  with  neutral  stability, 
to  increase  as  the  load  remained  constant.  Substantial  test  data 
are  not  available  to  verify  this  assumption  and  it  is  probable  that 
the  effects  of  non-linearity  in  material  properties  and  other  sources 
of  energy  dissipation  may  affect  the  critical  load  to  induce  debonding 
and  whether  it  would  then  remain  constant  or  not.  The  objective 
of  this  study  is  to  make  a  preliminary  estimate  of  such  effects 
upon  the  critical  load  and  on  the  adhesive  fracture  energy.  An 
experimental  model  similar  to  the  one  adopted  by  Malyshev,  et  al 
(3)  was  therefore  used  to  test  a  plexiglass-steel  combination.  The 
model  was  equipped  with  a  dial-gage  apparatus  to  measure  the 
central  deflections.  A  number  of  load-deflection  diagrams  were 
plotted  for  each  specimen  at  different  initial  diameters  from  which 
a  preliminary  estimation  of  variation  in  the  critical  load  was  made. 

The  details  of  the  test  procedure  and  test  results  are  given  in  this 
report. 

Test  Details 

The  specimen  details  and  the  test  set-up  are  shown  in  Fig.  1. 

The  specimens  were  made  up  of  1/8  in.  and  1/16  in.  plexiglass  plates 
of  2.25  to  2.4  in.  diameter,  bonded  to  3  inch  diameter  steel  plates 
using  "Epox  Cement"  manufactured  by  G.  C.  Electronics,  Rockford,  II. 
The  load  was  applied  to  the  plexiglass  plate  by  a  plunger  passing  through 
a  central  hole  in  the  steel  plate  as  shown  in  Fig.  1.  A  small  area  at 
the  center  of  the  specimen,  of  about  0.  5  in.  diameter,  was  covered 
with  a  circular  piece  of  scotch  tape  to  prevent  clogging  the  central  hole 
during  the  initial  bonding  of  the  specimen.  Moreover,  this  area  at  the 
center  of  the  specimen  where  the  st  eel  and  plexiglass  are  unbonded, 
acts  like  a  crack  initiator.  Proper  precautions  were  taken  such  as 
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cleaning  the  surfaces  with  alcohol,  uniform  mixing  of  the  epoxy- 
cement,  constant  curing  time,  etc.  The  manufacturer's  suggestions 
were  followed  for  the  mixing  operation.  The  specimens  were  lab-cured 
for  24  hours  after  which  they  were  tested  using  a  Baldwin  Universal 
Testing  machine. 

The  loads  and  the  central  deflections  were  recorded  during 
the  test.  A  special  dial-gage  apparatus,  shown  in  Fig.  2,  was  used 
to  record  the  central  (under  the  load)  deflection.  (The  recorded  read¬ 
ings  were  converted  to  deflection  readings  using  the  proportionality 
constants:  deflection-gage  readings  x  5).  For  a  few  specimens  the 
loads  were  applied  continuously  until  complete  debonding  occured  but 
for  the  majority  of  the  specimens  a  number  of  load-deflection-crack 
diameter  readings  were  taken  prior  to  complete  debonding.  To 
accomplish  this  the  test  was  stopped  and  the  specimen  was  unloaded 
when  the  load  remained  steady  or  started  to  decrease.  Indication  of 
a  decrease  in  the  loading  rate,  (or  when  the  load  seemed  to  remain 
steady)  was  observed.  The  diameter  of  the  actua±  debonded  area  as 
seen  through  the  plexiglass  plate  then  was  measured  using  a  scale. 

The  specimen  was  set  under  the  plunger  again  and  the  loading  process 
was  repeated.  The  specimen  was  again  unloaded  and  the  diameter 
of  debonded  area  was  recorded.  This  procedure  of  loading,  unloading, 
and  measurement  of  the  crack  extension  was  repeated  a  number  of  times 
prior  to  complete  debonding.  Tabulations  of  the  load-deflection-crack 
diameter  readings  for  most  of  the  specimens  that  were  tested  in  the 
program  and  the  load-deflection  curves  are  in  Ref.  4.  For  the  few 
specimens  for  which  the  loads  were  applied  continuously  until  the 
final  break-off,  typical  load-deflection  diagrams  are  as  shown  in  Fig.  3. 
The  loads  could  be  seen  to  increase  linearly  to  a  peak  load,  followed 
by  a  sudden  drop  and  a  gradual  increase  combined  with  an  increase  in 
deflections.  Although  no  conclusive  statements  can  be  drawn  at  this 
stage  due  to  the  limited  number  of  tests  performed,  there  is  an 
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indication  that  the  extent  of  the  load-drop  may  be  influenced  by 
the  thickness  of  the  specimen.  It  was  also  found  that  the 
variations  in  the  initial  peak -load  between  specimens  of  the 
same  thicknesses  are  small. 

Fig.  4  shows  typical  load-deflection-crack  diameter 
diagrams.  The  initial  and  final  crack  diameters  are  indicated 
on  each  load-deflection  diagram,  the  details  of  which  are  given 
in  Ref.  4.  The  test  results  indicated  that  for  most  of  the 
specimens,  the  diameter  of  the  cracked  area  increased  with 
increase  in  load,  although  this  increase  (in  load)  was  quite 
small  for  a  few  specimens.  Due  to  the  preliminary  nature  of  the 
investigation  no  solid  conclusions  can  be  drawn  at  this  stage  but 
it  is  felt  that  the  shape  of  the  load  deflection  diagram  during  the 
crack  extension  stage  is  dependent  on  the  specimen  thickness. 

For  most  specimens  of  1/16  in.  thickness  the  final  break-off 
occurred  at  a  load  smaller  than  the  initial-diameter  peak  load, 
whereas,  for  specimens  of  1/8  in.  thickness,  the  final  break-off 
load  was  not  much  different  from  the  initial -diameter  peak  load. 

Fracture  Energy  Determination 

Malyshev,  et  al  (3)  determined  the  adhesive  fracture 
energy  for  centrally  loaded  bonded  specimens  as  follows: 

F or  a  circular  plate  loaded  at  the  center,  the  central 
deflection  w  is  given  by 

W  =  Pa^ 

16  TV  D  (1) 

where  P  is  the  central  load,  a  is  the  initial  radius  and  D  is  given  by: 

D  =  Eh^ 

12  (1  -  -O  (2) 

where  E  is  the  modulus  of  elasticity,  h  is  the  plate  thickness  and 
^  is  Poisson's  ratio. 


102 


For  an  increase  in  crack  radius  of  amount  da,  the  area  of  debond 
extension  is  approximately  equal  to  (2  fT  a)  da.  Using  the  energy 
balance  fracture  criteria 


(2  TV  a  da)  T  =  dU 

fi 

where  is  the  specific  adhesive  fracture  energy  and  U  is  the 
elastic  strain  energy. 

But,  U  is  given  by  one  half  the  applied  load  times  the 
equilibrium  displacement,  or 

u  -  1  f  p2j2  \ 

U  "  2  [  TSTTD-  / 

from  which  dU  can  be  calculated  and  inserted  into  (3)  to  give 

y  = _ P-L^ 

A  32 IT2D 

(Note  that  the  absenct  of  the  radius  in  (5)  indicates  neutral 
stability.  ) 


(3) 


(4) 


(5) 


From  this  preliminary  investigation,  the  critical  loads 
were  found  to  be  in  the  range  of  80  to  120  lbs  for  specimens  of  1/16  in. 
thickness  and  90  to  180  lbs  for  specimens  of  1/8  in.  thickness.  This 
shows  a  large  scatter  in  the  test  reulsts  for  specimens  of  the  same 
series,  of  the  order  of  about  100  percent.  But  for  a  given  specimen 
the  average  increase  in  the  critical  load  with  increase  in  crack 
diameter  of  about  100%  (0.  50  in.  -1,0  in.  )  was  found  to  be  about 
16.5  percent.  This,  when  compared  with  the  above  values  of  the 
experimental  scatter  in  the  critical  load,  can  be  regarded  as  small. 

Since  the  adhesive  fracture  energy  is  proportional  to  the 
square  of  the  critical  load,  a  small  change  in  the  critical  load  can 
strongly  affect  its  predicted  value.  But  considering  the  fact  that 
even  specimens  of  the  same  thickness  failed  at  critical  loads  that  were 
distinctively  different,  it  may  be  justified  to  neglect  the  effects  of  the 
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increase  in  critical  load  (with  increase  in  initial  crack  diameter) 
in  fracture  energy  determination.  However,  this  is  not  intended 
to  down-play  the  role  of  the  critical  load  in  fracture  energy 
determination  which  is  obvious  from  Eq.  5. 

A  revised  method  which  makes  use  of  direct  measure¬ 
ments  of  the  strain  energy  for  each  increase  in  crack  diameter 
is  suggested  in  Ref.  4.  Due  to  the  preliminary  nature  of  this 
investigation,  no  conclusive  statements  can  be  drawn  using  that 
method  at  this  stage.  Further  research  is  suggested. 

Discussions  and  Conclusions 

The  Malyshev-Salganik  model  to  determine  the  specific 
adhesive  fracture  energy  has  been  re-examined  and  verified  for  a 
plexiglass  to  steel  bond. 

The  preliminary  test  results  have  suggested  that  there  may 
be  a  load  with  change  in  crack  diameter  which  is  not  predicted  by 
the  elementary  strength  of  materials  analysis.  Further  investigations 
are  recommended  to  determine  the  effects  of  plate  thickness  and 
plate  diameter  on  the  adhesive  fracture  energy  and  load-deflection 
diagrams.  As  an  alternative,  direct  measurement  of  fracture 
energy  from  the  load-deflection  diagrams  is  also  recommended. 
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A  Bibliographic  Overview  of  Composite  Materials  with  Special 
Focus  upon  Recent  Publication  on  Fracture 


1. 0  Background 

A  retrospective  search  was  made  for  survey  or  research 
articles  (reports  and  journals)  covering  a  broad  aspect  of  the 
problem  rather  than  a  specific  area.  There  have  not  been  a 
great  many.  In  fact  much  of  the  early  material  appeared  as 
chapters  in  books. 

In  1965  the  National  Materials  Advisory  Board  (NMAB) 
published  a  "Report  of  the  Ad  Hoc  Committee  on  Interface 
Problems  in  Fibrous  Composites"  (1)  dealing  with  fiber  surface, 
the  process  of  composite  formation,  retention  of  properties 
in  service,  and  the  interface  during  composite  failure.  The 
objective  was  "to  display  potentially  critical  facets  of  the 
situation  in  order  to  identify  areas  in  which  additional  research 
is  needed  . .  . .  ".  A  second  report  came  from  NMAB  in  1965 
entitled  "Composites"  (2).  The  purpose  of  this  study  was  to 
examine  research  and  development  problems  which  limited  the 
further  advance  of  composites  technology  and  included  re¬ 
commendations  for  research  and  development  efforts.  In  1974 
the  NMAB  published  "Metal-Matrix  Composites:  Status  and 
Prospects"  (3),  a  report  of  the  Ad  Hoc  Committee  on  Metal-Matrix 
Composites.  The  stated  goal  of  this  study  was  to  determine  the 
role  of  metal-matrix  composites  particularly  in  DCD  applications. 
Specific  subjects  covered  were:  the  state-of-the-art;  the  advantages 
and  disadvantages  compared  to  monolithic -wrought  materials  and 
resin-matrix  composites;  the  potential  in  DCD.  NASA,  and 
commercial  applications;  fabrication  capability  and  the  status  of  the 
data  base;  cost  projections:  previous  service  experience  and  current 
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impediments  and  requirements  for  expediting  wider  use  of 
metal- matrix  composites. 

Also  in  1974  an  article  appeared  in  the  AIAA  JOURNAL 
entitled  "Composite  Material  Mechanics:  Structural  Mechanics" 
by  C.  W.  Bert  and  P.  H.  Francis  (4).  The  purpose  of  that  survey 
was  to  bring  together  and  review  some  of  the  principal  contributions 
in  structural  mechanics  to  the  field  of  composite  materials. 

Emphasis  was  given  to  the  macro -mechanical  structural  analysis 
of  various  structural  elements.  In  the  meantime,  in  1972. 

Fracture.  VII  appeared  with  the  chapter  on  "Fracture  Mechanics 
of  Composites"  by  H.  T.  Corten  (5). 

Continuing  the  concern  for  metal  matrix  composites 
applications  to  DCD  requirements,  the  Institute  for  Defense  Analyses 
Science  and  Technology  Division  published  the  proceedings  of  the 
First  MMC  Workshop  held  9-10  July  1975.  "Technological 
Development  of  Metal  Matrix  Composites  for  DOD  Application 
Requirements"  (6).  The  purpose  of  the  workshop  was  to  assess 
the  technology  of  metal  matrix  composites  for  DOD  application 
requirements  and  included  applications,  loads  and  environments, 
analyses  and  design  and  manufacturing  and  fabrication.  Also  in  1975 
Illinois  Institute  of  Technology  prepared  a  report  on  "Fracture  Tough¬ 
ness  of  Fiber  Composite  Materials"  (7)  for  the  Air  Force  Office  of 
Scientific  Research.  Fracture  toughness  test  methods  were 
evaluated  and  fracture  toughness  parameters  measured  for  random 
fiber  composites  having  epoxy  resin  matrices. 

"Fracture  Mechanics  of  Composite  Materials"  by  C.  K.  H. 
Dharan  (8)  published  in  the  Journal  of  Engineering  Materials  and 
Technology  in  July  1978  attempted  to  present  a  survey  of  the  field  in 
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its  current  state  and  to  demonstrate  the  approach  taken  by  various 
investigators.  "Research  Priorities  for  Advanced  Fibrous  Composites" 
a  report  by  K.  J.  Baumann  and  J.  L.  Swedlow  (9)  was  sponsored  by 
NASA/Lewis  and  recommends  research  priorities  in  advance 
fibrous  composites  for  the  aerospace  industry.  A  state-of-the-art 
review  is  presented  with  service  experience  and  technical  literature. 
The  bibliography  is  very  long  and  very  complete  and  precludes  the 
need  to  duplicate  the  effort.  It  is  also  very  current  and  probably 
includes  the  majority  of  the  composites  literature  through  1979 
featuring  in  particular  the  Journal  of  Composite  Materials.  The 
Journal  of  Composite  Materials  has  been  reviewed  for  1980  and  the 
relevant  articles  have  been  added  in  the  attached  bibliography. 

Inasmuch  as  the  primary  immediate  purpose  of  this  effort 
has  been  to  examine  tbe  fracture  of  composite  materials,  a 
selected  bibliography  of  relevant  citations  is  attached.  The  more 
general  problem  of  adhesive  fracture  can  of  course  be  viewed  as 
a  fracture  in  composite  media  and  for  this  reason  six  basic 
historical  references  to  the  bimaterial  problem  have  been  included 
in  the  collection  of  special  fracture  references.  References  have 
also  been  included  to  the  related  Fiber  Pull-Out  Problem  and  the 
Interfacial  Crack  Problem  bibliographies  prepared  for  the  last 
annual  report  (SETEC  SF  79-5,  October,  1979). 
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FOREWORD 


Each  of  the  Services  within  the  Department  of  Defense  has  established 
agencies  having  the  responsibility  of  advancing  the  state  of  the  art  in  those 
technologies  having  military  value.  Allocations  of  funds  to  specific  R  &  D 
programs  are  made  in  accordance  with  formally  adopted  "requirements"  with  due 
consideration  to  cost-benefit  ratio,  technical  feasibility,  etc. 

Before  research  in  a  given  area  starts,  technical  feasibility  is  often 
quite  uncertain,  but  is  viewed  optimistically  if  the  potential  benefits  are 
large.  If  after  a  number  of  years  of  effort  rate  of  progress  is  judged  to  be 
slackening  as  a  result  of  the  law  of  diminishing  returns,  the  priority  of  that 
field  of  endeavor  is  lowered,  and  the  corresponding  objective  may  even  be 
dropped  as  a  military  requirement.  In  that  case,  even  if  a  technical  break¬ 
through  is  subsequently  discovered,  no  funds  to  pursue  it  are  available  because 
no  allocation  was  made  in  the  budget.  When  this  happens  the  judgment  that 
future  progress  will  be  of  limited  value  becomes  a  self-fulfilling  prophecy. 

A  case  in  point  may  be  the  development  of  a  micromechanical  theory  for 
the  tensile  strength  of  fibrous  composites.  During  the  1960's  and  early  70's 
research  of  this  type  was  extensively  supported  by  the  three  services  (especially 
the  Air  Force),  NASA,  and  the  Atomic  Energy  Commission.  In  spite  of  this  a 
rational  failure  criterion  analogous  to  Griffith  crack  instability  theory 
failed  to  emerge.  Today  funding  in  this  field  is  at  a  very  low  ebb.  Even 
though  a  breakthrough  appears  to  be  on  the  horizon,  it  may  not  be  pursued 
because  it  was  not  anticipated  and  therefore  is  not  budgeted. 

According  to  an  ancient  Chinese  proverb  "It  is  better  to  light  a  candle 
than  curse  the  darkness".  Translated,  those  who  perceive  an  unsatisfactory 
situation  have  a  duty  to  call  it  to  the  attention  of  others  in  a  position  to 
do  something  about  it.  This  thought  was  the  motivation  of  the  following  paper. 
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INTRODUCTION 

Structural  efficiency  was  for  many  years  limited  by  an  empirical  law  which 
states  that  the  strength-to-weight  ratios  of  all  of  the  common  structural 
materials  tend  to  be  roughly  the  same.  This  was  also  true  of  stiffness-to-weight 
ratios.  A  few  decades  ago  it  began  to  be  recognized  that  certain  fibers,  such 
as  glass,  boron,  kevlar,  and  carbon  or  graphite,  could  be  prepared  with  many 
times  the  specific  strength  or  specific  stiffness  values  attainable  in  bulk 
materials.  As  a  result,  a  very  substantial  development  effort  was  undertaken, 
especially  within  the  military  services,  to  capitalize  on  these  properties. 

In  spite  of  the  extensive  use  now  being  made  of  fibrous  composites  there 
is  a  widespread  feeling  that  the  current  applications  fall  far  short  of  the 
ultimate  potential  of  such  materials  and  are  not  coinnensurate  with  the  very 
large  development  effort  carried  out  during  the  last  twenty  years.  There  are 
a  number  of  factors  contributing  to  this  unfortuante  situation.  Among  them 
are.  practical  matters  such  as  joining  difficulties  and  reduced  reliability 
sometimes  resulting  from  manufacturing  imperfections.  But  there  can  be  little 
doubt  that  one  of  the  most  important  causes  is  lack  of  confidence,  due  to  the 
fact  that  the  fundamentals  of  damage  accumulation  in  fibrous  composites  during 
the  loading  process  and  the  mechanics  of  ultimate  failure  are  as  yet  poorly 
understood. 

Since  the  basic  purpose  of  high  performance  composites  is  to  capitalize 
on  the  extraordinary  tensile  properties  of  the  fiber  employed,  the  first  and 
most  important  task  in  rectifying  this  deficiency  is  to  achieve  an  understanding 
of  the  tensile  strength  of  unidirectionally  reinforced  composites.  The  severity 
of  the  shortcomings  of  presently  available  theory  are  probably  known  only  to  a 
few  specialists  in  this  field.  Those  required  to  make  strength  calculations 
try  to  compensate  for  a  lack  of  basic  understanding  by  making  large  empirical 


corrections  to  the  theory  and  using  substantial  safety  margins. 


The  next  section  will  briefly  review  and  critique  published  theories 
for  the  tensile  strength  of  unidirectionally  reinforced  fibrous  composites. 
The  succeeding  section  will  outline  a  new  theoretical  approach  which  offers 
a  prospect  of  eliminating  empiricism  and  providing  a  sound  basis  for  greatly 
improved  design  estimates  of  composite  strength. 
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REVIEW  OF  PREVIOUS  WORK 


The  theory  for  damage  accumulation  and  fracture  of  a  "loose  bundle",  or 
collection  of  independent  fibers  of  equal  length  loaded  in  parallel,  was 
worked  out  decades  ago  by  Daniels  and  Coleman  (1-2).  The  subject  is 
by  now  quite  well  understood.  Such  is  not  the  case,  however,  with  fibrous 
composites.  In  many  composites  a  relatively  soft  matrix  is  used  to  bind  very 
strong  fibers  together.  The  matrix  does  not  appreciably  contribute  directly 
to  the  tensile  strength,  but  it  contributes  indirectly  by  transferring  stress 
between  fibers.  When  a  fiber  breaks  in  such  a  composite  its  load  is  trans¬ 
ferred  to  neighboring  fibers  by  the  matrix.  But  the  stress  tends  to  equalize 
with  distance  from  the  cross  section  at  which  the  break  occurred,  and  returns 
to  normal  in  a  distance  0.5  £  on  each  side  of  the  break.  This  distance  is 
variously  referred  to  as  "ineffective  length"  and  "stress  transfer  length". 

The  net  effect  is  to  localize  the  damage. 

Many  papers  containing  approximate  theoretical  solutions  for  the  tensile 
strength  of  such  a  composite  have  been  published.  Harter  (3)  has  reviewed  the 
work  of  the  major  contributors  to  this  field  through  1976.  In  the  present  paper 
the  treatment  of  prior  work  will  be  limited  to  a  very  brief  discussion  of  a 
representative  sample  of  these  contributions.  Suffice  it  to  say  at  this  point 
that  for  nearly  two  decades  theoretical  analyses  have  generally  been  based  on  a 
model  which  regards  a  composite  as  a  chain  of  very  short  bundles  which  we  shall 
refer  to  as  "minibundles". 

In  1962  Gucer  and  Gurland  (4)  proposed  a  model  for  a  particulate  composite 
consisting  of  many  slices  connected  in  series.  Rosen  (5)  adapted  this  approach 
for  fibrous  composites.  Each  slice  became  a  minibundle  of  fibers.  The  strength 
of  each  mlnlbundle  was  computed  by  loose  bundle  theory,  and  the  failure  of  the 
chain  of  minibundles,  each  having  the  length  5,  was  then  estimated  using  weakest 
link  theory. 
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Using  Che  geometrical  model  of  Rosen,  Zweben  (6)  investigated  the 
Influence  of  load  concentrations  caused  by  fiber  breaks.  Zweben  and  Rosen  (7) 
then  discussed  crack  growth  in  3-D  unidirectional  composites  on  the  basis  of 
the  chain  of  minibundles  model,  taking  into  account  stress  concentrations  in 
nearest  fibers,  and  looking  for  a  Griffith-type  instability.  Perhaps  as  a 
result  of  the  attempt  to  estimate  probabilities  for  all  possible  events,  the 
analysis  became  too  complicated  to  carry  through  for  3-D  composites.  They  did, 
however,  analyze  the  2-D  case  (single  ply  tape)  in  some  detail,  working  out 
the  expected  number  of  multiple  breaks  of  various  orders  as  a  function  of 
applied  stress.  Nevertheless,  they  did  not  succeed  in  establishing  a  theoreti¬ 
cally  based  failure  criterion,  stating  that  "one  can  only  say  that  the  expected 
failure  stress  will  be  bounded  by  the  expected  stress  at  the  first  break,  o^, 
and  the  expected  stress  at  which  N  adjacent  fibers  will  be  broken,  o  ”.  As  a 
practical  matter,  they  proposed  using  the  occurrence  of  the  first  multiple 
break  as  a  failure  criterion. 

An  improvement  in  regard  to  approaching  physical  reality  is  achieved  by 
taking  into  account  the  fact  that  in  many  composites  the  stress  gradient  in 
the  broken  fiber  get  so  high  that  the  matrix  is  unable  to  supply  the  shear 
stress  needed  to  establish  equilibrium.  This  appears  to  have  been  done  first 
by  Scop  and  Argon  (8).  As  a  result  either  the  matrix  yields  or  debonding 
between  fiber  and  matrix  sets  in.  In  either  event,  6  Increases  beyond  the 
elastic  transfer  length  and  becomes  a  function  of  applied  stress.  A  minibundle 
of  fibers  of  length  £  becomes  weaker  as  £  increases  due  to  the  Ueibull  sise 
effect,  and  failure  occurs  when  the  estimated  bundle  failure  stress  for  length 
£  is  equal  to  the  stress  required  to  produce  that  value  of  £.  It  has  been 
found  that  an  empirically-based  several-fold  Increase  in  £  must  be  assumed  in 
order  to  bring  theory  and  experiment  into  acceptable  agreement  in  the  case  of 
carbon-carbon  (9), 
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Harlow  and  Phoenix  are  placing  the  problem  of  composite  damage  and 
failure  on  a  more  rigorous  mathematical  foundation  than  has  theretofore  been 
available.  Their  approach  is  based  on  the  chain-of -minibundles  model  and 
seeks  to  establish  a  converging  sequence  of  upper  bounds  to  composite  strength. 

Some  appreciation  of  the  analytical  complexity  that  is  involved  in  a  rigorous 
treatment  can  be  obtained  by  noting  that  an  investigation  resulting  In  five 
papers  extending  over  a  three  year  period  (10-14)  has  been  required  to  solve 
the  problem  of  a  purely  elastic  single  ply  tape,  using  a  highly  simplified 
local  Load  sharing  rule.  It  is  the  opinion  of  the  present  writer  that  while 
this  approach  is  much  too  cumbersome  for  application  to  most  problesis  of 
practical  interest,  it  is  valuable  as  a  standard  to  help  assess  the  accuracy 
of  approximate  techniques  capable  of  a  wider  range  of  application.  Such  a 
use  will  be  made  in  the  next  section. 

In  spite  of  the  fact  that  almost  all  micromechanical  theories  of  composite 
tensile  strength  in  the  last  decade  and  a  half  have  been  based  on  the  chain-of- 
minibundles  approach  there  are  some  serious  deficiencies  which  are  intrinsic  to 
the  chain-of -bundles  model  Itself.  One  is  that  the  use  of  such  a  model  requires 
the  analyst  to  use  the  same  "ineffective"  length  for  a  multiple  fiber  break, 
regardless  of  order,  as  for  a  single  isolated  break.  This  assumption  is  par¬ 
ticularly  bad  when  debonding  occurs,  since  the  higher  the  multiplicity  of  the 
break,  the  higher  the  stress  concentration  factor,  and  therefore  the  lower  the 
debonding  threshold.  Another  deficiency  is  that  Interaction  between  cracks 
(multiple  fiber  breaks)  is  assumed  to  occur  only  if  the  cracks  are  in  the  same 
link  of  the  chain.  Thus,  the  chaln-of-bundles  model  can  treat  only  planar 
fracture  and  not  the  jagged  breaks  that  are  usually  observed  (see  for  instance. 

Fig.  1).  We  note  in  passing  that  the  occurrence  of  a  jagged  break  must  mean 

128 

>! 


i  \ 


I 


Chat  such  a  break  occurs  at  a  lower  stress  than  a  planar  break,  so  that  an 
accurate  theory  for  a  planar  break  must  be  unconservative.  Other  deficiencies 
Include  the  Inability  to  date  to  develop  a  satisfactory  failure  criterion 
except  in  cases  of  academic  interest. 


MULTIPLET  THEORY  OF  COMPOSITE  FAILURE 


A  new  theory  for  damage  accumulation  and  failure  is  proposed  herein  which 
rejects  the  chain-of-bundles  model  as  an  unnecessary  Intellectual  impediment, 
and  concentrates  instead  on  the  creation  and  growth  of  multiple  fiber  fractures. 
For  convenience,  it  will  be  referred  to  here  as  the  multiplet  theory  of  failure. 
For  the  sake  of  both  brevity  and  clarity,  only  a  simplified  version  of  the 
theory  will  be  discussed  in  this  section.  Complications  required  for  more 
accurate  modelling  of  composite  damage  will  be  mentioned  later. 

Consider  a  composite  containing  N  fibers  each  having  a  length  L  held 
together  by  a  matrix.  Damage  resulting  from  loading  is  assumed  to  consist  of 
breaks  in  the  fibers.  There  will  be  single  isolated  breaks  that  we  shall 
designate  as  singlets,  pairs  of  breaks  we  will  call  doublets,  as  well  as 
triplets,  quadruplets,  or  in  general  "i-plets". 

Every  singlet  Is  surrounded  by  n^  nearest  neighbors,  each  of  which  is 
subject  to  an  over stress  described  by  a  stress  concentration  factor  b^ 

(relative  to  the  nominal  fiber  stress)  at  the  cross  section  in  which  the  break 
is  located.  The  overstress  extends  over  a  distance  of  6^,  but  diminishes  with 
distance  from  the  cross  section  containing  the  break  (a  point  neglected  in  most 
previous  treatments) .  A  doublet  is  considered  to  have  n2  nearest  neighbors 
each  subjected  locally  to  a  stress  concentration  factor  b^,  with  the  varying 
overstress  applied  over  a  length  The  subscript  i  will  designate  the 

corresponding  quantities  for  an  i-plet. 

We  assume  in  the  present  section  that  fiber  failure  conforms  to  Weibull's 
2-parameter  representation,  i.e.,  the  cumulative  probability  of  failure  Pf  at 
stress  o  is  given  by 

Pj(o)  *  1  -  exp  (-Zk cm)  (1) 

when  the  stress  o  is  applied  uniformly  over  fiber  length  l;  k  and  a  are  the 
Weibull  parameters.  If  P^  <<  1,  this  expression  can  be  approximated  by 
Pf  -  kZam 
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If  there  are  N  fibers  each  of  length  L,  the  number  of  singlets  that  will 
have  been  created  by  the  time  the  stress  rises  to  o  is  given  by 

Ql  -  HP*  S  NLkcr*  (3) 

Ue  note  that  NL  is  the  total  length  of  fiber,  so  from  (3)  it  follows  that 
the  total  number  of  breaks  depends  only  on  the  total  length  of  fiber,  and 
not  on  the  dimensions  of  the  composite  per  se. 

The  first  failure  occurs  when  Qj  *  1  or  P{  ■  1/N.  From  (3)  the 
corresponding  stress  is 


The  second  break  occurs  when 


and  the  i ' th  when 


Here  o  is  nominal  stress,  i.e.,  the  uniform  stress  that  would  be  present  in 
all  the  fibers  in  the  absence  of  breaks. 

If  we  plot  the  logarithm  of  the  number  of  singlets  created  by  the  time 
the  stress  reaches  a,  0^,  against  log  a,  it  follows  from  (3)  that  we  obtain 
a  straight  line  of  slope  m,  as  shown  schematically  in  Fig.  2.  Ue  now  inquire 
into  the  value  for  Q.,,  the  number  of  doublets  that  have  been  created.  Note 
that  the  number  that  exist  at  stress  c  is  in  general  less  than  the  number  that 
have  been  created .  The  number  of  singlets  existing,  q^,  is  the  number  created 
less  the  number  that  were  turned  into  doublets,  i.e., 

qx  - -oL  -  q2  (6) 
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It  is  assumed  here  that  every  doublet  came  from  a  singlet,  every  triplet 
came  from  a  doublet,  etc.  Thus 


qi  =  Qi  “  Qi+1  (7) 

A  singlet  becomes  a  doublet  when  one  of  the  neighboring  fibers  breaks  in 
its  overstressed  region.  Such  a  fiber  has  a  probability  of  failure  given  by 

T  f4l/2  a  1  fV2  a 

F#  •  1  -  exp  -1c  I  o  dx  ;  2k  a  dx  (8) 

*  1  -V2  J  J° 

If  we  assume  the  stress  a  varies  linearly  from  k^o  to  a  over  the  length  6^/2, 


we  obtain 


A«r  T 

pf  ■ 21  Jc  [V  -  Tl(br1>5J 


bfl-i 

where  A  »  5  — -  (11) 

bi(br1)(m+1) 

is  the  "effective  length"  of  the  over-stressed  section  of  fiber  for  the 
purpose  of  calculating  P^.  The  ratio  A^/ 6^  is  shown  in  Fig.  3  for  various 
values  of  m  and  b^.  We  note  that 

\  \  when  bi  *  1  (12) 


■*  when  »  1  (13) 

These  limiting  relationships  are  familiar  from  Weibull  theory.  Equations 
(10)-(13)  apply  equally  well  when  i  is  used  as  the  subscript  instead  of  1. 


132 


Except  for  (8)  and  (9),  all  of  the  references  cited  assumed  <5  is  a  constant, 
independent  of  i,  and  is  equal  to  a  fixed  X  for  all  and  all  m. 

The  probability  that  a  given  singlet  becomes  a  doublet  is  (from  (10)) 


“  ^(b^a)  n^l^ 


(14) 


The  number  of  doublets  created  in  loading  to  stress  <J  thus  becomes 


Qj  “  Qj^hjO) 


(15) 


Generalizing  this  result. 


Q1+l  “  QjkOsjO)  n^ 

We  can  also  express  Q  ^  in  absolute  terms.  From  (15) 


(16) 


RLk0 


bi  -1*1 

ML 


(HLk) 


Similarly 


3  _  3m 

Q3  -  (NLk)  C 


and 


i-i 


j-l  (NL)l_i 


(17) 


(18) 


(19) 


From  (19)  we  see  that  a  plot  of  fnQj^  vs.  ina  is  a  straight  line  of 
slope  ira.  From  (16),  Qi+^  “  when 


k(b^cr)  niXi  -  1 


(20) 
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If  we  were  to  assume  that  b,  n  and  A  are  all  constant,  independent  of  i,  then 
all  Q's  would  be  equal  when  a  takes  the  value  implied  by  (20).  Actually 
n^  and  all  increase  as  i  increases.  Thus  the  value  of  a  at  which  Qi+1  =  Q^» 
namely 


keeps  decreasing  as  i  increases.  This  situation  is  illustrated  schematically 
in  Fig.  2. 

Now  since  each  multiplet  comes  into  being  by  fracture  of  a  fiber  adjacent 
to  a  multiplet  of  the  next  lower  order,  Q  ^  £  Q^.  Consequently  the  portion 
of  the  line  extending  beyond  its  intersection  with  has  no  physical 

meaning,  and  is  therefore  omitted.  The  envelope  of  the  lines  is  drawn  in 
heavily,  because  it  has  a  special  significance.  Over  the  stress  range  within 
which  an  i-plet  lies  on  the  envelope,  it  is  unstable.  That  is,  as  soon  as  it 
is  created  it  will  immediately  become  an  (i+l)-plet,  which  will  immediately 
become  an  (i+2)-plet,  etc.,  resulting  in  fracture  of  the  composite.  The 
failure  stress  is  the  lowest  stress  at  which  any  unstable  i-plet  is  present. 
This  is  the  stress  at  which  the  envelope  curve  (or  more  precisely,  broken  line), 
intersects  the  horizontal  line  Q  =  1,  or  £nQ  =0.  We  note  in  passing  that  only 
those  portions  of  a  Q-line  to  the  left  of  this  intersection  have  any  physical 
significance;  clearly  higher  stresses  than  the  fracture  stress  cannot  be 
applied  to  the  composite. 

The  envelope  plays  a  role  somewhat  analogous  to  that  played  by  the 
failure  line  in  a  conventional  Weibull  plot.  In  such  a  plot  y  ■  £n  in  P  ^  Is 
plotted  against  x  *  ina .  The  abscissas  are  thus  identical.  The  ordinates 
are  closely  related  since  for  sufficiently  small  Pf, 
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£n  InP”1  "  f.n[-£n(l-Pf)  ]  -  £nPf  (22) 

If  the  survival  probability  of  a  volume  of  a  material  is  given  by  the 
Weibull  expression 

P#  -  expl-k^a®]  (23) 

Then  the  Weibull  plot  becomes  a  straight  line  of  slope  m  as  shown  in  Fig.  4. 
The  line  for  volume  is  another  line  of  slope  m.  The  vertical  separation 
of  the  two  lines  is  tn(V^/V^) .  The  mean  fracture  stress  is  given  approxi¬ 
mately  by  the  intersection  of  the  fracture  line  with  the  line  corresponding 
to  Pg  =  0.5,  i.e.,  -Cn  £nP2^  =  -0.37.  The  change  in  fracture  stress  due  to 
the  size  effect  can  be  found  by  observing  how  the  intersection  point  changes 
as  the  fracture  line  is  shifted  vertically  a  distance  appropriate  to  the 
volume  change. 

In  the  case  of  the  composite,  it  follows  from  (19)  that  «  NL.  Since 
for  fixed  fiber  ratio,  NL  <*  V,  the  line  for  each  multiplet  in  Fig.  2  (and 
therefore  also  the  envelope)  moves  up  or  down  a  distance  Zn(V^/V^)  when  the 
volume  is  changed  from  to  V  The  effect  on  the  fracture  stress  is  found 
by  seeing  how  the  intersection  of  the  envelope  with  the  line  In  Q  =  0  changes 
with  such  a  vertical  shift. 

Another  way  of  showing  the  size  effect  is  to  plot  £na^  vs.  £nV^  where 
af  is  the  mean  fracture  stress.  The  result  for  a  Weibull  material  (i.e.,  one 
obeying  (23))  is  a  straight  line  of  uniform  slope  (-1/m),  as  shown  in  Fig.  5a. 
In  the  case  of  the.  composite,  according  to  the  present  theory  such  a  plot 
leads  to  a  broken  line.  Each  line  segment  corresponds  to  a  line  segment  of 
the  envelope  shown  in  Fig.  2  and  has  slope  (-1/im).  The  result  is  shown 
schematically  in  Fig.  5b.  This  finding  differs  sharply  from  those  for  most 
previously  published  theories.  For  instance.  Fig.  6a  taken  from  Zweben  (7) 
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shows  that  such  a  plot  leads  to  straight  lines  for  all  the  theories  he 
considered.  Fig.  6b,  taken  from  McKee  and  Sines,  shows  a  bilinear  relation 
for  such  a  plot.  Since  their  theory  predicted  a  6"  cable  (5  x  10^  volume 
elements)  would  fail  by  instability  of  a  septuplet,  the  present  theory  would 

O 

predict  at  least  6  line  segments  to  the  right  of  the  abscissa  10  in  Fig.  5b. 
The  flattening  out  of  the  curve  for  large  composites  predicted  by  the 
present  theory  may  be  very  important;  if  confirmed  by  additional  research, 
it  means  that  very  large  composite  structures  can  be  substantially  stronger 
than  has  heretofore  been  realized. 

In  many  cases  we  seek  not  the  mean  fracture  stress  of  the  composite  but 
the  stress  corresponding  to  some  prescribed  low  probability  of  failure.  This 
is  the  probability  that  an  unstable  crack  is  present  at  the  stress  under  con¬ 
sideration.  To  find  the  stress  at  which  =  0.01,  for  example,  all  that  is 
needed  is  to  read  off  the  stress  at  which  the  ordinate  of  the  envelope  is 
Q  =  0.01,  or  £nQ  =  fn  0.01  =  -4.61. 

The  preceding  discussion  points  out  a  number  of  similarities  between 
the  failure  envelope  for  composites  and  the  failureline  of  a  simple  Weibull 
material.  There  are  also  important  differences.  It  is  not  uncommon  to  hear 
numerical  values  cited  for  the  Weibull  modulus  m  of  a  composite.  This 
implies  that  the  failure  envelope  is  a  straight  line,  which,  as  we  have 
shown,  is  not  so.  One  can  identify  an  effective  modulus  for  calculating 
the  change  in  strength  in  going  from  one  size  to  another.  But  strictly 
speaking,  the  two  sizes  have  to  be  specified.  One  cannot  even  speak  of  a 
modulus  over  a  range  of  sizes,  since  as  we  shall  find  in  the  next  section, 
a  more  accurate  analysis  gives  a  smooth  curve  instead  of  a  broken  line  for 
the  failure  envelope.  The  effective  modulus  is  thus  the  slope  of  the  chord 
extending  between  the  two  points  on  the  failure  envelope  that  correspond  to 


the  two  sizes  under  consideration.  In  view  of  this  any  discussion  of  the 
Weibull  modulus  of  a  composite  should  be  carefully  qualified. 


MATHEMATICAL  CONSIDERATIONS 

It  will  be  obvious  to  critical  readers  that  the  engaging  simplicity  of 
the  theory  outlined  in  the  previous  section  is  largely  a  result  of  the 
numerous  mathematical  approximations  employed.  Most  obvious  perhaps,  is 
the  use  of  Eq.  (2).  Others  include  use  of  instead  of  -  Qi+^  for  the 
number  of  i-plets  actually  present,  neglect  of  the  fact  that  the  estimated 
probability  of  fracture  of  an  over stressed  fiber  element  should  really 
take  into  account  the  fact  that  it  survived  some  lower  stress,  neglect  of 
edge  effects,  etc.  An  accurate  estimate  of  the  errors  involved  in  such 
approximations  in  cases  of  practical  interest  would  involve  the  prohibitively 
difficult  task  of  finding  an  exact  solution.  A  rough  idea  of  the  errors 
involved  can  be  obtained,  however,  by  using  the  theory  outlined  to  solve  a 
2-D  problem  treated  exactly  (or  nearly  exactly)  by  Harlow  and  Phoenix. 

To  do  this  we  assume  with  them  that  the  probability  of  failure  of  a 
fiber  of  length  6  is  given  by 

Pf(6)  =  1  -  exp  (-(a/ °(j)m)  (24) 

or  using  our  approximation 

Pf(6)  =  (f-)m  (25) 

We  adopt  their  local  load-sharing  rule  for  a  single  ply  tape,  which  leads 
to 

b1  -  1.5  b3  -  2.5 

b3  =*  2.0  b^  =  3.0 

With  them  we  assume  that  n^  •  2,  and  that  X  -  6  *  constant.  Using  Eq.  (25) 
for  the  probability  P^(6)  that  any  arbitrarily  chosen  fiber  segment  of  length 
6  is  fractured,  together  with  Eq.  (16)  rewritten  in  the  form 
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(26) 


p<+1(6)  -  P  (fi)  2  b  “  (-f-)®1 
1+1  1  1  0  £ 

for  the  number  of  (i+l)-plets  divided  by  the  total  number  of  elements  NL/6, 
we  obtain  in  no  more  than  two  or  three  minutes  the  table  below. 

ln(o/a&)  •  -2  £n(a/a5)  -  -1 


In  P1 

-10 

-5 

In  P2 

-17.28 

-7.28 

In  P3 

-23.12 

-8.12 

In  P4 

-27.85 

-7.85 

In  ?5 

-31.66 

-6.66 

The  data  in  this  Table  are  displayed  graphically  in  Fig.  7.  These 
results  are  to  be  compared  with  those  given  by  Harlow  and  Phoenix  for  the 
same  problem.  Fig.  8.  In  making  the  comparison  due  account  must  be  taken 
of  differences  in  notation:  a  =  x;  m  ■  p;  P  «  W. 

The  most  conspicuous  difference  between  the  two  figures  is  that  Harlow 
and  Phoenix's  failure  envelope  is  a  smooth  curve  instead  of  a  broken  line. 
Obviously  it  should  be  a  smooth  curve.  This  follows  immediately  from  the 
fact  that  when  Eq.  (1)  is  used  for  P^,  0^  approaches  ^  asymptotically 
instead  of  crossing  it.  But  perhaps  the  most  important  point  to  note  is 
that  for  a  given  ordinate  the  abscissas  differ  by  less  than  0.05,  which  means 
the  failure  stresses  obtained  by  the  two  methods  will  generally  differ  by 
only  a  few  percent.  This  is  considerably  less  than  the  usual  experimental 
scatter  and  also  considerably  less  than  the  uncertainties  that  must  be  taken 
into  account  in  most  practical  problems.  These  uncertainties  include  the 
stress  at  which  debonding  commences,  the  length  that  is  debonded  for  each 
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multiplet  as  a  function  of  stress  level,  the  friction  between  fiber  and 
matrix,  etc.  One  of  the  primary  virtues  of  the  present  approach  is  that 
it  provides  an  analytical  framework  for  taking  all  such  effects  into 
account  as  soon  as  they  are  known,  without  encountering  the  obstacle  of 
mathematical  intractibility.  For  instance,  if  debonding  occurs,  b^,  6^, 
and  all  become  functions  of  the  applied  stress  a,  and  the  £n  vs  tn  a 
plots  become  curves  instead  of  straight  lines.  If  the  functional  relations 
b^o)  and  X^(o)  are  known,  however,  the  curves  can  be  constructed  in 
only  a  few  minutes  each,  and  their  envelope  becomes  the  failure  envelope 
for  the  problem  under  consideration. 
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RELATION  TO  GRIFFITH  AND  WEIBULL  CONCEPTS 

In  this  theory  there  are  two  levels  of  defect.  In  the  first  level 
are  pre-existing  fiber  defects  as  in  Weibull  theory  and  rupturing  of 
individual  fibers  is  governed  by  the  weakest  link  concept.  The  i-plets 
discussed  in  the  previous  sections  constitute  the  second  level  of  defect. 
Although  some  singlets  may  pre-exist  as  a  result  of  manufacturing  (and 
their  presence  could  be  readily  incorporated  into  the  theory)  the  prime 
emphasis  here  is  on  the  multiplet  structure  resulting  from  a  series  of 
fiber  ruptures  as  the  load  increases.  Thus  the  ultimate  failure  is  a 
second-level-defect  phenomenon,  in  which  the  "cracks"  or  multiplets  are 
not  pre-existing  but  are  caused  by  the  application  of  stress. 

It  is  the  multiplet  model  which  relates  to  Griffith's  crack  instability 

\ 

concept.  The  analogy  insists  in  fact  that  an  i-plet  becomes  unstable  at 
the  stress  at  which  the  Q1  curve  intersect  the  Qi+1  curve.  There  are  also 
significant  differences,  however.  For  instance,  Griffith  like  Weibull, 
thought  of  the  cracks  as  pre-existing,  whereas  here  they  are  a  result  of 
load  application.  More  important,  in  his  theory  all  cracks  of  a  given 
size  have  the  same  strength.  In  the  present  theory  there  is  a  range  of 
strengths  for  a  multiplet  of  any  given  rank.  If  only  one  is  present,  the 
strength  is  the  expected  strength  of  the  weakest  of  the  n^  fiber  segments 
of  length  X^  surrounding  the  multiplet.  If  M  are  present,  the  strength  of 
the  weakest  of  the  multiplets  is  determined  by  the  strength  of  the  weakest 
of  the  Mn^  surrounding  fiber  segments. 
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CONCLUDING  DISCUSSION 


The  theory  proposed  herein  is  believed  to  combine  the  best  features  of 
the  various  preceding  treatments,  and  in  so  doing  has  many  similarities  to 
them.  There  are,  however,  some  unique  features.  These  include:  a  better 
definition  of  "crack  instability"  as  the  concept  pertains  to  composites; 
a  quantitative  representation  of  the  defect  structure  of  a  composite  at  any 
given  stress  level,  which  can  be  used  to  establish  the  stress-strain  relation 
and  a  recognition  of  the  difference  between  the  elastic  transfer  lengths  and 
the  "effective  length"  X  in  the  context  of  weakest  link  theory.  Perhaps 
even  more  significant,  it  can  incorporate  important  effects  that  cannot  be 
treated  in  the  conventional  chin-of-minibundles  approach,  such  as  the 
variation  of  stress  transfer  length  with  multiplet  order,  the  interaction 
between  multiplets  in  different  minibudnles,  etc.  As  a  result,  the  present 
approach  is  believed  to  model  physical  reality  more  closely  than  its  pre¬ 
decessors.  In  addition,  it  provides  a  methodology  which  allows  the  various 
complications  discussed  to  be  incorporated  and  their  effects  analyzed  in  an 
orderly  fashion. 

To  apply  the  new  approach  to  practical  situations,  however,  several 
advances  in  our  understanding  are  required.  We  need  to  know  the  local 
elastic  stress  distribution  in  the  fibers  in  the  neightborhood  of  each 
miltiplet  order  that  is  likely  to  occur.  The  stress  at  which  matrix  yield¬ 
ing  or  debonding  begins  needs  to  be  known,  as  well  as  the  changing  influence 
of  these  factors  on  stress  distribution  as  the  stress  continues  to  increase. 
These  matters  require  an  understanding  of  the  inelastic  behavior  of  the 
matrix,  the  frictional  shear  stress  between  fiber  and  matrix  in  the  debonded 
region,  the  inception  and  growth  of  cracks  in  the  matrix,  etc.  The  investi¬ 
gation  of  such  questions  appears  to  be  the  next  logical  step  in  achieving  a 
useful  and  accurate  model  for  cumulative  damage  and  ultimate  failure  of 
fibrous  composites  in  tension. 
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Fig.  6  -  Zweben  (a)  and  McKee  and  Sines 
(b)  Predictions  of  Composite  Size-Strength  Relation 
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Probability  of  Failure 


Dimensionless  Load,  x/x8 

Ckl 

Fig.  8  -  Convergence  of  the  characteristic  distribution  function  Ww  "'(x) 
to  its  limit  W(x)  as  k  increases.  Fiber  elements  have  a 
Weibull  strength  distribution.  (Taken  From  Harlow  and  Phoenix) 
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The  Correlation  of  Viscoelastic  Behavior  with 
Fracture  Properties 
Co-Principal  Investigator 
Donald  J.  Plazek 


The  viscoelastic  behavior  of  elastomeric  and  epoxy  network  polymers  with 
known  molecular  structural  variations  is  being  determined  for  correlation  with 
the  fracture  and  static  mechanical  response.  The  preparation  along  with  static 
and  fracture  testing  is  being  carried  out  by  the  Akron  University  team  of  in¬ 
vestigators. 
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Following  the  establishment  of  the  crosslinking  procedures  for  the 
polybutadiene,  PB,  and  styrene-butadiene,  SBR,  elastomers  at  the  Polymer 
Science  Institute  of  Akron  University,  initial  samples  designated  PB-1  and  2 
and  SBR-1  were  provided  to  us  at  the  University  of  Pittsburgh  for  the  char¬ 
acterization  of  their  viscoelastic  behavior.  The  study  of  the  physical 
properties  of  these  elastomers  was  started  early  this  summer  (1980)  both 
in  Akron  and  in  Pittsburgh.  To  date  the  studies  of  the  viscoelastic 
behavior  have  been  solely  that  of  torsional  creep.  The  shear  creep  compli¬ 
ance  curves  at  temperatures  between  about  -50°  and  40°C  have  been  determined 
on  the  three  materials  PB-1  and  2  and  SBR-1.  PB-2  has  a  crosslink  concen¬ 
tration  that  is  about  1/2  of  that  of  PB-1.  The  indicated  equilibrium 
compliance  of  the  polybutadiene  with  the  larger  concentration  of  elastic 
chains,  the  PB-1,  is  indeed  lower  by  a  factor  of  0.72.  If  the  concentration 
of  chains  between  crosslinks  were  truly  twice  as  great  and  if  there  were  no 
or  relatively  few  effective  entanglements  a  factor  of  1/2  in  the  compliance 
would  be  expected.  A  comparison  of  the  reduced  creep  compliance  curves 
which  have  been  determined  are  shown  in  Figure  1  where  it  can  be  seen  that 
the  difference  in  the  curves  increases  with  the  passage  of  time.  The 
entanglement  compliance  plateaus  are  clearly  visible  because  of  the  charac¬ 
teristic  low  molecular  weight  per  entangled  chain  unit,  M  ,  for  polybutadiene, 
1  1 
^2000.  The  for  natural  rubber,  1,4  cis*  polyisoprene  is  5800  and  at 

comparable  crosslinking  levels  no  separate  entanglement  plateau  can  usually 

2 

be  distinguished.  Using  2.8  x  10  cin  /dyne  for  the  plateau  compliance 
seen  in  Figure  1  a  molecular  weight  between  junction  points  of  about  640  is 
indicated.  This  low  number  clearly  indicates  that  both  crosslinks  and 
entanglements  are  involved  in  determining  the  number  of  active  network 
chains.  The  data  displayed  in  the  figures  shown  here  have  been  reduced  to 
represent  the  response  over  an  enhanced  time  scale  range  at  a  reference 
temperature  of  0°C.  Time-temperature  reduction  procedures  were  utilized  in 
the  data  analysis. 

The  PB-2  data  are  shown  at  the  experimentally  expected  response  times 
for  0°C,  but  it  was  necessary  to  shift  the  PB-1  data  by  1.05  decades  of 
time  to  shorten  time^  of  response  to  achieve  the  superposition  seen  at  the 
short  time  end  of  the  reduced  time  scale.  The  slower  response  of  the  more 
highly  crosslinked  PB-1  network  suggests  that  the  glass  temperature,  T^,  for 
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the  PB-1  elastomer  is  higher  than  that  of  the  PB-2.  A  difference  in  T  in 

8 

this  direction  is  expected  from  the  difference  in  cross linking  level,  but 
this  has  not  been  confirmed  at  present  due  to  unexpected  behavior  observed 
in  linear  dilatometry  measurements  of  the  polybutadiene  sample. 

Since  most  of  the  specimens  are  being  cemented  into  the  testing  instru¬ 
ments  with  a  cyanoacrylate  adhesive,  absolute  compliance  levels  are  being 
checked  by  making  measurements  on  cylinders  with  different  dimensions  to 
insure  that  the  end  conditions  of  the  specimens  are  not  causing  appreciable 
errors . 

Figure  2  shows  the  retardation  spectra  L(lnr)  which  have  been  obtained 
from  the  reduced  compliance  curves.  Two  maxima  are  observed  reflecting 
the  presence  of  two  groups  of  viscoelastic  mechanisms  which  contribute  to 
the  two  plateau  compliance  levels  seen  in  Figure  1.  The  curve  characteriz¬ 
ing  the  response  of  the  looser  network  PB-2  shows  larger  values  in  the 
terminal  response  region  reflecting  its  higher  compliance.  Utilizing  these 
spectra  the  dynamic  mechanical  properties  of  these  two  samples  have  been 
calculated  and  are  presented  in  Figure  3.  Although  it  can  be  seen  that  the 
loss  tangent.  Tan  6,  of  the  looser  network  is  higher  at  low  frequencies, 
it  is  clear  that  this  is  so  only  because  the  storage  modulus,  G',  is  lower 
since  the  loss  moduli  are  virtually  identical  once  the  shift  in  time  scale 
of  a  decade  is  accounted  for.  This  result  is  somewhat  surprising  and  con¬ 
stitutes  a  prediction  that  the  fracture  behavior  of  these  materials  should 
be  nearly  the  same. 

Creep  measurements  are  also  partially  completed  on  a  series  of  three 

Viton  elastomers  whose  crosslinking  densities  are  0.127,  0.251  and  0.471. 

This  series  of  elastomers  provides  a  contrasting  molecular  structure  since 

they  are  copolymers  of  vinylidene  fluoride  and  hexafluoropropylene .  The 

measurements  to  date  on  these  fluornated  elastomers  definitely  show  more 

pronounced  long  time  losses  as  the  level  of  crosslinking  decreases.  This 

2 

is  in  accord  with  earlier  findings  on  natural  rubber. 

Two  creep  apparatuses  are  currently  being  used  for  the  measurements 
discussed  above.  A  third  one  is  being  modified  so  that  some  of  the  measure¬ 
ments  on  the  polybutadiene  and  the  styrene-butadiene  elastomers  can  be 
carried  out  at  temperatures  near  their  T  's.  A  low  temperature  liquid 
thermostat  is  being  constructed  which  should  make  measurements  possible 
down  below  -80°C,  if  necessary. 
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Storage,  G',  and  loss,  G",  moduli  of  the  polybutadienes  plotted  logarthmically  as  functions 
of  the  reduced  angular  frequency  scale  along  with  linear  values  of  the  loss  tangent.  Tan 
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FEASIBILITY  DEMONSTRATION  OF  MELT  FRACTURE  SET-UP 


J.  T.  Lindt  and  Y.-S.  Wang 


Review  of  Objectives 

1.  To  design  an  apparatus  for  studies  of  melt  fracture 

2.  To  generate  flow  conditions  under  which  melt  fracture  occurs 

3.  To  identify  and  measure  the  characteristic  variables 

4.  Interpretation  of  data 

Summary  of  Progress 

The  above  objectives  have  been  set  for  a  period  of  one  year.  During 
the  7-month  period  covered  by  this  communication  an  experimental  set-up  has 
been  developed  and  flow  regimes  associated  with  melt  fracture  have  been 
generated.  The  first  part  of  the  experimental  program  has  centered  on  the 
nature  of  pressure  fluctuation  associated  with  melt  fracture  of  polyethylene 
In  general,  the  pressure  fluctuations  show  a  fairly  regular  periodicity. 

To  analyze  the  nature  of  the  pressure  fluctuations,  a  faster  recording 
system  is  required. 

Description  of  the  Experimental  Set-up 

To  generate  a  flow  regime  associated  with  melt  fracture,  an  apparatus 
has  been  designed  and  subsequently  built  as  shown  in  Fig.  1.  The  experi¬ 
mental  set-up  can  be  divided  into  three  parts  fulfilling  the  following 
functions: 

1.  Generation  of  homogeneous  melt, 

2.  Expelling  the  melt  at  a  steady  rate, 

3.  Generation  of  melt  fracture. 

The  melt  is  generated  in  a  plasticating  1"  NRM  extruder. 

The  hydraulic  system  to  discharge  the  melt  through  the  test  cell  con¬ 
sists  of  two  cylinders  mounted  back  to  back.  The  right  cylinder  serves  as 
a  reservoir  of  melt  while  the  left  one  contains  hydraulic  oil  and  is  to 
drive  the  expelling  unit.  Both  cylinders  have  a  3  1/4"  bore  and  12"  stroke 
(manufactured  by  Tomkins-Johnson  Company) .  The  pressure  required  is 
generated  by  a  hand  pump  (model  160-6,  Templeton  Kenley  and  Company)  con¬ 
nected  to  an  oil  accumulator  (R4-20-5Co,  Hydrodyne  Industries) . 
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The  test  cell  designed  to  generate  and  observe  melt  fracture  is  shown 
in  Fig.  2.  It  consists  of  two  Identical  profiled  aluminum  plates  that  bound 
the  flow.  The  distance  between  the  plates  hg  is  adjustable  (so  far  hg  ■  2mm 
has  been  used).  The  side  plates  fitting  the  vertical  slots  have  been  made 
of  stress-free  glass  for  optical  measurements  of  stress;  a  "heavy-duty"  set 
of  aluminum  plates  has  also  been  made  to  be  used  during  the  developmental 
stages  of  the  project,  or  in  measurements  dedicated  to  hydrodynamic  pressure 
determinations. 

The  instrumentation  includes  Dynisco  pressure  transducers  (model  PT 
420/411),  and  Omega  iron-cons tantan  thermocouples;  the  outputs  were  measured 
using  Houston  Instruments  two-channel  recorsers  (model  BS  239-51) .  Locations 
of  the  sensors  are  shown  in  Fig.  2. 

Experiments  with  LDPE 

During  the  developmental  stage  a  low  density  polyethylene  was  used. 

Although  most  of  our  effort  during  this  initial  period  was  devoted  to 

designing  and  developing  the  experimental  set-up,  a  number  of  tests  have 

5  2 

been  run.  Melt  fracture  has  been  observed  at  a  stress  level  10  N/m  in 
agreement  with  some  older  data. ^  A  set  of  data  illustrating  our  experi¬ 
mental  conditions  is  shown  in  Table  1. 


Table  1 

EXAMPLE  OF  FLOW  MEASUREMENTS  (ARCO  LDPE) 


FLOW 

RATE 

[g/min] 

MELT 

TEMPERATURE 

[°C] 

PRESSURE 

GRADIENT 

PSI/cm 

SHEAR 

STRESS 

N/m2 

REMARKS 

18.5 

141 

88 

6.1  x  104 

Shark  Skin 

23 

141 

106 

7.3  x  104 

30.5 

141 

111 

7.6  x  104 

_  »»_ 

39 

141 

120 

8.4  x  104 

Melt  Fracture 

46 

141 

114 

7.8  x  104 

49 

141 

111 

7.6  x  104 
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The  appearance  of  "shark  skin"  on  the  surface  of  extrudate,  and  the 
gross  distortions  of  it,  were  accompanied  by  pressure  fluctuations  in  the 
test  cell.  Although  they  seemed  to  show  a  fairly  regular  amplitude  and 
periodicity,  some  additional  signal  conditioning  and  faster  recording  capa¬ 
bilities  are  required  for  a  successful  analysis  of  the  nature  of  the 
pressure  fluctuations  which  are  believed  to  play  an  important  role  in  the 
propagation  of  flow  instabilities  under  consideration. 

Recommendations  for  Future  Research 

In  view  of  the  premature  termination  of  this  project  after  eight  months 
devoted  to  demanding  development  of  the  experimental  set-up,  the  following 
suggestions  are  made  for  a  possible  follow-up  in  the  future: 

1.  To  analyze  the  frequency  and  amplitude  of  the  pressure  fluctuations 
as  functions  of  the  stress/flow  rate  level,  of  the  shape  of  the 
inlet  regions  (c.f.  Fig.  2),  and  of  the  molecular  parameters  of  the 
given  polymer. 

2.  To  perform  a  theoretical  analysis  of  the  oscillatory  flow  involved. 


Reference 

Vlachopoulos  and  T.  W.  Chan,  Journal  of  Applied  Polvmer  Science, 
21,  1177  (1977). 
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FRACTURE  IN  FULLY  PLASTIC  BODIES 


Introduction 

Fracture  during  hot  working  processes  is  an  increasingly  important 
problem  because  of  the  use  of  highly  alloyed  high  temperature  materials 
in  modern  defense  and  industrial  systems.  Ductile  fracture  during  cold 
working  processes  is  well-known  to  consist  of  void  initiation  at  inclusions 
followed  by  growth  and  coalescence  along  planes  of  maximum  shear.  In  hot 
working  processes,  fracture  mechanisms  are  complicated  by  dynamic  structural 
changes  at  high  temperature  and  such  mechanisms  are  yet  to  be  elucidated. 

In  a  previous  report,  a  test  system  was  described  for  the  evaluation 
of  fracture  at  hot  working  at  controlled  temperatures  and  strain  rates. 

The  system  can  perform  upset  compression  and  bend  tests  at  temperatures  up 
to  1200°C  and  strain  rates  up  to  20  sec  ^ .  In  addition,  a  quenching  apparatus 
has  been  devised  which  permits  locking  in  of  the  hot  worked  metallographic 
structure  for  determination  of  mechanisms  of  high  temperature  fracture. 

Since  fracture  in  complex  alloy  systems  is  complicated  by  the  large 
number  of  phases  present,  it  is  useful  to  start  by  examining  a  simple,  single¬ 
phase  material  and  then  increasing  the  complicity  of  the  alloy  and  fracture 
mechanisms.  For  this  purpose,  a  35  Co,  35  Ni,  20  Cr,  10  Mo  alloy  was  chosen 
because  at  temperatures  above  900°C,  it  is  single  phase  fee.  The  alloy  is 
also  vacuum  cast  so  it  has  a  negligible  amount  of  inclusions  to  complicate 

the  fracture  process.  Finally,  the  alloy  also  has  a  very  low  stacking  fault 
2 

energy  (25  erg/cm  )  so  that  it  readily  undergoes  dynamic  recrystallization. 

The  cobalt-nickel  alloy,  known  as  MP35,  was  tested  in  bending  and  friction¬ 
less  compression  at  temperatures  from  1010°  to  1205°C,  and  at  strain  rates 
from  0.01  to  10  sec  Bend  tests  were  performed  to  various  increments  of 
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deformation  to  permit  tracking  of  the  fracture  mechanisms  and  homogeneous 
compression  tests  were  performed  to  monitor  the  dynamic  recrystallization 
and  flow  stress  behavior  of  the  material. 

Fracture  testing 

Qualitatively,  the  fracture  behaviour  of  MP35  alloy  in  plane  strain 
bending  can  be  represneted  in  a  matrix  of  temperature  and  strain  rate.  Fig.  1. 

The  total  bending  strain  is  0.6  in  each  test.  At  a  given  temperature,  the 
effect  of  strain  rate  in  arresting  fracture  can  be  seen  immediately  along 
each  row.  The  effect  of  temperature  at  constant  strain  rate  can  be  assessed 
from  the  vertical  columns.  Note  that  two  regions  are  nearly  free  of  surface 
cracking:  one  at  the  low  temperature,  high  strain  rate  corner  and  the  other 
at  the  high  strain  rate,  high  temperature  corner  on  Fig.  1. 

Flow  Stress  Measurements 

The  flow  stress  as  a  function  of  strain,  strain  rate  and  temperature 
was  measured  by  performing  constant  strain  rate,  isothermal,  homogeneous 
compression  tests  on  cylinders.  Combinations  of  strain  rates  and  temperatures 
similar  to  those  used  for  the  bend  tests  were  selected.  Flow  stress  curves 
for  the  strain  rates  of  0.05,  1  and  10  sec  1  are  presented  in  Fig.  2,  3,  and 
4. 

The  peak  followed  by  a  drop  in  flow  stress  indicates  that  the  MP35  exhibits 
extensive  dynamic  recrystallization  over  the  range  of  temperatues  and  strain 
rates  investigated.  The  concurrent  very  limited  capacity  for  dynamic  recovery 
is  also  apparent.  This  behaviour  is  consistent  with  its  low  value  of  the 
stacking  fault  energy.  That  is  to  say,  the  alloy  first  work  hardens  very 
rapidly  with  very  little  dynamic  recovery  taking  place.  When  the  critical 
strain  for  dynamic  recrystallization  is  achieved,  it  undergoes  the  characteristic 
softening  until  the  steady  state  is  reached  and  the  flow  stress  remains  constant. 
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Metallographic  Results 

Optical  metallography  of  sections  of  the  bend  specimen  was  adopted  as 
the  main  tool  to  characterize  the  structure.  In  general,  sections  parallel 
to  the  plane  of  flow  at  mid-breadth  of  the  specimen  were  examined,  both 
unetched  and  etched.  Figure  5  is  a  typical  section  of  a  well  developed  crack, 
unetched.  The  90°  change  in  direction  of  the  cracks  is  typical  of  well  developed 
cracks,  following  lines  approximately  at  45°  with  the  free  surface.  Figure  6 
shows  etched  sections  of  specimens  bent  at  1175°C,  and  at  strain  rates  of 
0.01,  0.1,  1  and  10  sec  ^ .  This  figure  should  be  examined  in  conjunction  with 
Fig.  7,  which  shows  the  pertinent  flow  stress  data  at  1175°C  as  a  function 
of  strain  rate.  In  this  plot,  the  maximum  strain  achieved  in  bending  (0.6) 
has  been  indicated  by  the  vertical  dashed  line. 

Two  things  are  apparent.  First,  the  amount  of  cracking  diminished 
strikingly  as  the  strain  rate  increases,  as  observed  previously  in  Fig.  1. 

Second,  the  structure  of  the  matrix  also  changes  with  strain  rate,  in 
a  manner  that  is  consistent  with  the  flow  stress  data  for  this  temperature. 

Both  specimens  deformed  at  0.01  and  0.1  show  the  equiaxied  fine  grains 
characteristic  of  the  dynamic  recrystallized  steady  state.  It  can  be  readily 
seen  from  the  flow  stress  data  that  the  maximum  bending  strain  lies  well  beyond 
the  onset  of  the  steady  state  for  these  strain  rates.  On  the  other  hand,  for 
the  strain  rate  of  10  sec  \  the  maximum  bending  strain  lies  before  the  onset 
of  the  steady  state,  and  indeed,  a  partially  recrystallized  structure  can  be 
observed.  For  a  strain  rate  of  1  sec  \  an  almost  fully  recrystallized  structure 
can  be  observed,  again  consistent  with  the  position  of  the  maximum  bending 
strain  with  respect  to  the  onset  of  the  steady  state. 
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Mechanics  of  Plane  Strain  Bending 

It  will  be  noted  that  plasticity  theory  establishes  that,  for  the  case 
of  plane  strain  bending,  the  maximum  shear  stress  develops  along  lines  that 
start  at  45°  to  the  free  surface  of  the  specimen,  Fig.  8.  It  can  also  be 
established  that  along  these  lines  the  shear  strain  rate  can  reach  very 
high  values.  In  other  words,  plasticity  theory  allows  the  possibility  of 
localization  of  the  deformation  along  certain  lines.  In  the  case  of  plane 
strain  bending,  these  lines  are  at  45°  to  the  free  surface.  Moreover,  along 
these  same  lines,  the  shear  stresses  are  maximu. 

During  bending,  the  convex  surface  undergoes  tensile  stress  and  strain, 
and  the  magnitude  of  these  quantities  diminishes  toward  the  mid-plane  of  the 
test  specimen.  On  the  concave  side,  the  stresses  become  compressive.  In 
addition,  a  notch  is  placed  on  the  concave  side  to  increase  the  tensile  strains 
achievable  on  the  convex  surface. 

Microscopic  Response  and  Fracture 

From  a  microscopic  point  of  view,  as  the  bending  of  the  plane  strain 
specimen  proceeds  at  constant  temperature  and  strain  rate,  the  critical  strain 
for  dynamic  recrystallization  will  be  achieved  in  the  zone  of  maximum  strain 
being  considered.  The  strain  is  maximum  at  the  point  in  the  free  surface 
opposite  to  the  punch,  and  diminishes  both  outward  along  this  surface  and 
inward  along  the  radial  direction.  As  the  critical  strain  is  reached  within 
this  zone  the  new,  small,  equiaxied  dynamically  recrystallized  grains  begin 
to  form  along  the  existing  boundaries  of  the  hot  worked  grains.  Grains  outside 
this  zone  remain  hot  worked  to  the  strain  reached  at  any  particular  point. 

A  polished  and  etched  section  of  a  specimen  fast  quenched  at  this  stage  is 
shown  in  Fig.  9.  Note  also  that,  by  chance,  some  of  the  old  grain  boundaries 
just  covered  by  the  new  grains  will  be  lying  at  approximately  45°  to  the  free 
surface,  or  along  the  lines  of  maximum  shear  stress. 
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Consider  also  a  hot  worked  grain  boundary  at  45°  to  the  free  surface. 

Just  before  the  new  grains  nucleate,  the  level  of  the  maximum  shear  stress 
along  this  boundary  is  given  by  the  peak  in  Fig.  7.  As  the  new  grains  nucleate 
and  cover  the  boundary,  they  can  only  support  a  level  of  stress  given  by 
the  steady  state  region  in  Fig.  7.  Such  an  imbalance  of  the  stresses  or, 
in  other  words,  this  localized  loss  of  load  carrying  capacity,  is  expected 
to  lead  to  the  localization  of  plastic  flow  along  the  slip  lines. 

Thus,  the  mechanism  of  fracture  in  this  particular  alloy  can  be  characterized 
in  Fig.  10.  Shear  stresses  during  bending  are  superimposed  on  regions  of 
localized  dynamic  recrystallization  along  prior  grain  boundaries.  The  softer, 
recrystallized  grains  permit  localization  of  shear,  with  subsequent  generation 


of  voids  and  cracks. 


Initiation  of  fracture 


In  order  to  determine  the  nature  of  fracture  initiation  leading  to  gross 
cracking  along  dynamically  recrystallized  grains,  a  series  of  specimens 
were  bent  to  strains  of  0.2,  0.35,  0.45,  0.55,  and  0.60  respectively 
at  a  temperature  of  1110°C  and  strain  rate  of  0.05.  The  specimen  bent  to 
a  strain  of  0.2  just  before  the  peak  strain  of  0.25,  and  fast  quenched  provided 
the  most  useful  information.  Figure  11  shows  that  indeed  very  little  dynamic 
recrystallization  has  taken  place  along  prior  grain  boundaries.  Also,  only 
a  few  tiny  starting  microcracks  (long  arrow)  can  be  found.  More  important, 
considerable  subsurface  damage  in  the  form  of  microvoids  along  grain  boundaries 
(short  arrows)  can  be  observed.  Higher  magnification  views  of  these  voids 
(Fig.  12)  showed  that  they  formed  only  at  these  grain  boundaries  just  covered 
by  the  soft  layer  of  dynamic  recrystallized  grains.  Very  few  if  any  of  these 
voids  were  found  along  neighboring  grain  boundaries  along  hot  worked  grains 
at  which  dynamic  recrystallization  had  not  taken  place. 

No  cavities  were  found  at  the  original  boundaries  between  hot  worked  grains. 
Instead  extensive  cavity  formation  was  found  at  those  grain  boundaries  at 
which  dynamic  recrystallization  took  place.  This  is  interpreted  on  the  basis 
of  the  argument  advanced  before:  The  high  stresses  developed  by  deforming  the 
hot  worked  grains  are  sufficient  to  promote  extensive  sliding  when  the  old 
boundaries  are  replaced  by  the  soft  layer  of  recrystallized  grains. 

A  final  comment  on  the  role  of  inclusions  follows.  Although  the  volume 
fraction  of  inclusions  is  very  low,  some  square  shaped  inclusions  were  found. 

Fig.  13,  both  along  grain  boundaries  and  within  the  grains.  In  both  cases  they 
lead  to  microvoid  formation  at  rather  low  strains  (less  than  0.2,  from  the  specimen 
bent  to  a  maximum  strain  of  0.2  at  1110°C/0.05  sec  S .  The  inclusions  appear 
to  be  very  hard,  since  their  shape  is  not  altered  by  deformation.  A  single  one 
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along  a  recrystallized  grain  boundary  was  found  cracked.  Fig,  14,  This  is 
taken  as  indicative  of  a  high  stress  concentration  induced  by  the  sliding, 
as  argued  before. 

Conclusion 

A  mechanism  for  the  first  stage  of  the  fracture  process,  void  nucleation 
has  been  described.  This  occurs  by  grain  boundary  sliding  induced  by  dynamic 
recrystallization.  Inclusions  play  a  very  minor  role  due  to  their  negligible 
volume  fraction.  The  voids  at  this  stage  are  truly  microscopic  and  their  presence 
does  not  give  a  complete  explanation  of  the  macroscopic  crack  observed  and 
their  features. 

By  chance,  at  certain  points  in  the  specimen  and  at  about  45°  to  the  surface, 
a  layer  of  soft  recrystallized  grains  containing  voids  will  develop.  Superimposing 
on  this  layer  the  maximum  stress  acting  along  the  slip  lines  (Fig.  10),  one 
obtains  a  configuration  that  leads  to  the  localization  of  the  plastic  flow  along 
these  directions. 
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Figure  3.  Flow  stress  of  MP35N  as  a  function  of  temperature  at  con¬ 
stant  strain  rate  of  1  sec.“l 
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Figure  7.  Flow  stress  of  MP35N  at  1175°C,  as  a  function  of 
strain  rate 
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Elementary  mechanics  of  plane  strain  bending, 
a:  specimen  (w/h  =  8)  and  principal  axes,  b:  definition 

of  stresses,  c:  state  of  stress  near  surface,  d:  Mohr's 
diagram  for  plastic  strains,  e:  illustration  of  stresses 
and  strains  at  surface  of  specimen,  valid  for  deviatoric 
stresses  and  total  or  deriatoric  strains 


Figure  11.  Plane  strain  bending  specimen  of  MP35N  bent  to  a  strain 
of  0.2  just  before  the  peak  strain  of  0.25  (1110°C/0.05 
sec  .“*•) 
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Introduction 

It  is  well  known  that  processing  substantially  alters  the  thermal 
behavior  of  ceramic  materials.  Reduction  in  particle  size  increases 
the  driving  force  leading  to  densification  while  increasing  the  specific 
surface  area  which  leads  to  increased  reactivity  with  the  environment  or 
other  phases  present.  A  well  known  example  is  the  addition  of  fine  quartz 
to  a  porcelain  body.  The  same  is  true  for  the  nonoxide  ceramics  such  as 
SiC.  which  have  a  tendency  to  react  v:  .  the  environment  result  in.-  in 
tenancious  oxide  layers  even  at  low  temperature.  This  reaction  is  extremely 
difficult  to  avoid  and  even  mare  difficult  to  reverse.  Consequently 
attempts  to  sinter  silicon  carbide  by  reducing  the  particle  size  to  the 
submicron  level  have  had  only  moderate  success  due  to  the  substantial  in¬ 
crease  in  oxide  contamination.  Partial  sintering  results  from  the  isola¬ 
tion  of  oxide  islands  within  the  paths  for  mass  transport.  The  obvious 
solution  is  to  introduce  a  new  phase  that  will  react  with  the  oxide  to 
yield  a  liquid  phase  at  high  temperature  and  increase  the  sinterability . 
Unfortunately,  the  reactions  that  occur  are  reversible  and  lead  to  a 
decrease  in  high  temperature  creep  resistance.  By  careful  choice  of  the 
phase(s)  added.  Carborundum  is  able  to  sinter  t- SiC  to  high  density  while 
maintaining  desirable  mechanical  properties. 

One  solution  to  the  problem  of  utilizing  very  fine  powders  that  are 
contaminated  with  oxide  is  to  reduce  the  oxide  thus  exposing  a  clean, 
highly  active  surface  of  SiC.  Our  approach,  which  to  the  best  of  our 
knowledge  has  not  been  attempted,  has  been  to  expose  the  contaminated 
powder  to  a  low  temperature  plasma  (glow  discharge)  followed  by  sintering 
at  higher  temperatures . 
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Research  Objectives 

The  objective  of  this  research  is  to  study  the  influence  of  the  pre¬ 
paration  of  simple  shape  products  of  nonoxide  silicon  based  ceramics  on 
the  reproducibility  of  the  final  product  using  fracture  mechanics  as  the 
major  indicator.  Efforts  are  being  directed  towards  consistency  of  the 
samples  and  not  towards  optimization  of  the  properties.  Powders  con¬ 
sidered  in  this  project  include  Carborundum's  alpha-Silicon  Carbide  Premix, 
Lonza  Chemical's  1200  mesh  black,  1200  mesh  green,  and  UF  Carbogran,  and 
Aremco's  alpha-Silicon  Carbide.  Transverse  rupture  specimens  will  be 
cut  from  larger  disks  prepared  in  a  traditional  dry  press/sinter  cycle 
or  pressed  directly  to  size.  Uniformity  is  insured  by  cutting  specimens 
from  the  larger  specimens. 

Processing  variables  to  be  considered  in  this  work  include: 

Powder  characteristics; 

particle  size  and  distribution, 
crystallographic  modifications  of  the  powders, 
particle  morphology, 
particle  surface  topography. 

Compact  ion  variables ; 

ultimate  applied  pressure, 
initial  (green)  density, 
binders/lubricants , 
delaminations. 

Sintering  variables; 

heating  rate, 

soak  time  at  temperature, 

surface  activation  using  low  temperature  plasma, 
plasma  variables  (initiation  temperature,  duration, 
type  of  gas ) . 

The  Carborundum  premix  powder  is  to  be  used  as  a  control  since  it  can 
bo  sintered  to  high  density  yielding  excellent  mechanical  properties.  It 
will  be  used  as  received  to  allow  a  better  comparison  of  the  other  powders 
before  and  after  processing  changes. 
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Status  of  Research 


Characteristics  of  the  Lonza  and  Carborundum  powders  were  quantitative- 
vely  determined  early  in  this  reporting  period.  Particle  size  and  dis¬ 
tribution  were  determined  by  direct  measurement  using  SEM  of  dispersed 
powders  and  is  given  in  Table  II  of  the  attached  paper.  At  least  95%  of 
each  powder  is  less  than  5y  m  as  compared  with  the  Carborundum  powder  hav¬ 
ing  >98%  less  than  1  um.  The  Carborundum  powder  is  spray  dried  and  as 
such  has  excellent  flow  properties  and  good  compaction  characteristics. 

Determination  of  the  crystallographic  polytype  makeup  of  silicon 
carbide  is  complicated  due  to  several  overlapping  x-ray  diffraction  lines 
for  the  major  phases.  Hoxcever,  computer  analyses  at  the  x-ray  intensities 
allows  a  semi-quantitative  determination  of  the  major  polytypes;  i.e, 
i  —  SiC  (15R  6H,  4H)  and  S-SiC  (3C).  The  results  for  the  Lonza  powders 
are  given  in  Table  I,  attached  paper;  the  Carborundum  powder  is  11%  15%R, 
86%  6H,  1.5%  4H,  and  1.5%  3C.  This  powder  is  advertised  as  pure  a  -  SiC 
but  does  contain  the  6  -  SiC  poly type.  The  Lonza  powders  have  a  similar 
crystallographic  makeup  which  allows  for  comparisons  between  all  powders 
as  a  function  of  processing. 

During  the  previous  reporting  perxod  work  centered  on  designing  and 
building  an  induction  vacuum  furnace  capable  of  2100C  (Fig.  1).  Simple 
disks  (13mm  x  5mm)  were  pressed  to  58%  density  and  utilized  to  isolate 
problems  with  the  heading  schedule,  compaction  problems,  and  SEM  analysis. 
Disks  were  prepared  by  first  adding  3-4%  PVA  to  the  powder  as  a  binder/ 
lubricant  and  pressed  uniaxially  at  100,000  psi. 

During  this  reporting  period  experimentation  with  different  levels 
of  binder  and  proper  aging  prior  to  pressing  has  resulted  in  disks  13mm 
and  31mm  in  diameter  being  compacted  to  60%  density  using  70-75,000  psi. 
Delaminations  during  pressing  have  been  eliminated.  Higher  initial  den¬ 
sities  are  not  desired  for  two  reasons;  delaminations  during  firing,  and 
any  enhancement  resulting  from  processing  changes  are  easier  to  detect  with 
medium  to  low  densities. 

The  most  significant  processing  variable  considered  to  date  is  the 
removal  of  particle  surface  contamination  by  exposure  to  a  low  temperature 
plasma  during  the  sintering.  Details  are  presented  in  the  attached  paper. 


Current  efforts  are  directed  towards  a  maximization  of  the  plasma  activation 
leading  to  higher  final  densities. 

It  has  been  demonstrated  that  the  plasma  variables  are  more  influential 
than  the  powder  characteristics,  including  prior  mechanical  treatment.  Each 
powder  has  been  pressed  to  60%  density  and  exposed  to  a  plasma  for  30  min. 
at  temperatures  from  800-1700C.  Morphological  changes,  as  observed  using 
the  SEM,  begin  at  1585-1600C  in  as  little  time  as  3  min.  Activation  on 
an  atomic  scale  presumably  occurs  at  a  lower  temperature.  Using  the  SEM 
only  the  onset  of  diffusion  resulting  from  such  activation  is  detected. 

Larger  disks  (31mm  x  5mm)  have  been  pressed  and  sintered  using  the 
Brew  High  Temperature  furnace.  This  allows  sintering  at  2200C  in  inert 
atmospheres  or  vacuum.  Considerable  changes  are  being  made  to  the  furnace 
to  allow  for  the  sintering  of  larger  disks  (51mm)  necessary  for  prepara¬ 
tion  of  transverse  rupture  specimens.  Disks  of  Carborundum  powder  have 
been  sintered  to  92%  density  at  2150C  in  argon.  These  will  serve  as  the 
control  for  comparisons  with  the  powders  activated  in  the  plasma  followed 
by  sintering  to  the  same  temperature. 

The  Lonza  powders  are  now  being  reduced  in  size  using  an  air-jet-mill 
grinder  to  a  level  comparable  with  the  Carborundum  powder, 
minimize  grain  size  effects  during  fracture  testing. 
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I .  INTRODUCTION 

Sintering  of  ceramic  materials  in  a  glow  discharge  has  been  demonstrated 

1-4 

to  be  a  viable  consolidation  process.  Rapid  sintering  rates  have  been  reported 
in  addition  to  the  development  of  favorable  microstructural  features.  More 
importantly,  a  glow  discharge  plasma  may  enhance  normal  sintering  by  activation 
of  powder  surfaces. 

The  driving  force  for  sintering  is  directly  dependent  on  the  surface 

energy  of  the  constituent  powder.  Factors  that  alter  surface  energy  generally 

affect  the  sinterability  of  the  powder.  In  any  powder  array  the  excess  surface 

energy  can  be  lowered  by  faceting,  change  in  shape,  or  by  the  formation  of  necks 

between  adjacent  particles.  For  an  isolated  particle,  or  one  isolated  by  a 

second  phase  such  as  an  oxide  layer,  the  surface  energy  can  be  further  changed 

by  the  presence  of  the  oxide  or  second  phase.  Non-oxide  ceramics  such  as  SiC 

and  Si.N .  are  difficult  to  sinter  due  to  a  high  degree  of  covalency  and  the 
3  4 

presence  of  tenacious  oxide  layers  that  are  thermodynamically  stable  in  an 
oxygen  containing  environment.  The  use  of  submicron  powders  increases  powder 
activity  but  also  increases  the  tendency  towards  low  temperature  surface 
oxidation.  Addition  of  sintering  aids  to  SiC  contaminated  with  oxides  layers 
has  been  successful  in  sintering  to  high  densities,  most  probably  due  to 
formation  of  intergranualar  liquid  phases.  It  is  desired  to  sinter  SiC  to 

Presented  at  the  82nd  Annual  Meeting,  The  American  Ceramic  Society,  Chicago,  IL., 
April  27,  1980;  Basic  Science  Division,  Paper  27-13-80. 
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high  densities  without  additives  since  residual  phases  may  affect  the  high 
temperature  mechanical  properties. 

2- 

Rapid  sintering  rates  of  oxide  ceramics  have  been  achieved  using  plasmas. 

In  some  oxides,  partial  melting  of  the  surface  may  lead  to  a  modified  liquid 
phase  sintering  condition.  However,  it  appears  the  major  cause  of  enhanced 
sintering  is  the  rapid  attainment  of  temperatures  sufficiently  high  such  that 
diffusive  mass  transport  mechanisms  leading  to  densif ication  overwhelms  those 
leading  to  surface  changes  only.  Alternatively,  a  plasma  could  be  utilized 
to  break  up  oxide  surface  layers  or  segregated  contaminants  creating  a  highly 
active  surface.  The  sintering  process  is  then  reduced  to  that  for  a  pure 
system  having  active  surface  diffusive  mechanisms.  Plasma  activation  can  be 
used  in  combination  with  conventional  heating  allowing  for  better  temperature 
control  than  sintering  in  a  plasma  alone. 

The  surface  condition  of  any  powder  is  of  major  importance  in  the  sintering 
process.  Total  removal  of  oxide  layers  from  SiC  must  be  accomplished  if  the 
initial  stage  of  sintering  is  to  lead  continuously  to  final  densification. 

U 

7  8 

It  has  been  demonstrated  by  Heath  and  Evans,  Porter,  and  Hara  and  Akechi 
that  oxide  layers  or  irregular  surfaces  can  create  inhomogeneous  necks  during 
the  very  initial  stage  of  sintering.  These  inhomogenities  may  be  oxide  islands 

or  pores  which  effectively  block  diffusive  mechanisms  leading  to  densification. 

10 

Addition  of  sintering  aids  such  as  B^C  or  C  to  SiC  contaminated  with  oxide 
layers  may  result  in  binary  or  ternary  liquid  phases  at  the  neck  surfaces 
rather  than  within  the  neck,  and  therefore  enhanced  sintering  behavior. 

This  paper  will  describe  the  sintering  of  pure  a-SiC  powders  activated 
by  a  plasma  during  the  early  stages  of  densification. 
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II.  EXPERIMENTAL  PROCEDURE 

Three  grades  of  commercial  SiC  were  investigated;  denoted  as  1200  mesh  black*, 
1200  mesh  green*,  and  ultra-fine*.  X-ray  diffraction  (XRD)  of  the  as-received 
powders  indicated  predominantly  a-SiC  with  6H  the  major  polytype  (Table  I) . 

The  particle  size  distribution  of  dispersed  powders  is  given  in  Table  II. 

A  substantial  fraction  of  each  is  submicron  with  the  ultra-fine  containing  the 
highest  fraction  of  fines  and  exhibiting  the  best  compaction  characteristics. 

Powder  compacts  were  prepared  by  first  adding  2-3%  polyvinyl-alcohol 
as  a  binder  and  aged  for  four  days  prior  to  pressing.  Compacts  1.25  cm  D  X 
0.64  cm  were  pressed  uniaxially  at  800  MPa  yielding  a  green  density  of  58-60% 
of  theoretical. 

The  plasma  activated  sintering  (PAS)  apparatus  is  shown  schematically 

in  Fig.  1.  A  graphite  susceptor  serves  as  sample  support  and  heat  source 

with  a  W-W26Re  thermocouple  passing  through  the  susceptor  base.  Alumina 
+ 

wool  ,  held  in  place  by  0.005"  Mo  foil,  serves  to  insulate  the  fused  quartz 
furnace  tube.  Both  end  plates  are  aluminum  and  the  bottom  plate  is  water 
cooled . 

With  a  pellet  in  place,  the  unit  was  evacuated  to  5  pm  and  held  overnight. 
Heating  slowly  to  500-601C  accomplished  binder  burnout  followed  by  a  more 
rapid  temperature  increase  (40-50C/min. )  to  1000C.  A  glow  discharge  was 
initiated  by  bleeding  in  Ar  which  when  coupled  with  the  field,  ignited  the 
plasma.  Following  a  5-6  min.  hold  at  1000C  with  discharge,  the  Ar  was  turned 
off  and  the  temperature  raised  rapidly  to  the  sintering  temperature,  1850-1950C. 
Total  elapsed  time  was  50  min.  from  room  temperature. 


*Lonza  Chemical  Company,  Fair  Lawn,  MJ 
+Saffil,  ICI  Americas,  Wilmington,  DE 
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At  the  sintering  temperature  a  second  discharge  was  initiated  lasting 
up  to  10  min.  Following  the  hold  at  the  sintering  temperature,  power  was 
turned  off  and  the  entire  assembly  allowed  to  cool  under  vacuum.  Due  to  the 
intense  radiation  at  the  sintering  temperature,  it  was  difficult  to  observe 
the  glow  discharge  directly.  However,  an  induced  potential  is  present  across 
the  top  plate  and  grounded  graphite  susceptor  and  monitoring  of  this  potential 
allows  for  detection  of  the  discharge. 

Macroshrinkage  was  determined  from  height  and  diameter  measurements  of 
the  pellets  and  microshrinkage  from  SEM  observations  of  fractured  surfaces. 

III.  RESULTS 

The  microstructure  of  the  black  1200  mesh  SiC  sintered  at  1900C  without 
and  with  plasma  activation  is  shown  in  Fig.  2(a)  and  2(b)  respectively.  With¬ 
out  the  plasma  the  particles  have  retained  the  as-received  angular  morphology. 
Although  macroshrinkage  (determined  as  change  in  pellet  dimensions)  was  0.2%, 
little  if  any  microshrinkage  occurred  as  evidenced  by  the  lack  of  necks  between 
particles.  In  contrast,  when  sintered  with  the  discharge  the  morphology  is 
altered  substantially.  The  sharp,  angular  particles  have  become  rounded  and 
smooth  and,  more  importantly,  substantial  neck  growth  has  occurred.  Measured 
macroshrinkage  has  increased  from  0.2  to  1,0%  and  local  x/a  ratios  are  in  the 
range  0.2  to  0.5. 

Increasing  the  sintering  time  without  plasma  activation  resulted  in 
crystallographic  faceting  as  shown  in  Fig.  3  for  1200  mesh  black  SiC.  The 
faceting  reduces  the  surface  free  energy  and  therefore  the  available  driving 
force  for  densification  and  no  microshrinkage  is  observed.  It  is  apparent  from 
Fig.  3  that  the  hexagonal  polytypes  are  the  predominant  phases  present  at  1900C 


When  the  sane  material  is  sintered  for  longer  periods  of  time  in  addition  to 
using  the  plasma  as  described  above,  the  faceting  is  suppressed  and,  as 
expected,  significant  microshrinkage  occurs  between  the  smooth,  rounded 
particles  (Fig.  4). 

The  1200  mesh  green  SiC  sintered  with  and  without  plasma  activation 
exhibits  microstructural  features  similar  to  that  of  the  black  SiC.  Some 
degree  of  crystallographic  faceting  has  occurred  and  necks  have  formed  without 
the  plasma  activation  (Fig.  5(a)).  The  1200  mesh  green  is  a  finer  material 
having  60%  of  particles  less  than  1.0  pm  diameter  and  therefore  would  be 
expected  to  sinter  more  than  the  coarser  black  material.  When  sintered  with 
plasma  activation,  (Fig.  5(b)),  the  morphology  was  altered  as  it  was  with 
the  1200  mesh  black.  Substantial  neck  growth  occurred  with  high  local  x/a 
ratios  (>0.4)  and  several  areas  where  smaller  particles  have  coalesced  with 
the  larger  particles.  Although  the  microshrinkage  is  high  in  localized  areas, 
little  macro shrinkage  occurred  as  would  be  expected  at  1900C. 

The  microstructure  of  the  ultra-fine  SiC  sintered  at  1900C  for  30  min. 
without  a  plasma  activation  period  exhibits  faceting  similar  to  that  in  the 
1200  mesh  black.  Fig.  3,  but  to  a  lesser  extent,  Fig.  6(a).  With  the  plasma 
activation,  surface  alteration  occurred,  Fig.  6(b),  with  substantial  neck 
growth  and  areas  of  high  raicroshrinkage. 

X-ray  diffraction  analyses  of  all  the  materials  sintered  with  and  without 
plasma  activation  indicated  mainly  a-SiC  polytypes,  i.e.,  15R,  6H,  and  4H. 

This  is  consistent  with  published  data  indicating  that  the  cubic  B-SiC,  3C, 

1 1 

becomes  unstable  at  T  >  1900C. 

IV.  DISCUSSION 

Exposure  to  a  glow  discharge  during  sintering  significantly  altered  the 
crystallographic  makeup  as  compared  to  the  material  sintered  without  the 
discharge.  XRD  of  the  material  sintered  without  the  discharge  indicated  an 
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increase  in  the  15R  and  6H  polytype  and  a  decrease  in  the  4H  and  3C.  When 
sintered  with  the  discharge,  the  15R  decreased  to  less  than  that  in  Table  I, 

3C  further  decreased  to  0%,  while  the  4H  and  6H  increased  to  a  greater  level 
than  in  Table  I.  Phase  transformations  in  SiC  are  not  well  understood  and 
many  depend  strongly  on  impurities  and  the  atmosphere.  In  this  work,  the 
presence  of  a  glow  discharge  would  primarily  affect  the  surface  indicating 
that  the  polytype  transformations  are  dependent  on  surface  activity  and 
structure . 

Silicon  carbide,  like  other  non-oxide  ceramics,  is  thermodynamically 
unstable  in  an  oxygen  containing  environment  with  surface  contamination  by 

V 

tenacious  oxide  layers  the  most  common  result.  Total  removal  of  such  films 
and  layers  is  difficult  but  necessary  if  normal  sintering  is  to  proceed 
leading  to  densification.  In  this  work,  the  as-received  SiC  powers  were 
sintered  without  a  glow  discharge  with  no  attempt  made  to  clean  the  surfaces 
or  destroy  the  oxide  contamination.  Very  little  shrinkage  occurred  even  on 
a  microscopic  level  with  the  exception  of  the  1200  mech  green  SiC  (nominally 
purer  than  the  black  SiC) .  This  is  consistent  with  observations  hy  others 
that  SiC  is  difficult  to  sinter  without  sintering  aids. 

Crystallographic  faceting  is  observed  in  the  sintered  powders  indicating 
significant  surface  or  vapor  transport  mechanisms.  It  also  suggests  an  incubation 
time  for  faceting  greater  than  3  minutes  at  1900C.  As  the  contaminated  powder 
compacts  are  heated  in  vacuum,  the  surface  energy  will  have  a  tendency  to  lower 
via  faceting  or  sintering,  i.e.  reduction  in  surface  free  energy.  Simultaneously, 
the  oxide  layers  experience  a  greater  tendency  towards  volatization  with 
increasing  temperature.  Without  the  glow  discharge,  the  volatization  will  be 
delayed  allowing  the  faceting  to  proceed,  hence  reducing  the  surface  energy 
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and  lowering  the  available  driving  force  for  densification.  Even  partial 
volatization  of  the  oxide  films  could  disrupt  normal  sintering  by  creating 
islands  of  oxides  within  small  necks  created  by  surface  transport.  Blockage 
of  the  transport  paths  leading  to  shrinkage  would  allow  further  lowering  of 
surface  energy  via  faceting. 

When  exposed  to  a  glow  discharge  at  a  temperature  lower  than  the  sintering 
temperature,  it  is  believed  that  the  surfaces  are  effectively  cleaned  by  the 
impinging  ions,  which  would  allow  volatization  of  the  oxide  layers.  The 
particle  surfaces  would  be  left  in  a  highly  active  state  with  a  considerable 
defect  density.  Adjacent  particles  could  then  lower  their  surface  energy  by 
surface  transport  phenomena  leading  to  neck  formation.  Surface  transport  is 
obviously  high  as  verified  by  the  rounding  of  the  particles.  The  resulting 
necks,  although  not  sufficient  to  cause  densification,  do  provide  homogeneous 
regions  for  bulk  diffusive  mechanisms,  volume  or  grain-boundary  diffusion, 
which  lead  to  densification. 

V.  CONCLUSIONS 

Plasma  activated  sintering  (PAS)  has  been  shown  to  provide  apparently 
clean,  active  surfaces  in  SiC  powder  compacts.  These  active  surfaces  have 
a  higher  tendency  towards  surface  energy  reduction  via  diffusive  sintering, 
leading  to  densification,  than  powders  sintered  without  a  glow  discharge. 

Although  experiments  have  been  conducted  only  with  low  density  compacts 
to  eliminate  plasma  permeability  problems,  it  is  expected  that  higher  green 
densities  could  be  used. 
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TABLE  II.  PARTICLE  SIZE  DISTRIBUTION 


Size 

Range 

1200  Mesh 
(Black) 

1200  Mesh 
(Green) 

Ultra-fine 

1-3  ym 

>1  um 

22% 

59% 

65% 

1-5 

65 

36 

34 

5-10 

10 

4 

>1 

1C-15 

2 

>i 

>1 

>15 

>1 

>i 

>1 
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TABLE  I. 

X-RAY  DIFFRACTION 

ANALYSIS  OF 

SiC  POWDERS* 

SiC 

1200  Mesh 

1200  Mesh 

Ultra-fine 

Polytype 

(Black) 

(Green) 

1-3  ym 

a-SiC 

15R 

10  w/o 

5  w/o 

13  w/o 

6H 

71 

81 

74 

4H 

10 

5 

11 

3-SiC 

3C 

9 

9 

^2 

*Lonza  Chemical  Co.,  Fairlawn,  NJ 
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Fig.  1 


Schematic  of  plasma  activated  sintering  apparatus. 


Fig.  2 


Fig.  3 


Fig.  4 


Fig.  5 


Fig.  6 


Microstructure  of  1200  mesh  black  SiC  sintered  in  vacuum  for 
3  min.  at  1900C;  (A)  without  plasma  activation,  and  CB )  with 
plasma  activation  for  6  min.  at  1000C  and  3  min.  at  1900C. 

Microstructure  of  1200  mesh  black  SiC  sintered  in  vacuum  for 

9  min.  at  1900C  without  plasma.  Note  the  prevalence  of  hexagonal 
morphology , 

Microstructure  of  1200  mesh  black  SiC  sintered  in  vacuum  for 
30  min.  at  1900C  with  plasma  activation.  Faceted  surfaces  are 
not  present  and  significant  morphological  alteration  has  occurred. 
Microstructure  of  1200  mesh  green  SiC  sintered  in  vacuum;  (4) 
for  6  min.  at  1895C  without  plasma  activation,  and,  CB )  for 

10  min.  at  1900C  with  plasma  activation. 

Microstructure  of  ultra-fine  SiC  sintered  in  vacuum  for  30  min. 
at  1900C;  (A)  without  plasma  activation,  and,  (b)  with  plasma  activation. 
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FIG.  1  Plasma  ^intoririK  Apparatus 
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SIC  SINTERED  3  MIN.  AT  1900C  WITH  Ar  GLOW  DISCHARGE 


FIG.  2(B) 


FIG.  3 


SIC  (1200  Mesh  Black)  SINTERED  9  MIN.  AT  1900C 
WITHOUT  GLOW  DISCHARGE 
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fig.  5(A)  SIC  (1200  Mesh  Green)  SINTERED  6  MIN.  AT  1895C 
WITHOUT  GLOW  DISCHARGE 


FIG.  5(B)  SIC  SINTERED  10  MIN.  AT  1900C  WITH  Ar  GLOW  DISCHARGE 
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ABSTRACT 


Rubbery  and  glassy  network  polymers  have  been  prepared  and 
characterized  by  various  physical  and  chemical  methods.  Fracture 
studies  have  been  undertaken  on  these  polymers  using  a  range  of 
loading  rates  and  temperatures.  Initial  efforts  to  describe 
local  inhomogeneities  in  cured  epoxy  resins  by  electron  micros¬ 
copy  have  been  reasonably  successful.  The  feasibility  of  model 
inhomogeneous  glassy  network  polymers  with  precisely  controlled 
structure  has  been  established,  with  the  intent  of  comparing  the 
fracture  behavior  of  these  materials  with  analogous  epoxy 
structures.  Also,  a  series  of  end- linked  polybutadienes  and 
poly(butadiene-styrene)  copolymer  networks  have  been  prepared 
and  characterized.  These  polymers  provide  known  variations  in 
cross-link  density  and  backbone  flexibility.  Viscoelastic 
properties  of  some  of  the  rubbery  polymers  have  been  determined 
at  both  the  University  of  Akron  and  the  University  of  Pittsburgh 
in  dynamic  and  creep  experiments.  Correlation  of  mechanical 
loss  with  the  time  dependent  fracture  energy  in  each  material 
system  is  underway. 
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I.  EPOXY  NETWORKS 
Introduction 

The  objective  of  this  research  is  the  quantitative  correlation 
of  possible  microphase  heterogeneity  with  time-dependent  behavior, 
including  fracture.  This  investigation  has  followed  several  paths, 
all  of  which  are  linked  by  a  unifying  concept;  a  relationship 
^ists  between  microstructural  variables,  mechanical  loss  character¬ 
istics  ,#and  time- dependent  fracture  energy.  Most  of  the  effort  has 
been  spent  developing  the  diverse  techniques  and  analyses  required 
for  a  successful  completion  of  this  project. 

Network  Synthesis 

Variations  in  reactant  stoichiometry  and  process  history  are 
expected  to  produce  measurable  variations  in  structural  hetero¬ 
geneity.  Considerable  latitude  in  thermal  processing  variation  is 
possible  using  systems  with  low  initial  viscosities  and  long  gel 
times . 

Since  fracture  and  fatigue  properties  are  particularly  sensitive 
to  specimen  defects,  the  production  of  void-free  castings  has  been 
a  high  priority  objective.  It  has  been  found  that  the  degree  of 
care  required  to  produce  defect- free  samples  is  determined  mostly 
by  the  reactant  system's  viscosity. 

Aliphatic  amine  cured  systems,  such  as  triethylene  tetramine/ 
Epon  828,  gel  in  a  relatively  short  period  of  time  at  room 
temperature.  These  systems  have  viscosities  greater  than  50  poise 


hobsono  Fios  bumc-not  nuo 


Vv/ 


2 


at  room  temperature.  In  addition,  the  viscosity  increases  as 
the  reaction  proceeds  to  the  gel  point.  Entrapped  air  produced 
during  reactant  mixing  must  be  removed  under  vacuum.  After 
debulking,  simply  pouring  the  mixture  into  a  mold  appears  to 
produce  void-free  castings.  However,  inspection  of  these  castings 
with  a  polarized  light  microscope  reveals  that  a  large  number  of 
voids  are  present  which  cannot  be  detected  with  the  unaided  eye. 

It  has  been  found  that  these  defects  can  be  eliminated  by  mold 
filling  under  vacuum.  This  extra  care  would  be  unnecessary  if  the 
viscosity  could  be  lowered  by  increasing  the  reaction  temperature; 
however,  the  pot  life  is  considered  to  be  too  short  even  at  room 
temperature.  These  same  systems  can  be  handled  quite  easily 
without  any  special  precautions  (at  room  temperature)  if  the 
mixture  contains  a  non-reactive  diluent  such  as  dimethyl  sulfoxide 
Void-free  pre-swollen  castings  have  been  prepared  by  pouring  the 
diluted  reaction  mix  into  molds,  without  the  need  of  any  vacuum. 

In  contrast,  the  aromatic  amine  cured  systems,  such  as 
methylene  dianiline/Epon  828,  have  long  pot-life  at  elevated 
temperatures  and  are  therefore  particularly  easy  to  work  with. 

No  special  precautions  are  required  to  produce  void-free 
specimens  with  these  low  viscosity  mixtures.  With  the  above 
combination,  for  example,  the  reactants  are  mixed  at  60°C  and 
poured  directly  into  preheated  molds.  Air  entrapped  during  the 
pouring  step  quickly  dissipates. 

Two  closed  mold  configurations  have  been  used.  One  has 
permitted  the  casting,  in  vacuum,  of  flat  plates  of  dimensions 
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4"  x  6"  x  1/4"  using  the  vertical  position  shown  in  Figure  1. 

The  second  mold,  shown  in  Figure  2,  allows  the  casting  of  four 
compact  tension  specimens  at  one  time.  Both  molds  are  teflon 
coated.  In  addition,  the  white  posts  seen  in  Figure  2  are  produced 
from  teflon  rod. 

Network  Characterization 

The  average  molecular  weight  between  crosslink  points,  Mc, 
is  typically  evaluated  in  elastomers  via  stress-strain  data 
plotted  in  a  Mooney-Rivlin  fashion.1  Since  epoxy  resins  have 
high  glass  transition  temperatures,  and  this  method  requires  data 
obtained  in  the  "rubbery"  regime,  network  degradation  is  of  some 
concern.  A  significant  achievement  has  been  the  preparation  of 
epoxy  networks  in  the  presence  of  a  high  boiling  diluent,  resulting 
in  a  pre-swollen  material  with  Tg  near  room  temperature.  A  typical 
Mooney-Rivlin  plot  of  an  epoxy  network  synthesized  in  the  presence 
of  dimethyl  sulfoxide  is  seen  in  Figure  3. 

The  Mooney-Rivlin  equation  is  based  on  the  rubber  elasticity 
theory  for  networks  composed  of  statistically  distributed  and 
freely  rotating  Gaussian  chains.  In  the  example  shown  in 
Figure  3,  the  Mc  value  of  ^615  g/mole,  corresponds  to  two  monomer 
segments  of  DGEBA  between  crosslinks,  or  less  than  30  main  chain 
atoms  between  junction  points.  Over  40%  of  the  monomer  molecular 
weight  resides  in  the  bulky  benzene  rings,  which  reduce  network 
chain  flexibility  and  rotation. 
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Rubber  elasticity  theory  cannot  be  expected  to  apply  to  such 
a  tight  network  composed  of  short  strands  with  low  flexibility. 

Yet,  the  Mc  shown  in  this  example  is  not  too  far  removed  from  the 
theoretically  determined  Mc  of  ^380  g/mole. 

Tg  is  being  measured  via  DSC  and  dilatometry.  A  typical 
dilatometric  curve  is  seen  in  Figure  4  for  a  network  synthesized 
in  the  presence  of  dimethyl  formamide. 

Structural  Heterogeneity  Characterization 

Several  investigators  have  shown  that  a  characteristic  nodular 
morphology  can  be  seen  on  carbon/platinum  self-shadowed  replicas 
of  epoxy  network  fracture  surfaces  and  solvent  etched  surfaces. 

The  nodules  are  thought  to  be  intrinsic  structural  features  as 
their  size  was  found  to  vary  with  reactant  stoichiometry.  It  has 
been  proposed  that  the  nodules  are  sites  of  high  crosslink  density 
embedded  in  a  matrix  of  lower  crosslink  density  material.  Koutsky 
and  Mijovic2  have  shown  a  correlation  between  nodule  size,  extent 
of  plastic  deformation,  and. fracture  energy;  thus  there  is  reason 
to  explore  this  relationship  further. 

The  technique  of  self- shadowed  surface  replication  has  been 
advanced  to  a  high  degree,  as  is  evidenced  by  micrographs  shown 
later  as  Plates  1- 

Work  is  in  progress  examining  the  rate  sensitivity  of  nodular 
size  obtained  from  etched  fracture  surfaces. 

It  is  not  clear  at  this  time  whether  nodules  are  intrinsic 
structural  features  or  a  result  of  the  technique  used  to  obtain 
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micrographs  of  them.  For  example,  Mijovic4  prepared  samples  of 
aliphatic  amine  cured  epoxies  varying  in  stoichiometry  on  the 
epoxy  excess  side,  and  subsequently  determined  their  fracture 
energies  and  resulting  morphologies  at  one  temperature  and  cross¬ 
head  speed.  Networks  formed  with  different  stoichiometric 
reactant  ratios  will  have  different  crosslink  densities  and  glass 
transition  temperatures,  and  therefore  different  rate  sensitivities 
and  abilities  to  deform  under  applied  loads.  Herein  lies  one  as 
yet  unexplored  possibility  for  a  nodule  size- fracture  energy 
relationship. 

We  intend  to  examine  the  fracture  surface  morphology  of 
networks  differing  in  their  chemical  constitution  and/or  thermal 
processing  history  via  the  replication  technique.  Measurable 
structure  variations  will  be  correlated  with  time-dependent 
mechanical  properties. 

Time-Dependent  Fracture  Characterization 

Fracture  energy  is  being  examined  over  a  wide  temperature 
range  using  single  edge  notched  rectangular  tensile  specimens. 

The  fracture  behavior  is  being  examined  within  the  framework  of 
Andrews’  generalized  fracture  theory,3  wherein  the  cohesive 
fracture  energy,  T,  per  unit  area  of  fracture  surface  is  given  by 

T  -  kcW0ic  CD 

where  k  is  a  parameter  dependent  on  specimen  dimensions,  strain 
level,  and  crack  orientation,  c  is  the  crack  length,  and  W0jC  is 
the  critical  input  strain  energy  density.  A  linear  relation  is 
predicted  between  W0iC  and  c* 1 .  A  plot  of  this  relation  for  the 
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1/1  Epon  828/MDA  specimen,  whose  fracture  topography  is  seen  in 
the  enclosed  micrographs,  is  shown  in  Figure  5.  The  slope  of  the 
best  fit  line  yields  a  value  for  T  of  345  J/mz . 

The  behavior  of  this  material  at  different  test  temperatures 
and  the  same  composition,  but  of  different  thermal  history, 
is  currently  under  investigation. 

Fracture  specimens  have  been  prepared  by  casting  in  an  open 
mold  using  a  teflon- coated  plate  and  an  o-ring  dam  of  about  6" 
diameter.  Samples  are  cut  from  the  sheets  obtained  and  edge 
notches  inserted  at  temperatures  above  the  sample  Tg ,  using  the 
sliding  fixture  shown  in  Figure  6. 

Mechanical  Loss  Characterization 

Attempts  to  characterize  loss  modulus  and  tan  6  via  the 
Rheometrics  Mechanical  Spectrometer  using  an  eccentric  rotating 
cantilever  beam  (ERCB)  experiment  have  been  unsuccessful.  We 
are  unable  to  measure  the  low  forces  associated  with  the  loss 
modulus  in  this  low  energy  dissipating  material.  We  are  looking 
into  the  possibility  of  a  torsional  mode  in  conjunction  with  an 
instrument  available  through  the  Lord  Corporation.  A  master 
curve  of  the  storage  modulus  of  a  1/1  Epon  828/MDA  sample  obtained 
via  the  ERCB  mode  is  shown  in  Figure  8. 

To  facilitate  data  reduction  in  this  program  a  computer 
routine  has  been  written  which  shifts  any  viscoelastic  function 
via  the  WLF  transform.  Figures  7  and  8  were  obtained  using  this 
program. 
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Fatigue  Crack  Propagation 

Compact  tension  specimens  are  being  stockpiled  and  testing 
should  be  started  within  the  month. 

K'e  intend  to  examine  the  sensitivity  of  fatigue  behavior 
toward  epoxy  network  thermal  processing  variations. 


II.  MODEL  INHOMOGENEOUS  GLASSY  NETWORKS 


A  method  of  dispersing  preformed  cross-linked  particles  of 
known  size  and  spherical  shape  (gel  balls)  in  a  polymer  matrix, 
and  of  molding  the  resulting  mixture  has  been  found.  Particles 
may  be  purchased  or  produced  from  emulsion  polymerized  latex 
using,  for  example,  about  II  divinyl  benzene  copolymerized  with 
styrene.  This  discovery  makes  possible  a  study  of  the  properties 
of  model  systems  containing  nodular  structures  (ca.  10-200  nm 
diameter)  like  those  believed  to  be  present  in  epoxy  resins. 

For  example,  the  effect  of  "gel  balls"  of  known  size  and  content 
on  physical  properties  of  rigid  polymers  can  be  determined.  The 
morphology  of  systems  known  to  contain  these  particles  can  be 
examined.  Techniques  believed  to  make  the  dispersed  gelled  particl 
visible  can  be  established.  These  studies  should  lead  to  added 
understanding  of  the  epoxy  resins  themselves. 

The  experimental  procedures  have  been  evolved  through 
considerable  trial-and-error  experimentation.  Simply  mixing  the 
gelled  particles  with  linear  polymer  was  unsuccessful.  Molded 
disks  were  inhomogeneous  and  not  well-fused.  Absorbed  monomer 
containing  initiator  in  the  particles  was  polymerized  before 
molding.  This  led  to  molded  disks  that  appeared  to  be  well- fused 
both  to  the  eye  and  when  examined  through  a  polarizing  film.  The 
modified  particles  could  be  mixed  with  pure  polymer  and  the 
mixture  could  then  be  molded  to  give  disks  with  variable  gel 
particle  content.  Some  pre-polymerization  of  the  monomer  is 
desirable.  Molding  in  the  presence  of  absorbed  monomer  was  an 
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improvement  over  simple  mixing  of  the  particles  with  polymer. 
However,  control  of  variables  such  as  particle  concentration, 
molecular  weight  of  polymer,  and  distribution  of  particles  in 
the  polymer  matrix  was  comparatively  poor.  Details  of  experi¬ 
ments  in  which  the  particles  were  either  Styragel  or  Merrif ield 
resin  and  the  polymer  matrix  was  po lystyrene  while  the  monomer 
was  styrene  follow. 

Materials : 

Polystyrene :  Dow  Chemical's  Styron  678  -  26  Clear. 

Merrifield  Peptide  Resin:  Pierce  No.  24610,  1%  crosslinked  - 
chloromethylated  styrene-divinylbenzene  resin,  chloromethylation 
level  -  1-2  meq./g  resin,  <44-75  microns. 

Styragel  Type  500  %:  Part  No.  27126,  Column  Packing  material. 
Batch  1632,  37-75  microns. 

Styrene :  Eastman  1454,  Lot  No.  J6A. 

Procedure : 

The  polystyrene  was  dissolved  in  benzene  and  reprecipitated 
from  methanol.  After  washing  with  methanol  it  was  dried  and 
degassed  in  a  vacuum  oven  at  room  temperature.  The  resins  were 
washed  with  hexane  and  petroleum  ether  and  dried  and  degassed 
along  with  the  polystyrene.  After  drying  the  polystyrene  and 
resins  were  kept  under  nitrogen  until  used.  The  styrene  was 
distilled  under  nitrogen  before  use.  B.P.  -  36°C/8mm.  Purified 
AIBN  (azo-bis - isobutyronitrile)  obtained  from  C.-C.  Kuo  was  used 
as  received. 
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Modifications  of  the  resin  surfaces  were  carried  out  under 
nitrogen  in  screw  cap  vials  with  teflon  liners.  All  manipulations 
were  done  in  a  glove  bag  filled  with  nitrogen.  The  mixtures  in 
each  vial  are  listed  in  Table  1.  A  vortex  stirrer  was  used  to  mix 
thoroughly.  When  the  styrene  monomer  was  polymerized  before 
molding,  the  vials  were  heated  in  an  oven  at  12S°C  for  2.5  hours 
and  kept  at  room  temperature  for  three  days  before  use. 

Molding  was  done  in  a  one- inch  compression  mold  using  about 
three  grams  of  sample.  Two  Pasadena  presses,  one  at  160°C,  and 
the  other  with  the  platens  filled  with  continuously  circulating 
cold  water,  were  used  for  molding.  The  materials  to  be  molded 
were  first  premolded  between  teflon  coated  aluminum  foil  at  a 
pressure  of  5  for  30  sec.  (All  pressures  are  on  a  5-inch  ram.)  A 
typical  molding  cycle  consisted  of  10  minute  preheat  of  the  mold, 
adding  of  appropriate  pieces  of  the  premolded  material,  pressure 
up  to  indicated  10  for  10  seconds  and  release  three  times  to 
minimize  bubble  formation,  molding  at  160°C  for  two  minutes  at 
a  pressure  indication  of  5,  rapid  transfer  to  the  cold  press, 
and  cooling  at  the  same  pressure  for  eight  minutes.  The  modi¬ 
fied  resins  were  free  flowing  powders  and  could  be  used  without 
premolding.  The  results  of  the  molding  experiments  are  summarized 
in  Table  2. 

Electronmicrographs  (Figure  9)  of  thin  sections  of  the  molded 
specimens  prepared  from  the  samples  in  which  the  styrene  was 
polymerized  in  situ  indicated  that  the  gelled  particles  of  both 
the  Styragel  and  Merrifield  Resins  are  clearly  distinguishable 
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from  the  polystyrene  matrix.  The  micrographs  shown  in  Figure  9 
were  taken  by  transmission  electronmicroscopy .  Other  micro¬ 
graphs  of  the  above  resins  and  of  another  one,  presumably  with 
smaller  particle  size,  were  taken  by  scanning  electronmicroscopy. 
The  purposes  of  these  studies  were  twofold: 

1.  to  verify  that  the  sizes  of  the  resins  were  as 
stated  by  the  manufacturers  and  to  get  some  idea 
of  the  particle  size  distribution;  and 

2.  to  determine  if  the  particles  swelled  significantly 
during  polymerization. 

From  the  photographs  taken  so  far  we  are  able  to  conclude  that 
indeed  the  average  size  of  the  new  resin  particles,  stated  to  be 
<  37p ,  are  smaller. than  those  of  the  other  two  resins,  which  are 
stated  to  be  in  the  37- 75y -range.  The  distribution  of  particle 
sizes  of  the  Merrifield  Resin  appears  to  be  lower  than  that  of 
the  Styragels,  which  are  clearly  not  uniform. 

Scanning  electronmicrographs  of  the  surfaces  of  fractured 
molded  specimens  are  planned  for  the  future.  We  are  awaiting  a 
specimen  holder  that  will  allow  us  to  calibrate  the  scope  and 
make  quantitative  size  comparisons  possible.  Alternately,  or 
in  addition,  a  Zeiss  Particle  Analyzer  available  in  our 
Institute  will  be  used  to  study  particle  size  and  particle  size 
distribution. 

Other  smaller  sizes  of  gelled  particles  are  being  purchased 
and  still  others  will  be  synthesized. 
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III.  RUBBERY  NETWORKS 

The  objectives  of  this  task  are  two- fold:  to  prepare  a 
viscoelastic  network  polymer  applicable  to  a  variety  of 
engineering  studies,  and  to  correlate  time-dependent  fracture 
energy  with  molecular  structure.  The  networks  are  prepared  with 
varying  average  molecular  weight  between  cross-links  (Mc) ,  and 
have  two  different  backbone  structures.  While  the  early  phases  of 
this  task  are  examining  molecular  effects  on  fracture  energy  in 
rubbery  networks,  the  fracture  behavior  of  glassy  network  polymers 
such  as  epoxy  resins  are  also  being  studied  above  their  Tg,  if 
such  networks  can  be  characterized  quantitatively  by  methods 
indicated  in  section  I  above. 

Two  prepolymers  from  ARCO  have  been  selected  for  early 
study.  ARCO’s  R-45HT  is  a  dihydroxy-terminated  liquid  poly¬ 
butadiene  of  reported  average  molecular  weight  3000.  It  is 
colorless,  noncrystallizing  and  commercially  available  in  large 
quantities.  ARCO's  CS-15  is  a  dihydroxy- terminated  liquid 
styrene-butadiene  copolymer  of  reported  average  molecular  weight 
3400.  It  is  noncrystallizing,  but  a  little  hazy  and  available 
in  limited  quantities. 

Both  tetrafunctional  and  trifunctional  curing  agents  were 
tried  in  preparing  the  networks;  however,  only  the  trifunctional 
isocyanate  linked  structures  succeeded.  Tetrafunctional  curing 
agents  such  as  1 , 2 , 4 , 5 -benzenetetracarboxylic  anhydride  and 
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1 , 4 , 5 , 8-naphthalenetetracarboxylic  anhydride  are  both  high  melting 
materials  and  not  soluble  in  the  prepolymer.  These  drawbacks  make 
the  efficient  mixing  of  these  materials  practically  impossible. 

The  trifunctional  curing  system,  toluene  diisocyanate  and 
2, -ethyl- 2- (hydroxymethyl) -1 , 3-propanediol ,  consists  of  low 
melting  reactants.  The  curing  reaction  of  the  mixture  is  slow 
enough  to  allow  time  for  good  mixing.  Chain  extended  networks  of 
this  type  have  also  been  prepared  using  1,4-butane  diol. 

Several  samples  of  these  trifunctional  networks  have  been 
characterized.  The  measurements  of  their  fracture  energies  have 
been  started,  and  some  of  the  preliminary  results  will  be  given 
later  in  this  report. 

Characterization 

Ten  samples  prepared  by  casting  on  a  teflon  plate  have  been 
characterized  to  be  used  in  fracture  studies.  Table  3  shows  the 
Mc,  Tg  and  sol  fraction  characteristics  of  the  replicate  castings. 
The  Mooney-Rivl in  parameter,  C2 ,  is  indicated  as  well.  Series 
TBJX  was  prepared  from  R45-HT;  series  TBzX  and  TB3X  are  chain- 
extended  R45-HT.  A  second  set  of  samples  of  similar  stoichio¬ 
metries  was  prepared  and  characterized  using  a  different  batch 
of  ARCO  R45-HT.  These  are  shown  in  Table  4.  Sample  TSi  is  an 
isocyanate  cross-linked  CS-15;  another  sample  designated  as  PBi 
is  a  high  molecular  weight  polybutadiene  cured  with  dicumyl- 
peroxide,  to  be  used  for  comparison  purposes. 

The  number  average  molecular  weight  between  cross-links,  Mc, 
was  computed  from  the  Ci  intercept  values  taken  from 
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'^^FTOpS  . 

Sampl^"'^ 

mm 

IHH 

HR9 

%  Sol 
Fraction 

TBn 

3,340 

2.72 

KOI 

0.8 

TBl2 

3,400 

B 

-79.6 

tb13 

3,400 

2.82 

-78. S 

tb14 

3,300 

2.65 

TB2  1 

5,200 

' 

2.55 

-79.0 

1.5 

TB  2  2 

5,170 

2.70 

-79.0 

tb23 

5,800 

2.80 

tb31 

8,300 

2.0 

-79.5 

3.0 

■m 

1 

1.52 

-79.0 

tb33 

9,700 

1.30 

Table  3 

Characterization  of  Networks 
for  Fracture  Studies 


a  :  from  Mooney-Rivl in  equilibrium  tensile  plot 
b  :  measured  by  DuPont  DSC  (Model  990)  at  5°C/min. 
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C  xl  0 '  5 
m2 


C2  xl 0 ‘ 5 


TB  i 

tb2 

4.25 

2.2 

1.30 

3.0 

2,750 

5,250 

3,100 

6,200 

-80 

-80 

0.954 

0.945 

32.9 

32.9 

33.4 

35.6 

PBi 


1.30 


3,700  8,300 


14.30 


15.0 


Table  4 

Characterization  of  Samples  TBj,  TB2,  TSj  and  PBj 
a  :  From  Mooney- Rivl in  plot 

b  :  Calculated  data  from  reported  prepolymer  Mn. 

c  :  Measured  by  DuPont  DSC  Model  990  at  the  speed 
of  10°C/min. 

d  :  Determined  by  pycnometer,  ASTM  D297. 
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equilibrium  stress-strain  Mooney-Rivlin  plots,  as  shown  in 
Figure  11.  Somewhat  surprisingly,  the  Mc  values  calculated  from 
this  plot  match  very  well  with  those  of  their  prepolymers  deter¬ 
mined  from  the  chemical  composition  (Table  4) .  Although  the 
average  hydroxyl  functionality  of  the  prepolymer  is  given  as 
slightly  greater  than  2,  it  is  well-known  that  the  functionality 
distribution  is  rather  large.  The  Mooney-Rivlin  plots  of  TB2 
swelled  in  hexane  still  had  a  noticeable  C2  value  (Figure  12) 
and  the  Cj  values  are  larger  than  those  determined  in  unswollen 
samples  by  approximately  201. 

The  true  stress  vs.  strain  plots  for  the  more  tightly  linked 
samples  show  linear  behavior  to  relatively  high  strains  (X^1.8), 
as  shown  in  Figure  13.  The  initial  modulus  of  rubber-like 
materials  should  be  between  6  Ci  and  6(Ci+C2).5  The  values  obtained 
from  these  samples  indicate  that  E0  is  almost  equal  to  6(Ci+C2) 
(Table  4) . 

Sample  TBi  is  currently  the  best  candidate  for  a  standard 
viscoelastic  network  polymer.  It  can  be  made  quite  transparent 
by  casting  on  Ferro  type  plates  instead  of  teflon-coated  plates 
used  frequently  in  these  studies.  It  is  linearly  viscoelastic 
below  extensions  of  about  15%  (Figure  14),  The  loss  modulus 
appears  to  be  independent  of  the  prestrain  up  to  strains  of  0.4 
(Figure  15)  while  tan  <5  increases  somewhat  linearly  with  strain 
(Figure  16).  Some  preliminary  fracture  studies  of  the  series 
TBii-TBm  also  indicate  that  it  is  rate  dependent  up  to  70°C. 
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Data  accumulated  thus  far  indicate  that  the  trifunctional 
end-linked  networks  behavior  similarly  with  randomly  cured 
rubbers.  However,  future  studies  over  a  wider  range  of  material 
variables  and  test  conditions  will  be  required  before  this 
conclusion  may  be  generalized. 


Fracture  Studies 

Samples  of  R45-HT  have  been  prepared  in  a  vertical  mold, 
providing  uniform  thickness,  with  and  without  chain  extender. 
Samples  TB11,  TB23,  and  TB3j  have  been  characterized  to  have  Mc 
values  of  3,340,  5,800  and  9,700,  respectively. 

Figure  17  shows  the  tear  energy  of  TBii  at  different 
temperatures  and  crosshead  rates.  The  WLF  procedure  is  used 
to  shift  the  data  to  produce  a  master  curve  using  the  following 
constants : 


log  aj 


-  31 . 1 (T-Tp) 
78.0  ♦  T-Tg 


A  preliminary  master  curve  of  the  tear  energies  of  TBii,  TB23  and 
TB33  is  shown  in  Figure  18.  It  may  be  seen  that  at  the  higher 
speeds  or  lower  temperatures,  the  tear  energy  curves  of  TBi!, 

TB23  and  TB33  separate  quite  a  lot,  while  at  low  speed  and  high 
temperature  they  converge. 
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Discussion 

The  major  efforts  in  this  task  during  the  first  year  have 
been  to  develop  convenient  methods  for  uniform  sample  preparation, 
refine  the  characterization  methods  for  end-linked  elastomers 
and  to  begin  tear  studies  on  the  networks  produced.  A  limited 
amount  of  viscoelastic  data  has  been  obtained  in  dynamic  and  stress 
relaxation  testing.  Future  plans  are  to  supply  samples  of  each 
of  the  end-linked  networks  to  Dr.  Plazek  of  The  University  of 
Pittsburgh  for  his  specific  studies  of  the  material  response. 
Various  viscoelastic  quantities  will  be  examined,  in  particular 
the  mechanical  loss  (G"  and  tan  6)  correlation  with  the  time- 
dependent  fracture  behavior  of  these  materials.  Several  samples 
from  earlier  networks  prepared  in  the  course  of  the  investigation 
have  been  tested  in  the  torsional  creep  mode  by  Dr.  Plazek.  The 
following  section  on  thermal  expansion  studies  raises  important 
questions  about  transitions  which  may  be  related  to  the  relatively 
high  concentration  of  urethane  links  in  the  various  networks. 
Subsequent  investigations  to  provide  information  on  the 
relationship  of  structure  to  volume  relaxation  or  phase- 
transitions  are  underway. 


251 


22 


IV.  Thermal  Expansion  Measurements 

Dilatometric  studies  were  started  on  the  polyurethane 
network  polymers  prepared  from  R-45HT  and  CS-15.  DSC  measure¬ 
ments  indicate  that  the  glass  transition  temperature  of  these 
polymers  O  -90°C)  is  well  below  the  freezing  point  of  mercury 
(-39°C)  and  of  most  other  suitable  confining  liquids.  Therefore, 
a  slightly  modified  thermal  expansion  apparatus  previously 
described  by  M.  L.  Dannis6*7  was  used  for  our  initial  studies. 
Since  no  confining  liquid  is  required,  this  apparatus  can  be 
used  in  any  desired  temperature  range  provided  a  suitable  cooling 
or  heating  chamber  is  available. 

Apparatus 

As  shown  in  Figure  19,  the  apparatus  is  a  simplified  form 
of  dilatometer  in  which  the  actual  expansion  or  contraction  of 
the  sample  is  measured  mechanically  in  one  chosen  direction  or 
dimension.  A  mirror  tripod  rests  one  leg  on  the  sample  and  the 
other  two  on  a  metal  cup  holding  the  sample.  A  circular  aluminum 
plate  beneath  the  leg  on  the  sample  prevents  indentation  of  the 
^  sample.  Movement  of  the  sample  relative  to  the  cup  tips  the 

tripod  and  mirror  and  thus  moves  a  spot  obtained  by  reflecting  a 
beam  of  light  from  the  mirror  onto  a  scale.  The  cells  were 

loaned  to  us  by  M.  L.  Dannis,  whose  advice  in  assembling  the 

I 

apparatus  is  gratefully  acknowledged. 

Temperature  of  the  sample  is  measured  by  an  iron- cons tantan 
thermocouple  inserted  close  to  the  midpoint  of  the  sample  and 
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connected  to  a  potentiometer.  The  thermocouple  is  made  of  No.  40 
gauge  wire  to  eliminate  false  movements  of  the  sample  arising 
from  the  springiness  of  heavier  gauge  wires.  An  ice  junction 
is  used.  Values  of  temperature  are  read  from  standard  tables 
relating  the  observed  millivolts  to  temperature. 

For  measurements  below  room  temperature,  the  assembled  cup 
containing  sample  and  tripod  is  flooded  with  liquid  nitrogen  in 
a  vacuum  flask.  After  the  excess  nitrogen  has  boiled  away, 
expansion  of  the  sample  on  warming  to  room  temperature  is 
followed  as  a  function  of  temperature.  The  rate  of  warming  near 
-100°C  is  about  0.3°C/min;  near  0°C  it  is  less  than  0.1°C/min. 

For  measurements  above  room  temperature  the  vacuum  chamber 
is  replaced  by  a  tubular  heater,  whose  temperature  is  increased 
by  gradually  increasing  the  applied  voltage  with  the  aid  of  a 
Variac.  The  heating  rate  is  about  0.25°C/min. 

Results 

Figure  20  shows  a  typical  differential  thermal  expansion 
curve  for  a  polybutadiene  containing  954  1 , 4-cis  units.  This 
figure  is  taken  from  data  reported  by  M.  L.  Dannis7  and  is 
included  as  an  aid  to  interpretation  of  the  data  presented 
below. 

Figure  21  shows  some  results  obtained  on  our  apparatus 
using  samples  of  two  well-known  polymers.  Firestone's  Diene 
35NFA  is  a  noncrystallizing  polybutadiene  containing  554  1 , 4-cis , 
354  1 , 4-trans  and  104  vinyl  units  and  was  cured  using  0.14 
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dicumyl  peroxide,  Hercules  Dicup.  The  data  obtained  with  this 
sample  shows  only  the  expected  Tg  at  approximately  -96.5°C. 

The  Pale  Crepe  sample  was  used  as  a  gum  rubber  (uncured)  and  the 
plot  obtained  using  it  shows  the  glass  transition  temperature, 
some  crystallization  and  the  onset  of  melting. 

Figure  22  shows  typical  results  obtained  using  sample 
TB24,  a  chain-extended  polyurethane  from  R-45HT  and  1,4- 
butanediol  having  Mc=6000.  Similar  data  have  been  obtained 
on  three  different  samples  of  the  same  material  and  on  repeated 
runs  on  the  same  sample.  The  curve  shows  some  unexpected 
features.  At  about  -90°C,  the  glass  transition  temperature 
expected  from  DSC,  a  contraction  occurs.  By  analogy  to  the 
results  given  in  Figure  20,  this  contraction  may  be  associated 
with  crystallization.  The  rest  of  the  curve  is  consistent  with 
such  an  interpretation.  However,  crystallization  of  this  sample 
was  not  expected  and  other  possible  interpretations  such  as  phase 
segregation  or  mechanical  relaxation  are  being  considered.  Similar 
measurements  with  the  other  samples  in  the  series  are  in  progress. 
The  results  on  the  polymer  prepared  from  the  copolymer  containing 
styrene  should  be  especially  significant  but  are  not  yet 
completed.  This  sample  should  not  crystallize. 

Discussion 

It  is  important  for  our  studies  to  establish  that  the 
transition  observed  at  ^  -90°C  is  real  and  if  it  is,  to  determine 
the  reason  for  the  transition.  Phenomena  such  as  crystallization 
or  crystallizability  have  very  significant  effects  on 
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time-dependent  fracture  of  elastomers.  Therefore,  several 
experiments  along  these  lines  are  planned. 

The  polyurethanes  being  used  in  this  study  are  of  variable 
stiffness.  The  possibility  of  mechanical  relaxation  cannot  be 
ignored.  The  dilatometric  measurement  will  be  repeated  in  a 
standard  dilatometer  using  an  appropriate  silicone  fluid  as  the 
confining  liquid.  If  the  same  transition  is  observed,  mechanical 
relaxation  will  be  eliminated  as  a  possibility. 

The  presence  or  absence  of  a  crystalline  phase  will  be 
determined.  The  polybutadiene  microstructure  of  the  R45-HT 
consists  of  60%  1 ,4-trans.  20%  1 ,4-cis  and  20%  vinyl  units.  It 
is  possible  that  the  1 ,4-trans  units  could  crystallize  at  the  low 
temperatures  used  in  our  experiments.  The  melting  temperatures 
reported8  for  1 ,4-trans  (99-100%)  polybutadiene  are  97  and  145°C 
for  Modifications  I  and  II,  respectively.  The  inflection  shown 
in  Figure  22  at  about  110°C  is  in  the  right  temperature  range 
for  crystallized  1 ,4-trans  units.  X-ray  studies  are  planned  to 
determine  if  the  sample  is  at  all  crystalline.  If  it  is,  the 
unit  cell  dimensions  will  be  compared  with  those  reported  for 
1,4-trans  polybutadiene. 

Other  studies  will  be  carried  out  as  required. 
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FIGURE  1.  Vacuum  Mixing  and  Casting  Apparatus  for  Epoxy  Samples 

FIGURE  2.  Four-Cavity  Mold  for  Epoxy  Compact  Tension  Specimen 
Casting 

FIGURE  3.  Mooney-Rivlin  Plot  of  201  DMSO  Pre-swollen  Epoxy 
Sample 

FIGURE  4.  Dilatometric  Determination  of  Glass  Transition 
Temperature  of  Pre-swollen  DMF-Epoxy  Sample 

FIGURE  5.  Critical  Strain  Energy  Density  versus  Reciprocal 
Crack  Length 

FIGURE  6.  Crack  Insertion  Tool  for  Fracture  Specimens 

FIGURE  7.  Representative  Computer  Plot  of  Dynamic  Storage 
Modulus  of  Epoxy  Material 

FIGURE  8.  Computer  Generated  Storage  Modulus  Master  Curve  for 
Epoxy  Material 

FIGURE  9.  Electron  Micrographs  of  Styragel  and  Merrifield 
Resin  in  Polystyrene  Matrix 

Figure  10.  Electron  Micrographs  of  Commercial  Gelled  Resin 
Particles 

FIGURE  11.  Mooney-Rivlin  Plot  of  Samples  TB i ,  TB2,  TS i  and  PBi 

FIGURE  12.  Mooney-Rivlin  Plot  of  Sample  TB2  ,  TS i  and  PB> 
in  Hexane. 
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FIGURE  13. 


True  Stress  vs.  Strain  of  Samples 


FIGURE 

FIGURE 

FIGURE 

FIGURE 

FIGURE 

FIGURE 

FIGURE 

FIGURE 

FIGURE 


14.  Stress-Relaxation  Data  for  Sample  TBj 

15.  Loss  Modulus  vs.  Statis  Deformation  e  for  Sample  TBi 

16.  Tan  6  vs.  Statis  Deformation  e  for  Sample  Bi 

17.  Temperature  and  Rate  Effects  on  Tearing  Energy  T 
for  Sample  TBi i 

18.  Master  Curves  of  Tear  Energy  vs  Effective  Strain 
Rate  for  Samples  TBii,  TB23  and  TB33 

19.  Schematic  Illustrations  of  Thermal  Expansion 
Apparatus 

20.  Typical  Differential  Thermal  Expansion  Curve  for 
His  cis  Polybutadiene  Rubber7 

21.  Thermal  Expansion  Curves  for  Diene  35NFA  and  Pale 
Crepe  Obtained  on  the  Apparatus  Described  in 
Figure  19 

22.  Typical  Thermal  Expansion  Curve  Obtained  using  Sample 
TB2„  and  the  Apparatus  Described  in  Figure  19. 
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Figure  7.  Representative  Computer  Plot  of  Dynamic  Storage  Modulus 
of  Epoxy  Material. 


i  iii.  OMi  i  irmo/'in  i 

Figure  8.  Computer  Generated  Storage  Modulus  Master  Curve  for 
Epoxy  Material . 


MERRIFIELD  RESIN 


Figure  11.  Mooney-Rivl in  plot  of  samples:  V  :  TBi,  0  :  TB2 ,  O  :  TS]  and 
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Figure  17.  Temperature  and  rate  effects  on  tearing  energy  T 
for  sample  TB  i  i . 


Figure  18.  Master  curves  of  tear  energy  vs.  effective  strain 
rate  for  samples  V  TBn,  O  :  TB2  3  and  □  :TB3 


EXPANSION 


TCMPCftATUftC  «N  -M  V  ON  l*ON 'CONSTANTA  N  THERMOCOUPLE 


Figure  20.  Typical  differential  thermal 
expansion  curve  for  high  cis 
polybutadiene  rubber.7 


(  0  C  ) 


Tampera  Jure 

Figure  21.  Thermal  expansion  curves  for  Diene  35NFA 
and  Pale  Crepe  obtained  on  the  apparatus 
described  in  Figure  19. 


described  in  Figure  19. 
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MICROGRAPHS 

Plate  1  Nodules  ranging  in  size  from  100-400  nm.  from  a 
4/3  TETA/EPON  828  surface  after  etching  24  hours 
with  1M  Cr03. 

Plate  2  Nodules  of  100  to  200  nm.  in  size  from  a  4/3 
TETA/EPON  828  impact  fracture  surface. 

Plate  3  An  SEM  micrograph  of  a  crack  initiation  region 
from  a  1/1  MDA/EPON  828  fracture  surface.  Note 
the  fibrillar  structures. 

Plate  4  Nodules  taken  from  surface  seen  in  Plate  3  after 

replication,  near  a  fibril. 

Plate  5  Replica  of  surface  seen  in  Plate  3.  The  end  of  a 

fibril  which  was  torn  from  the  surface.  Note  the 
directional  flow  of  material  towards  the  ridge. 

Plate  6  Replica  of  surface  seen  in  Plate  3.  A  fibril 
torn  from  the  surface  overshadowing  a  ridge 
behind  it. 

Plate  7  Replica  of  surface  seen  in  Plate  3.  A  fibril  above 
the  ridge  from  which  it  emanated. 


oeitlcs  r;ui;itt;>  in  size  Irvi  1 '  I  '>  -  4  Cl  'I  i,n  froi;i  a  i/.'i  'ITTf/iniN 
i  rfiitr  ;irii  i  i  liliii*  liners  villi  1'  t * 3  . 


of  a  crack  initiation  region  from  ;i  l/J  MI)A/ 
surface .  Note  the  fibrillar  structures. 
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